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NOTICES 


“When U.8. Govemment drawings, specifications, or other data are used 
for any purpose other than a definitely related government procurement operation, 
the Government thereby incurs no responsibility nor any obligation whatsoever, 
and the fact that the Government may have formulated, furnished, or in any way 
supplied the said drawings, specifications or other cata, {s not to be regarded by 
implication or otherwise, or in any manner licensing the t older or any other person 
Or corporation, or conveying any rights or permission to manufacture use, or scl] 
any patented invention that may in any way be related thereto. ": 


FOREWORD 


This report as submitted by Ultrasystems, Inc.,, Environmental and Applied 
Sctences Division, 2400 Michelson Drive, Irvine, California 92664, under Contract 
No. F04611-73-C-0638, Job Order No. 305909LZ with the Air Force Rocket Propul- 
sion Laboratory, Edwards, CA 93523, 


This report consists of three volumes. Volume I describes a computer 
procram for the prediction of Solid Propellant Rocket Motor Performance. Tne com- 
puter program described herein will be referred to as the SPP program, and describes 
the engineering analysis which was used in developing this computer program and 
the results obtained to date, 


Volume II of this repoit is a programming document of the computer program 
which was developed under this contract. It includes a subroutine-by-subroutine 
description of all of the elements of the SPP program, 


Volume III of this report {8 a Program User's Maaual which describes the 
input necessary to execute the SPP computer program and the information required 
to interpret the output. A sample case {8s also included. 


This report has been reviewed by the Information Office/DOZ and ts re- 
leasable to the National Technical Information Service (NTIS). At NTIS it wil) be 
available to the general public, {ncluding foreign nations. 


This report {s unclassified and suitable for public release, 
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computed as a series of independent efficiencies. The program currently treats 
the following losses: two-dimensional/two-phase (coupled), nozzle erosion, 
kinetics, boundary layer, combustion efficiency, submergence. The program 
predicts average uelivered pertormance, as well as mass flow, pressure, thrust, 
impulse, and specific impulse as functions of time and trajectory, 


In order to assess the validity of the SPP program, calculated results were com- 
pared to firing data for four different types of motors, ‘Vhile conclusive state- 
ments regarding the accuracy and range of validity of the SPP program cannot be 
made until additional verification efforts are conducted, the results of these four 
test cases were encouraging. These calculations also served to demonstrate 

the desirability of eliminating some of the present limitations of the program. 
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This ts Volume III of a three-part manual which describes a computer pro- 
gram for the prediction of Solid Propellant Rocket Motor Performance, The com- 
puter program described herein will be referred to as the SPP program, 
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Volume I of this report describes the engineering analysis which was used 
in developing this computer program. 


Volume II of this report is a programming document of the computer program 
which was developed under this contract. It includes a subroutine-by-subroutine 
description of all of the elements of the SPP program. 


Volume ITI of this report is a Program User's Manual which describes the 
input necessary to 2xecute the SPP computer program and the information required 
to interpret the outnut. A sample case is also included, 
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1, METHOD OF APPROACH 


The goal of the present effort was to develop a computerized analysis tech- 
% nique that could predict delivered specific impulse to within +2%, and delivered 
thrust coefficient and total impulse to within +5%, The method of approach used 
in the development of the Solid Rocket Performance Program (SPP) stressed modu- 
larity and flexibility, The program was designed so that the various elements 
which comprise the total performance prediction methodology could be utilized as 
an integrated set, inaividually, or in any user determined combination. In most 
cases both analytice! and empirical metods have heen incorporated into the pro- 
gram to accourt for the phenomena which causes motor performance to deviate 
from a specified “ideal.” The user, depending on his needs, can select either method. The | 


option ts »henomena oriented, so some effects may be treated empirically (simple, | 
short computation time) while others can be evaluated analytically (more calcula- } 
tion time, greater accurecy and generality). | 


Tne program has been structured in a modular fashion to allow future im- 
provements and mod.ficatlons to be easily incorporated, The program modules 
(with the exception of the master contro] module) are pre-existing computer 
programs which have ‘een adapted for the present purpose. The following pre- 
existing programs were used as the basis for the SPP program: 

One Dimensional Equilibrium Program (opz)' 1) 

One Dimensional Kinetics Frogrem (opk)! ° 

| Two Dimensional Two Phase Program (TD2P)' : 

Turbulent Boungary Layer Program (rer)! 4) 

Three Dimensional Grain Design Program (64101 yf s) 

Metor Ballistics Program (637 y! 6 


These programs have beer, combined w:th a contro! inoduis which permits them to 
be run teqether automatically. The extensive use of interne! communication be- 
tween modulcs has just about eliminated the need for redundent inputs by the 
user, 


The next two subsections wij] define the “!nsses" which have been con- 
sidered and the mann2or in which delivered performance is calculated. The balance 
of the section briefly lescribes ihe various modules, what they do, how they do 
it, and their tnter-rel.tionsnips,. 


ee Terminoiogy and Definition of Losses 


The basic method of approach used in the develonn.ent of the SPP computer 
program is to calculate the maximum or ideal ‘theoretical’ performance of a rocket 
motor uid to subtract from that maximum performance the individual losses or de- 
viations from the “theoretical” perfonnaance which are known to occur in a real 
rocket motor. This ncchod of approach is neither new to the field of rocket per- 
formance prediction or engineering analysis ir. general, It is, however, founded 
on our understanding of the basic physical processes involved in a real rocket 
motor and the interrelation of these processes with each other, 


The basic assumption involved in this type of approach is that certain "losses ” 
can be treated as essentially independent from the other losses which occur in 
the rocket motor, As noted in Section 1.2, all of the losses should not be treated 
independently. There ure some interactions strong enough to warrant considera- 
tion. With these ex-eptions in mind, the concept of treating the losses indepen- 
der.tly appears to be valid, considering the success of this approach and others 
of a similar nature. 


The term “Joss" as used here is really a misnomer, since it is used to 
describe a known physical mechanism or phenomena which causes the performance 
of wu "real" solid rocket motor to be Jess than the so-called "theoretical" perfor- 
mance, In the same sense the terminology "theoretical," or "ideal," is alsoa 
misnomer, because the wkcle methodology incorporated herein is based on our 
theoretical understan-iing of the processes which are prec: in real rocket motors, 
However, this terminclogy has gained wide spread acceptance in the engineering 
community and wil] be used in this report. 


ln, order to further define the terminology used here the following defini- 
tions will be used: 


@ Performance = unless atherwise specified, performance 
will bo in terms of specific impulse, I. ; 
delivered to a vacuum, P 

@Thacretica) Ian - the maximum possible I... which can be de- 


livered at the initial area ratio by a propel- 
lant at a given chamber pressure and total 
enthalpy. Referred to as Isp, ° 


@ Loss = the decrement from the theoretical I,,, which 
can be attributed to a physical phenomens 
net included in the calculation of Isp° 


@ Two phase tiow loss - the decrement in performance due to finite 
velocity and temperature differences be- 
tween the gas and condensed phase. 


@Two dimensional or - the decrement in performance which is due 
cs divergence lcss to the momentum of the rocket exhaust not 
being totally aligned with the axis of the 
motor, 
| 
@Finite Rate Kinetics - the decrement in performance which is due 
loss to incomplete transfer of latent heat to sen- 


sible heat caused by the finite time required 
for gas phase chemical reactions to occur, 


@Boundary Layer Loss - the decrement in performance due to viscous 
forces adjacent to the nozzle wall and heat 
transfer to the nozzle wall. ° 


as used here, the term combustion efficiency 
will refer only to the degradation of I,, due 
to a departure of chamber total temperature 
from the theoretically calculated total temper- 
ature. 


the decrement in performance due to propel- 
lant grain extending pass the nozzle inlet. 


@Erosion Loss > the decrement in performance due to the ero- 
sion of the nozzle throat. This loss encom- 
passes both area ratio effects (considcred 
only in the TD2P module) and throat roughness 
effects (not considered in the present analysis) 


@Combustion Efficiency 


@Submergency Loss 


eatin 


@Delivered Performance - the calculated performance of the rocket 
motor which includes all of the losses con- 
sidered herein. 

With our terminology now defined, we can now dascribe the rational and 
methods used by the SPP computer program to calcuia:. the delivered performance 


of a solid rocket motor, 


1.1.1 Calculation of Delivered Performance 


Within the framework already discussed (i.e., separation of individual 
losses), there are tn-ee basic methods which can be used to calculate delivered 
Tene The first method is to calculate discrete decrements in performance, AI sp’ 
and to subtract thee? losses from the theoretical specific impulsoa, I. 


p @ 
That is th 
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The second method is to calculate a series of efficiencies, n's, to be applied 
multiplicatively, i.e. 


EAP: - ETP 


‘spy "sPep TM oss ° (1-2) 
: Both of the above methods are the same through first order, thus, if the losses 
are small they yleld the same result. 


The third method was selected by the JANNAF Committee on Performance 
Standardization for liguid rocket engines and is a combination of the first two 
methods. That is, 21] of the losses except the boundary layer loss are treated as 
efficiences while the boundary layer loss is treated ss a decrement. Hence, 

‘Son ‘sPan Toss ~ Msper (1-3) 
This latter method >f apy.roach was selected for incorporation in the SPP code for 
two reasons. The first reason being that there was no compelling reason to de- 
viate from what has been established as a “standardized” method In the Industry, 
and secondly, that any other combination of the above methods yields substantially 
the same end result. 


The main p-cnlem associated with the separation technique is, of course, 
the selection of a Fase or reference value of perfcrmance to compare against the a 
model which incorporates the loss. Ideally, this reference calculation should be 
nade by the same computer code which is selected to make the loss calculation 
since this would minimize deviations Jue to numerical techniques, truncation and 
roundoff errors, However, as is generally the case, su~ an approach is not 
practical. For examnle, both the equations for the two $2 flow and chemical 
hinetics losses become “stiff” in the equilibrium limit result in unacceptably 
long and difficult caiculational procedures, Even if sophisticated numerical tect - 
niques are used to «alculote the “near” equilibrium case, extrapolation to the ful: 
equilibrium solutior, introduces errors which negate any advantage gained by using 
tho same computer code. Hence, the following criteria was adopted in ~siculating 
the reference performance for each of the losses (except tie boundary layer loss 
which is calculated directly). : 


The reference method should be the most exact method of cal- 
culation availoble. The conceptual differences between the 
reference anc ‘oss modcis should include only those mechan- 
isms which cause the loss. 
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TABLE I-! INTERACTION OF PHYSICAL PHENOMENA WITH 
PERFORMANCE LOSS CALCULATIONS 


PERFORMANCE {Divergence 
LOSSES Loss 


PHENOMENA 


One-Dimensional] 
Flow ne 


‘Jiscous And 
Heat Transfer 


Finite Rate 
Chemistry 


Multiphase Flow 


Incomplete 
Co.nbustion 


Logend: 
1, Primary Importance (could be + 0.2% effect on ign) 
2. Secondary Importance (probaLiy <0.2'4 effect on Lo! 
3. Generally Not :mportant 


Table 1-1 (from Ref. 9), shows the interaction of the physical phenomena 
considered here with each of the performance losses. The table should be read 
down and then to the left. For example, the main interaction between the diver- 
gence loss and the vavious physical phenomena Is that assoclated with multi- 
phased flow, while the boundary layer loss is chlefly affected by the nature of 
the specified boundary layer edge conditions, (i.e., from a one or two-dimensional 
flow field calculation), Table 1-2 is a condensation of Table 1-1, and indicates 
only those cross coupling effects that were deemed strong enough to warrant con- 
sideration. 


Table 1-2 Cross Coupling of Performance Losses 


Loss Mechanism Cross Coupling 


Finite Rate Chemistry (Kinetics) none 
Divergence Loss Two Phase Flow 
Twe Phase Flow none 
Boundary Layer Loss Non One-Dimensional Flow 
Combustion Eficlency Loss none 


Table 1-2 shows that while the kinetics, two phase flow and combustion 
efficiency losses can be treated independently, the divergence and boundary layer 
losses are coupled to non-one-dimensional and two phase f.ow phenomena, re- 
spectively. Hence, it was decided that the best course of action to take in the 
SPP code was to couple the two phase flow and divergence losses together and to 
u 23 the results of that calculation to obtain the edge conditions for the boundary 
layer calculation. 


With the above in mind, the following sub-sections describe the basic modules 
and assumptions used '!n the SPP code. 


1.2 Master Control Module (MC) 


The master control module (MC) controls the execution of the SPP computer 
code. The MC module selects, via user input, which of the five basic calculation 
modules (SDE, BAL, @i.K, TD2P, and TBL) are to be executed, and whether calcu- 
lated, input, or empirical, losses are to be selected for the calculation of delivered 
specific impulse, thrust, and total impulse. 
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This module handles or specifies all of the linkace communications bet- 
ween modules. The MC module selects which parameters (either calculated or 
input) are transmitted from one module tc another and tries to determine the best 
overall delivered performance consistant with the input to the SPP computer pro- 
gram, 


Figure 1-1 shows the basic input and output of the MC module, while Table 
1-3 shows the variables used to calculate delivered thrust, Tgp! and total impulse 
and indicates the module in which each is calculated, 


Table 1-3 Variables Used for Performance 
Prediction by the MC Module 


@DE/One Dimensional Equilibrium ee | va 
Ptr §PREQ 
BAL/Grain Desig and Ballistics ro. m,t*, Lt, y es 85 functions of t 


@DK/One Dimensional Kinetics loss a 
KIN 


TD2P/Two Dimensiona) Two Phase flow loss ae "TD2P as functions of g! and C 


TBL/Turbulent Boundary Layer loss al. 
Pre 


The vacuum delivered specific impulse is calculated as in equation !-3 
and is written In full nere, 


= J 


aoa Spy, "KIN "rbap © "cr ° "sup M (1-4) 


SPrpr 


The theoretical specific impulse, L, ° #8 base’ on the initlal nozzle 
“7 +t, 


exit area ratio. If a motor experiences throat erosion duriag the course of a firing 
the nozzle exit area ratio, and hence, tc specific impuise, varies with timc, 
This gives rise to an additional loss mechanism; one which was only mentioned 
but not discussed in Section 1.1. The serformance decrement assoctated with this 
phenomena is usual.y called the erosior 'oss. While this loss can be estimeted 
based on one-dimensional calculations, it was felt that a two-dimensional esti- 
mate of this effect weuld be superior. Thus. for eroding nozzles, Nrp2p in equa- 
tlon (1-4) is modified as follows: 

I (e) 


oF 
: Ther 
Mepep = Mop Ge (1-5) 
Prnap 
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*t 

where 

_ _ itty (¢ (t) dt 

1 Tp2P i a 
re) 


When defined in this way nypop includes the "erosion" loss in addition to 
the coupled two-phase and divergence losses. The erosion loss, as defined here- 
in, is equal to the ratio of the 2-D, 2-phase on at the time averaged expansion 


ratlo, ys (<), ‘o the specific impulse evaluated at the initial expansion ratio, 
TD2P 


ae (e " ). (Future modifications to the program should break the erosion loss 


out as a separate efficiency.) 


The combustion efficiency, CE is defined as follows: 


a cr 1.0, if the previous relation yields a value 

greater than unity. 
This definition prevents the mass flow effect relate! to the nozzle discharge cocffiviont, 
Cn, from being counted twice since it is an integral p«rt of the 2D-2 phase ioss, In 
the empirical performance prediction procedure "CE is taken to be equal to Nowe The 
emnirical determination of c* efficiency is discussed in Volume i of this report, 


The loss due to finite rate chemical! kinetics is treated as an efficiency 
factor, nen. The kirette efficiency is calculates as the .atia of an Isp calcula- 
ted by the ODK module, to an I Sp calculated by a specie) option in the ODE module. 
This special ODE option, referred to as the “Kestrictea Equilibrium" option causes 
the ODE program to compute an ay based on the more raostrictive physical assump- 
tions employed in the ODK calculation, This procedure is iecessary if certuin 
losses arc not to be counted twice. The aforementioned assumptions are descrip- 
ed in Section 1,5. 


The submergence efficiency, Nour i rased on an empirical correlation. enc 
is a function of the throat radius, nass fraction of condensed phase, length of svb- 
mergence to length of the {nterna] motor, and chamber rressure. This correlation 


is discussec in more detall !n Volume I of this report. 


The decrement in performance due to the formation of a turbulent boundary 
layer in the motor nozzle is transmitted directly to the MC module and is reported 
as a function of area ratio. The decrement in 5p is taken io be the value corres- 
ponding to the last (i.e. the maximum) area ratio input by the user. This decrement 
is assumed to be independent of erosion induced area ratio changes (it is more 
a function of the wettad length of the nozzle) and throat roughness. 


The equations used to calculate the "theoretical" performance losses, and 
the formulas used in the "empirical", or simplified, performance calculations are 
presented in Volume I and Volume II. 


1.3 Theoretical i 3 Module (@DE) 


The theoretical I5p module (ODE) is an outgrowth of the NASA Lewis chemical 
equilibrium computer program described in NASA SP-273 (Reference 1) and in the 
JANNAF standard TDK computer program manual (Reference 2). 


This module caiculates the theoretical, or maximum, I, which can be de- 
livered by a propellant at a given area ratio, The ODE module uses a free energy 
minimization technique along with the following assumptions to calculate the 
vacuum I sp° 

One dimensionz! flow 

Thermal and velocity equillbrium of the condensed phase with the gas 

phase flow 

This version of the ODE program has also been modified to calculate the 
transport properties nceded by the ballistics (BAL) module, two phase flow loss 
(TD2P) module, and the boundary layer loss (TBL) module. 


In addition, the program has been modified to calculate the reference con- 
dition for the loss in performance due to iinite rate chemical kinetics. This option 
is referred to as the “restricted cquillbrium” option and ls discussed in more de- 
tail in section 1.5. 


Figure 1-2 illustrates the major functions of the ODE module as well as 
the inputs and outputs of the module. 


1.4 Grain Design and Ballistics Module 
When utilized as part of the overall performance prediction methodology 
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the purpose of the Grain Design and Ballistics (GD&BC) module is to provide values 
of pressure, mass flow rate, throat radius, L*, and C* efficiency, as functions of 
time, Time averaged values of these quantities are also computed and used in the 
other modules. Figuwe 1-3 illustrates how the GD&BC module functions in this 
mode of operation. Various thermodynamic and transport properties are internally 1 
transmitted to the GD&BC module from the ODE module. Lixewise, the nozzle 

geometry and ambient oressure history are obtained from the Master Control module. 

The final thrust versus time output of the programs in this integrated mode of oper- 

ation are based on the GD&BC calculated chamber pressure and mass flow rate; the 

delivered I sp is calculated from the combined results of the other modules. 


The GD&BC module can also be employed in a stand-alone mode, In this 
case, all of the parameters which would otherwise be supplied by the ODE and 
MC modules must be input directly into the Ballistics module. The thrust versus 
time curve in this mode of operation is based on the calculated mass flow rate, 
an empirical C* efficiercy, and an input value of nozzle (C,) efficiency. 


The essential features of the GD&BC module are: 


The Hercules lnc. 6410] grain geometry analysis program, omitting non- : 
essential elen.ents for considerable size reduction, and modified to be com- 

patible with the Lockheed Propulsion Company interior ballistics analysis 

program, ‘ 


An improved version of the Lockheed Fropulsion Comnany 637 interior bellis- 
tics analysis program, 


A correlation of C* effictency,. 
A correlation of nozzle throat erosica 


1.4.1 Grain Geomet-v Program 


The Grain Design Program examincs a grain in three dimensions and solves 
for the surface versus bum distance relationship. The grains are not limited to 
families of grains such as Star, slotted tube, cic., but will handle these and more 
complex graizs, A detailed users manual for the original Heicules 64101 program 
exists as Rof, (5), The following extracts are essential portions for present purposes. 


1.4.1.1 Method of Approach 


The Grain Design Program simulates crafting techniques used in develop- 
ing the surface area versus burn distance curve of a solid propellant motor, given 
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an arbitrary number of symmetrical angular sectors. The program uses a right- 
handed orthogonal :oordirate system in X, Y, Z where X is the longitudinal axis 
passing along the center lire of the motor and the Y and Z axes pass from the center 
of the motor to the outer boundary. 


Three basic input figures ar. used in various combinations to describe the 
initial propellant geometry. Thesr ‘igures are: frustrum of a right circular cone, 
right triargular prism, and a sphere, A union of the input figures is performed by 
computing the intersection of a line Y=C and each of the inpuc figures, The over- 
lapping points are vliminated and those left represent the propellant geometry at 
the particular X, Y and burn coordinate consistent with the symmetry of the pro- 
pellant. 


Each figure méy be input as consisting of void or »ropellant, where the 
void is burning outward and the propellant is buming inward. Since the burning 
is assumed normal to the surface and the surface area is equal to the derivative of 
the volume of the p opellant with respect to the normal, the surface is computed 
as the derivative of the volume of the propellant with respect to the bum c‘'stance, 
Volume can be comruted as the integral of the propellant areas at each X station 
with respect to X, and propellant area can be expressed as the integral of the 
lengths across the propellant with respect to Y. 


1.4.1.2 General Description of the Solution 
The besic mass balance equation 
Messrs (1-7) 


is tae basis for the solution. By eliminating time and density, equation (1-7) 
bucomes 
s=dv/db - v/xb (1-8) 


ard from the mean value theorem s will occur semewhere in the interval .b. The 
half-way point was a:bitvarily accepted. 

m= mass 

r = rate of bum 

0 = mass deasity 

s = burning surface area 

v = volume of propellant 

b = bum distance 


In the Solution for b versus S the procvain .>-s not make any direct analytic 


solution for a burning surface, but rather makes v.:! 1 :tions for the voiume of 
| propellant at each burn distance specified in the inp ., :-cords the change in volume 
| between two consecutive bum distances, and cdiv::. . »y the difference in the two 

4 burn distances. Thus, the user should keep in mir’ shot he is preparing input 


to simulate the chcnges in volume which wil] occi:: ii the motor, not to airectly 
compute buming s'inace areas. 


1.4.1.3 Method oz Determining Volumes: Confic = 4 Simulation 


The case is assumed full of propellant init::.°:. The voids in the case 
are simulated by figures - cones, cylinders, privir: . ~ spheres - which may over- 
lap and/or protrude outside the case,if needed. % = should simulate the void 


as it is inittally, Tie figures may be norma] ot’: 0. :sures, grain-fill inbum- 
ing figures, or nonowning figures and may have «.. © corners and eages,if neciad, 
The order of input can be important. Each figure 2+. > placed in any orientation 
anywhere in space, inside or outside of the 5:7: ‘ «oof these is tesrriped 
in detail later, Eac’ figure is located in throes ci: ral space by the follow.ne 
means. 
‘ . 
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Spheres: 
oa aes Fs | 


r 


Generally, the grain will be symmetric about the axis of the motor. If 
there are symmetries in the grain they should be reflected in the input in order to 


Save work and computer time. See Figure 1-4, The symmetries are automatically 
accounted for in the whole motor, 


Figure 1-4: 


Sample symmetries: (End views of motors) 


Six 
Eieht 
Z vA Z 
= Y 9 > 
SX = 6.0 SX = 3.0 SX = 2.0 
6 = 60° @= 120° Q= 180° 


(Tre X axis is perpendicular to the peper) 
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1.4.2 Interior Ba. listics 


The Interior Ballistics Program models the internal ballistics of solid rocket 
motors, including motor environmental effects on burning rate (erosion, radiation), 
effects of radial slots between segments on the port flow port flow’ nressure drop, 
and delayed burning (ignition) in deep narrow slots. Semi-empirical correlations 
for combustion ineff'clency and nozzle throat erosion also are included. This 
digital computer program was adapted from an existing model used to analyze nozzle- 
less solid rocket motors 6 - It was selected because it contains features lacking 
in many published .ntemal ballistics programs. The intent was to predict internal 
ballistics for motors ranging from small tactical motors, to lerge boosters, wherein 
aspects of all the atocementioned features are encountered, The major effort in 
the present program development involved updating certain model expressions, aid 
rendering the computational procedure compatible with the Grain Design 
Module and the balance of the SPP Program Modules. 
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1.5 Finite Rate chemical Kinetics Loss Module (@DK) 


In most rocket rotors there is a loss of performance due to the fact that 
a finite time is required for chemical reactions to procecd to the maximum amount 
of sensible heat release, Normally, this loss is caus: 3 by the incomplete recom- 
hination in the rocket nozzle of the highly energctic srecies which are formed curing 
.ombustion in the rocket motor. 


The @DK module only considers the losses whick are due to gas phase 
chemical kinetics, Other losses which include incomplete or partial solidification 
of condensed phases 3re ignored in the ODK module since adequate models of these 
processes have yet to be formulated, The ODK modulo wlso assumes that the flow 
is one dimensional, with the gas and condensed phases in dynamic and thermal 
equilibrium with each other. No transfer of mass between “18 phase species and 
the condensed phase species are allowed. The liquid phase, however, is allowec 
to undergo solidification when the gas temcrature reaches the melting temperaturc 
of the condensed phese. 


In order to properly calcusate the chemical kir.ctic performance efficiency, 
Tone 2 reference equilibrium 55 which is baseca on the seme assumptions as ODK 


should be used. In order to provide this reference value; the @DE module was modi- 
fied to ca'culate a "restricted equilibrium” specific impulse ‘spar’ This modifi- 


cation to the ODE mudule allows the available gas phase species to remain in 
chemical equilibriun with themselves while restricting any mass transfer to the 
condensed phases. dence, the kinetic performance efficiency, which Is defined as 


"KIN 5 sbODK Asp. 


would approe“h unity as the chemical kinetic reaction rate constants approach in- 


Another advantage of the use of the restricted equilibrium option for computing 
kinetic efficiencies is that only thermodynamically important species need be con- 
sidered in the ODK-2* calculations. Thus, trace species, and reactions involv- 
ing them, need not be included in this calculation. 


Figure 1-5 siiows the basic input, output and assumptions used in the ODK 
module. A comprehensive set of reactions for aluminized propellants is also in- 
cluded in the volume, (See Table 2-13) Appendix A of Volume I indicates how the 
reactions and species whe:-< selected for aluminized propellants. 


— 


1,6 __ Two Dimensional-Two Phase Loss Module (TD2P) 

Almost all rocket nozzles of practical interest have a !ois in performance 
due to a nonaxial component of velocity at the nozzie exit plane. This loss ts 
usuélly referred to as & divergence, or two dimensiona! flow, loss. If the pro- 
pellant (s metalized, there is an additional loss when the particles cannot main- 
tain a state of dynamic and thermal equilibrium with the gas flow. This loss 
is referred to as the two phase flow loss. 


There are two 4'stinct -oupling effects between the two diaensionel and 
two phase flow losses. The first is that if the gas and conucnsed phase are not 
in dynamic (velocity) equilibrium with each other then the gas and particle stream- 
lines are not coincident . This can have a noticeable effect on the divergence 
loss. The second coupling effect is on the mass flux through the nozzle. It is 
now well known that the mass flow through a nozzle for single phase two dimen- 
sional flow is less than the omount of mass flow that can be achieved at a given 
chamber pressure for an equivalent one dimensional nozzle fiow. This effect results in 
nozzle discharge coetficients being less thin unity. However, it is also wel! 
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known that in nonequilibrium two phase flow more mass flows through a nozzle 
at a given chamber pressure than would if the particles were in velocity and 
thermal equilibrium. By itself,this effect would produce Ch values greater than 
unity. For most metalized solid propellant motors the net result of these two op- 
posing effects is 4 discharge coefficient greater than unity. If the combustion 
efficiency were 160% (complete combustion) the C* effictency of the motor, n 
is given by 


foi! 


Nex™ [om , 


Thus, in most cases, Noxwould be less than unity for metalized propellants, and 
would always be greater than unity for non-metalized propellants*. As a result of 
the relationship between N, end Ch, the Ch effect must be corrected out of the 
empirically obtained c* efficiency {f one is to obtain a value for combustion 
efficiency whic’ does not count this factor twice. 


The TD2P mudule is a modified version of the computer program developed 
by Kliegel and Nickerson (Ref. 3). This module, which has been extensively 
modernized, calculates the losses which are usually the largest encountered in 
most solid rocket motors, These losses are, the two phase flow loss, erosion 
loss, and divergence loss, 


The basic invuts, charactcristics, and outputs of the TD2P module are 
shown in Figure 1-€, Since this is the most important loss module in the SPP 
code it is felt that the following comments are in order, 


*Most solid rocket motors have normalized radi! of curvature ~2 and, in the ab- 
sence of particles C,, tends to be .99<Cp <1, or, in other words, with perfect 
combustion n, woul bounded as follows, 1 ents <1.01. This small effect 
tends to be masked by measurement and data reduction efrors, and many times 
test engineers are loathe to report values of n* greater than unity due to the 
mistaken asswnption that such values are physically impossible to achieve, 
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1) The TD2P module assumes that the gas phase flow behaves as an ideal 
gas. The gas phase properties are calculated by matching ideal one di- 
mensional “equivalent” gas results to the output of the ODE module at 
the solidification point. In order that the condensed phase properties, 
and, hence, the results of the TD2P module, would not be artificially de- 
pendent on a:ea ratio they were taken as the actual values obtained from 
the @DE module. | 


2) The transport properties which are transmitted from the DE module 
are sclected to insure the greatest degree of accuracy in the throat region. 


3) The reference Ig, (one dimensional flow with zero lags, Igy ), to 


which the TD2P Isp is ratioed, is calculated using the same iddal Gas properties 
as are used in the TD2P analysis. By using a ratio calculated in this man- 

ner, instead of using the absolute I... number from TD2P, the effect of 
neglecting real gas effects is minimized, Also, extra care was taken to 
accurately caiculate the reference I... during solidification so as to avoid 

the use of difterent assumptions in fe two modules. 


4) The fact that the two phase flow loss tends to decrease with increased 
throat diameter (erosion) was intentionally ignored, since it was felt that 
the added loss due to throat roughness would more than compensate for this 
effect. 


5) The transonic analysis currently incorporated int> the TD2P module is 
restricted to running normalized throat radii of curvature of no less than 
approximately !.5 and inlet angles less than about 45°. In addition the 
program requires the throat geometry to be circular with equal upstream 
and downstream radii of curvature, These restrictions on the geometries 
which may be realistically treated are a major shortcoming of the TD2° 
module, 


6) The edge conditions which are generated by the TD2P module for the 
turbulent bourdary layer, TBL, module have been adjusted so that the sub- 
sonic and transonic regions will faired into the supersonic regions Ina 
smooth manner. The stagnation conditions transnmuited to the TBL module 
have also been adjusted so as to allow the TBL module to reproduce the 
TD2P edge conditions (velocity, pressure) as closely as possible. 


Lad Turbulen cary Layer Module (TBL 


The losses associated with the presence of a boundary layer result from 
a transfer of heet to tie nozzle walls and irom a reduction In momentum due to shear 
stress along the walls. These losses can be calculated by properly integrating the 
stress tensor along the wall (pressure plus skin friction), or by calculating the 
momentum deficit in the boundary layer at the nozzle exit plane. The second 
method is the easiest to apply, and is the JANNAF recommended procedure. 


The boundary layer development in the nozzle is currently cemputed using 
the Bartz-Silver integial method, as embodied in the JANNAF adcpted TBL 
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program . This is not the most accurate houndary layer program available, 
however, its virtues (speed and simplicity), coupled with the usually small 
magnitude of the loss (usually less tha 1% for all but small mctors), make i* a 
reasonable choice for the present purpose. 


The boundary layer solution requires a specification of the inviscid flow 
along the wall streamline. In the SPP program these edge conditions are supplied 
automatically by the TD2P module (if executed). This allows the interaction be- 
tween two dimensional and boundary layer effects to be properly accounted for. 

The TD2P module supplies the TBL module with the nozzle geometry, Mach number 
along the wall streamline, ideal gas properties, and stagnation pressure and temp- 
erature. This procedure allows the reinaining edge properties to be computed in- 
temally in tne TBL module (pressure, velocity, etc.) to match, as closely as pos- 
sible, the respective TD2P computed counterparts. 


The boundery layer loss is computed directly as an I sp decrement, rather 
than as an efficiency, using the following formula 


Pp 

- _ b 

Ne - Qreguy Q.cOS q), f eh. 
Par, 3 ™rp2pP 


in which the various quantities are a!l evaluated at the nuzzle exit plane. It 
should be pointed ovt that this equation gives the complete boundary layer loss. 
The reduction in fluid momentum due to shear and heat losses to the wall are both 
properly accounted for. 


In most cases the effect of the invicid-viscous interaction between the TD2P 
and TBL modules is small and does not warrent the subsequent iteration which re- 
quires the re-running of TD2P module using the new invicid wall contour supplied by 
the TBL module. However, in extreine cases, the above ;rocedure may be implemented 
by using the invicid wall] contours suppiled by the TBL module as input to second 
TD2P-TBL calculatior.. However, when noundary layer losses are of such importance 
that such effects skould be considered, it is doubtful that the TBL analysis, and 
the assumptions inco:poratec within, are valid. Therefore, is is suggested that in 
these cases, that a more accurate method or analysis of the turbulent boundary layer 
loss be used. 
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That is to say, when the losses are small (;,1/2%) even a poor character- 
ization of the loss (+ 50%) is acceptable. When the losses due to a turbulent 
boundary layer are large, then secondary effects (invicid-vicid flow interactions) 
are small compared to the primary boundary layer analysis loss. Hence, if the 
boundary layer loss is significant corrections made to the TBL module met! 2d of 
calculation are of minor importance since the overall method selected is of 
doubtful accuracy’. 


The input, ‘nternal linkages, characteristics, and output of the TBL module 
are schematically portrayed in Figure 1-7. 


2. INPUT DESCRIPTION AND LINKAGE 


The SPP Proqram is structured so that the various program elements which 
comprise the total performance prediction methodology can be utilized as an integra- 
ted set, as individual program. in a stand-alone mode, or in any user specified 
combination. The individual program modules have been combined with a control 
module which provides a means for inter-module communication. This allows 
the complete methodology to be executed automatically, with a minimum of user 
prepared input. Minimum is used here in the relative sense, as the amount of 
input required to perform all the analyses is not small. 


Table 2-1 is essentially a block description and flow chart of the input 
data requirements of the SPP Program, In order to exercise the total program 
capability in a single run, one would have to consider tne input requirements of 
all of the data sets listed in Table 2-1, The notes column of the Table indicates 
which data sets can be skipped if one desires to execute only a limited number 
of the program options. The subsection in which each of the data sets is de- 
scribed is also indicated in Table 2-1. 


Table 2-1 st.ould be used as a guide when preparing Input for a given pro- 
blem. It lists the data sets in the order in which they should appear in the data 
deck, and also shows the special cards which must appear in each set (first card, 
last card, etc.) if the program is to function properly. Table 2-1 is basically 
self-explanatory when used together with the detailed innout descriptions which 
follow. The following additional notes are parenthetical. 

1. Thermodynamic data will normally be on tape or on a mass storage 


permanent file. Card Input will be required only when species not 
residing on the tape are considered, 


2, Thermodynamic data below 300°K will rarely be required for solid rocket 
motor problems. 


3. The Geometry data set is not optional if one desires the post-module 
execution summary outputs, or the final performance summary tables. 


4. Figures 2-1 and 2-2 of Volume II can be very useful in understanding 
how and in what order data is transmitted between modules. 


Teme eee 


Section 2.11, following the input data descriptions, describes the internal 
(and external, via aunch cards, perm, files, etc.) linkages between the basic 
computational moduies of the program. 
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Table 2-1 Input Data Set Description 


Data Sct Section Expleining Notes 
Date Set 
catd col, no, 
1234516786910 
Tithe cards 2.1 1,2,3 
first card TITLE 


jst ond 


Thermodynamic Dats end 2.¢ 3,4 
ficst cord 1 G 


EDC ioe) CL ( 


Thetino.lynamic Data 2.2.1 1,3 
Below 300°K 


first card ei aa 
jan cord owt crus 


Geometry card 2.3 1,3 


Trajectory card aay : 2.4 1,3 
JRA} El CTARY 
jsecond card (STRAP 
t) 
Problem card 2.8 § 


Reactants card ane | 2.6.2 ¢ 


first cord ANTS 
leat card nok blank cera required 


Omit ond Insert cards 2.6.2 & 


2.6.3 ’ 


QGDE Nemciist alll 6 


first card WN 2.6.4 
second card 7.6.8 
JeSt card vi SENG 
BAL Namelist 2.7.3 ? 
rat ca 1) Rae 
last card END 
Ballistics input dete ’ 2.762 ? 


inet cord } SL) 


Reactents cards 2.6.1 e 
{test curd NTS 
blenk card required 
Omit and Insert cards 2,6.2 r) 
aed ete 
Init Gord a 
first cof ; i aes 2.6.3 


fant e:agd 


Table 2-1 (continued) 


GDE Namelisats 


tiret card NeeipetsCidCti BN ‘ 
fecond end =pOoT] SC C~dCiCié«i AS 
jant cord SC: ees) 


GDXK option for input of 2.8.1 9 
Initial Species Concentra- 
tions 


femcerd Species SSCS 
Jas, card none roquired 
Reactions FS 2.8.2 $ 


firet card ECE TL, ae 


Inert Species Option 2.8.3 9 


wuird Body Reectton 2.8.4 9 
Rate Ratios 

fietcard IRM BODY SS 

uaced thst BARD 
@oK Nomelist 2.8.8 : 9 


lust card Kes a 
lest card S00 nT 


TO2P Namelist 2.9 10 
isct cord is 
lestcondsSOSC*«S EWN 
TEL Kamotiat 2.14 82 
{Irst card on 
lost “ad 
Notoa 
! Optiar 1}, 
2 More than one .c.4a0u, 
3 May appeer tr any orice befare the PROBLEM card, 
4 Optianend, Ua heriusty: ie date already existe, 
§ Mundatury. 
6 Requires) andy af tue GE gy tloe is eelneted, 
? Sequtted andy sl ise WAL option te sebeeteJ, 
8 wequired only ire eit A oualibrium option selected, 
9 Required only if t's “De ag tion ia selected, 
16 Ke qutied anly if tc Th7P uption ts selected, 
h Mequtes.d only if tc TAL aption le selerted, 


if roguested by irput in the SPROB namelist. previously 
qoenereted linkey  urta ray be included efter the SPAGS 
aamelist set (see Sectiun 2.12), 
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2.1 TITLE ChROS 


This input permits labeling of runs with alphanumeric information. As 
many title cards as desired may be input in sequence. Card format is as follows: 

col 1-5 col 6-77 

TITLE any alphanumeric information 


Ee ee eR eR MTT NENT PLT TA SITE SINT; 2 TAM» MRRP ESR PRP ESET COT 
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2.2 THERMODYNAMIC DATA 


This data set is identical to the THERMO DATA described in Appendix 
D of NASA SP-272 (i.e. Reference 1), 


Using this data set, thermodynamic data curve fit coefficients may be 
read from cards. the curve fit coefficients are generated by the PAC computer 
program described in NASA TN D-4097 (i.e. Reference 7). 


The thermodynamic data (i.e. C. , etc.) are expressed as functions of 


Pr 
temperature using 5 least squares curve fit coefficients (a) 5) and two integration 


constants (a, _.) as follows: 


Co 
p 
—_ 2 a,+ aot + a,T®+a,T + aT* (2, 2-1) 
HS anT a,T a,T° a.T* a6 
WR go ata eee oe i 


~~ = 3, int +a,T+—2— + ad — + a, (2, 2-3) 


For each species, two sets of coefficients (a, _, and a} _,) are specified 
for two adjacent temperature Intervals, lower and upper respectively. For the 
data available in Reference 1 the lower temperature interval is 300° te 1000°K 
and the upper temperature interval is 1000°K to 5000°K. 


the input tormat required for this thermodynamic data {s defined in Table 
2-2, Data cards for the species H, HO2, H2, H2O0(L), H2O, H202, O, OH, O02. 
Al203(L), and ALZO3(s) are listed in Table 2-3 as exampics. Thermodynamic data 
coefficients for 421 species are available in Reference | and are supplied with 
the computer program, a list of these 42] species names is presented in Table 
2-4, 


Data Tape Generation and Usage: 


A computer run using thermodynamic data card input will generate e data 
tape on logical unit JANAF. This tape may then be saved and used at a later time, 


2-S 


The program writes the THERMO data card images on unit JANAF as read but with 
two minor exceptions. The THERMO code card and the card numbers in card 
column 80 are omitted. 


If thermodynamic data cards are not input, the program assumes the 
thermodynamic data is on logical wit JANAF. Logical unit JANAF is currently 
assigned a value of 25. 
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TABLE 2-2 FORMAT FOR THERMODYNAMIC DATA CARDS 


Contents 


3A4 
3F10.3 


1 to 6 
1 to 30 


THERMO 


Temperature ranges for 2 sets of co- 
efficients: lowest T, common T, and 
highest T 


Species name 
Date 
Atomic symbols and formula 


Phase of species (S,L, or G for solid, 
liquid, or vas, respectively 


Temperature range 
Integer 1 


Corfficients a! (i = 1 to 5) in equations 
(22-1) to (2.263) 


(for upper temperature interva)) 
Integer 2 


eee ee in equations (2.2-3)-(2,2-3 
f a5 for upper temperature intervai 
6 ay 85. and a, for lower) 


3A4 
2A3 
4(A2, F3.0) 
Al 


1 to 12 
19 to 24 
25 to 44 
45 


2F10.3 
115 
5(E15.8) 


15 
§(E15.8) 


15 
4(E15.8) 


Inteyer < 


Sea te-oit 2 equations (2.2-1)-(2.2- 
@,,%c,8-,4- for Jower temperature in- 
tear) © 7 


Integer 4 


Repeut cards numberec! 2 to 4 in ce 80 
for cach species 


®Gascous spccies and condensed specics with only one condensed phase can be 
in any order, Howover, the scts {cr two or more condensed phases of the same 
specics must be adjacent. If there urc more then two condonscd phases of a 
specics their ccts musi be either in increasing or decreasing order according 

to thelr temperuture intervuls. 


SABLE 2-3, SAMPLE THERMO DATA CARDS 


(Species AR, H, H,, H,0, N, 0, OH, O,) 


THERM 
300, 0006 1000,000 5§000,000 
AR L §/660AR 100 000 000 0G 300,000 §000,.000 
0,2500N0000E 01 0, Qe 0. 0. 
©0,74537507E 03 0, 63660006€ 01 9,25000000E 01 0, 
0. 0, ©0.7453749RE 03 0.63660006E 0} 
| J 9765H 100 000 000 0G 300.000 5000-000 
0,25000000E 01 0, Qo 0. 0. 
4: 2547 1627E 05=0,46011763E 00 9. 25000000E 01 0, Oe 
06 06 OcVS4TIGZTE 05~6,4601)742E 00 


H? J3/6iH 20 00 00 0G 300.000 $000.000 
O,DLOOLSOIE 03 0,51119664E°03 9.526442) NE=07=0.34909973E910 0, 36965345E-14 

00877380426 03-0,19629421E 01 He3057T4451E 0) 0,26765200E-02-0 , 580991 626005 
abated EM Sa la Ole cub ais oo ba 03-0, 22997056E =o 
H20 3036 100 000 06 300,000 
0,27167633€ 01 0, 4) et eet) et 4 0, aeenn2Ec0n8 40672145 Ea) 4 

, 20,29905R26E 05 0,66308871€ 01 A,40701275E 01=0,110806096-02 0.4615211806-05 
of 20637504508 O SOTO CIBC -3AuA BLMIVTIAGC OSC, 322700LSE OC 

N2 J 9765N_ 20 00 00 £06 *300, 060 5000,090 

0, 28963194E 01 0,15354RA6E-02-0,57235277E 206 0,99807393E~1090,652235SSE=14 

“NeSOSBHIBSE 03 OoHLEISIOVE 0) HeI6TAAZHIE 01°00) 20B81500E°02 6622240) 02E705 

=0563217559F.-0990022577:.53E 9120010611 SARE 04 0,235R04624E 0) 

0 J 6/62% «160 000 000 06 300-000 $000.000 
0,25420596E 0)-0,27850e 19E 0420.31 02R033F~08 0,45510674E=] 1=0,43680515E 915 
0029230R03E 05 Ocd49203LA0E 01 Ae294642R7E 01°0016381665E902 0024210316E°05 

$00 1602R432E708 Oo 3RGOEYASET]2 00291 47604E 05 0029639909E 0) 

OW J 34660 14 100 000 06 300.000 500020 
0,.29206627E 01 0,95931650F—030f,194461702E=06 0,13756666F =) 0 $° or e2238426e15 
0,39353R1SE 04 0,54423045E 01 §.30378963F 01-0, 107788565 =02 0.9683037AE=06 
0.38713972E-09=0, 2287109432 0236412H23E 04 0,493970009F 006 

02 J 97650 20 00 00 06 3002000 $000°000 
0.36Z19S3SE 01 0. 736182hSE9N3-H,1965222RE—06 0,36201550E 9} 0-0. 2H948627E 914 

@0,1201GR25E 04 0,36150960¢ 01 9,362959HSE 0120, 1B7H21A4F 02 0, 705565446005 

«0 ° 675351376008 0 «2155540361100, 106752268 04 0 o 43052770 0} 

END 


Die tO WMH SWS ODN OP WN FW FUN FW = 


oe ome 


iim aa Ti ae ai a in li a 


AL203(S) 
AL203(L) 
AR 

AR+ 

B(S) 

B{L) 

B 


B+ 


Table 2-4, 


BECLE 
BUCLz (S) 
BECL? (L) 
BECL2 
BLF 

BEF 2 (1:3) 
BEF 2(S) 
BEF 2 (L) 
BEF2 
BEH 
BEH+ 
BEN 


SPECIES WITH THERMODYNAMIC DATA PROVIDED 


BEO(S) 
BEO(L) 
BrO 
BLOH 
BEOH+ 


C2H 
C2HF 
C2H2 


CS+ 
CSCL(S) 
CSCL(S) 
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FECL2(S) 
rECL2(L) 
FECL2 
FECL3(S) 
FECL 3 (L) 
FECL3 


FEO2H2(S) 


PEO2H2 


FEO3 113 (S) 


FE203 (S) 
FE304 (S) 
H 

H+ 

H- 

HALO 
HBO 
HBO+ 
HRBO2 


H2 
H2Q(S) 
HZO(L) 
HzO 
HzO?2 
Hes 
HeB206 
HI 

Hi+ 
K(o) 
Kd.) 

K 
K+ 
KCL(S) 
KCL (L) 
KCL 
KF(S) 
KF(L) 
KF 
KF2~ 
K2r2 
KO 


KOH 


- KOH(S) 


KOH(L) 


LICL(S) 
LICL(L) 
LICL 
LIF(S) 
LIF (L) 
LIF 
LIF2- 
LIFO 
LIH(S) 
LIH(L) 
LIH 
LIN 
LIO 
LIO- 
LIOH(S) 
LIOH(L) 
LIOH 
LION 
LI2 
LI2CL2 
LI2F2 
LI20(S) 
LI20(L) 
L120 


MGF 
MGF2(S) 
MGF2(L) 
MGF2 


: 
| 
| 
| 


Table 2-4. (cont'd) 


MGH e) 
MGN O+ 
MGO(S) O- 
MGO(L) OH 
MGO OH+ 
MGOH OH- 


MGOH+ O2 | 
MGO2H2 O2- 


N P 
NF P(S) 
NF2 P+ 
NF3 PCL3 
NH PF3 
NH2 PFS 
NH3 PH 
NO PH3 
NO+ PN 
NOCL PO 
NOF PS 
NOF3 P2 
NO2 P4 
NO2- S(S) 
NO2CL S(L) 
NO?F S 
N2 S+ 
N2C SF4 
N2H4 SF6 
N20 SH 
N204 SN 
WA (S) SO 
NA(L) SOF2 
NA $02 


NACL S1(S) 
Nar (S) SI(L) 
NAF(L) SI 
NAF SI+ 
NAF2- SIC 
NAH SIC2 
NAO SICL 
NAO- SICL2 


NAOH SIF 
NA2 SIF2 
NA2CL2 SIT3 
NA2F2 SIF4 
NE SIH 
NE+ S1#H4 


2.2.1 THERMODYN/MIC DATA BELOW 300°K 

ror low temperature calculations it may be necessary to extend the 
curve fit data in the Therniodynamic Data file (see Section 2.2). The lower 
temperature limit, T,. in the Thermodynamic Data supplied with the program ts 


pes .° 
T= 300° K 


Thermodynamic Data bclow the temperature, Ty. may be input by data cards as 
described below. 


col ly 


L@W T CPHS 


Directive for start of 
low temperature CPHS 
tables (co! 1 through 10). 


; 

eas cers an bts tetea aerel ae eon 12 character species 
name, left justified 
followed by the integer, 
n, punched in column 2], 
The integer n must be such 
that 1 $n 53 and represent 
the number of Thermodynamic 
Data points to be input for 
this species, 

rk Cp Hi SP l First 'hermodynamic Data 

point for the above species, 
input 4F 10.0, I5. 
card n+2 P He Sr 2 nth (1sns3) nth Thermodynamic Dieta 
T point for the sbove species, 
Input 4F 10.0, IS. 


Repea’ cards 2 through n + 2 ahove for each species to be input. 
Temperature must be T ) <T,<T,<T Le 


END LYW T CPHS end directive (col 1 


through 14) 


fn example of this input is given in Table 2-5 which shows ao card sting extend- 
ing the Thermodynamic Lata for an 0,/H, propellant to 100°K. Data in Table 2-5 
is taken directly from tne JANAF tables (Reference 8), except for Argon which is 


taken from NASA SP-3001, | 
The quontity Hy is defined as 
r° : 
HR (HP -H} 9g) = Ale , cal/mole 
298 
Cp , cal/mole - deg K 
T 


Se , cal/mole - deg K 
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TABLE 2-5 . LOW TEMPERATURE Cp ; te Sr DATA FOR AN O,/ Hy PROPELLANT 
: a 


LOw YT CPRHS 

AR 2 

100,90 4.9681 ~984,.5 31,556 | 
200,0 4.9681 487.7 34,999 2 
sl 2 

100.0 4.968 511186 212965 l 
200,0 4,968 51614, 25.408 2 
H2 2 | 

100.0 5.393 =1265. 24.3R7 | 
200,0 6,518 =567,0 98.520 2 
H20 2 

100,0 7.96) ~59376,9 364396 1 
200.0 7,969 ~58581,9 41,916 2 
N2 2 

100.0 T2074 ~12387,0 98-113 | 
200, 6.989 684, 62.984 2 
n ? 

100.0 5.666 5R4796 32 0466 1 
200.0 5.434 59036, 362340 4 
OH 2 | 

100,0 7.567 7879. 95.852 1 
200,0 7,309 8623, 41.02} 2 
02 2 

1000 62958 ©13816 @12e395 1 
2000 6096) =6856¢ 460e2)8 2 


END Lor 1 cPHS 
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240 Geometry Directive and Input 


The GEQMETRY directive allows the user to input the nozzle geometry in 
the SGE@M namelist set and have that information transmitted to the appropriate 
calculation modules. Any or all of the values input via the GE@METRY directive 
may be overridden in the input to a specific calculation module, 


The format of the GE@METRY directive and description of the input is as 
follows: 


The word (SE@METRY must start in card column 1 (Only the first 4 charact- 
ers are checked). The card following the GE@METRY card should contain $GE@M 
in card columns 2 through 6. The cards following the $GE@M card may be in any 
order up to the SEND card which specifies the end of the geometry input data. 

All cards in the namelist input set must start in card column 2 or greater. The 
input items are 


Item Lescription Units Assumed 
Value (s) 

SGE@M namelist name 

AS UP(1) = supersonic area ratios at none 0.0 
which information will be 
printed out 

NAS UP = number of entries in the none 0 
ASUP array =< 40 

ASUB(1) = subsonic area ratios at none 0.0 


which Information will be 
printed out 


NASUB = number of entries in the none 0 
ASUB array <10 
NN@Z = number of nozzles none l 
RSI = initial throat radius, nt. in 0.0 
RCURV = normalized throat radius of none 0.0 
curvature, ¢ S/t*. 
*esp@T = nozzle erosion rate, r*. mils/sec 0.0 
*TBURN = motor burn time, used to sec 0.0 


calculate r* as a function 
of time, i.e., r*(t)ert+r*t. 


*mhese values will be overridden by the Grain Design & Ballistics Module cel- 
culated values if BAL=! or 2 in SPROB Namelist. 
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Item Description Units Assumed 


Value(s) 
IWALL = wall specification option none 0 
flag 
l cone option (input THETA, 
RCURV, THETAI, RSI, and 
EPS) 
2 parabolic nozzle contour 


option (input THETA, RCURV, 
THETAI, RSI, RMAX, ZMAX) 


3 circular arc nozzle contour 
option (input-same as in 
IWALL=2) 

4 nozzle contour (spline) option 


(input THETA, RCURV, THETAI, 
THE, RSI, RS, ZS, NWS) 


**Note** - Certain of the calculation modules (@DK and TD2P in particular) con- 
tain more extensive nozzle geometry specification options than are presented 
here. However, the options cited above are applicable to all of the calculation 
modules in the SPP code. 


alae 2 


-- oe 


The following items are illustrated in Figure 2.1 on the next page. 


THETA = nozzle attachment angle DEG. 0.0 

THETAI = nozzle inlet angle DEG. 0.0 

THE nozzle exit angle (input DEG. 0.0 
if IWALL=4 only) 

EPS nozzle expansion ratio (in- none 0.0 
put if IWALL=1 only) 

RMAX normalized radius at the none 0.0 
nozzle exit plane (input if 
IWALL=2 or 3) 

ZMAX normalized axial position none 0.0 


of the nozzle exit plane (in- 
put if IWALL=? or 3) 


RS (2) = table of norinalized wall radil none 0.0 
downstream of the nozzle 
tangency point. RS(1) is com- 
puted a8 rn. (Input if IWALL=4) 


ZS (2) = table of normalized axially none 0.0 
positions downstream of the 
nozzle tangency point. 2ZS(1) 
is computed at 2, (input if 
IWALL=4) 
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Units Assumed 


named —y 


carer 


NWS = number of entries in the RS, none 0 
ZS tables, Includes the first 
entry. NWS <20. (Input if 
IWALL=4) 


RZNGRM = Optional normalizing factor none 1.0 
for the RS, ZS table. For ex- 
ample, if RS & ZS were input 
as dimensional numbers, RZNORM 
would be the throat radius in those 
units. 


nee 


RI = normalized inlet wall radius none 0.0 
of curvature. Only used if 
GDK calculation selected. 
It is not necessary to input 
RI even if an @DK calcula- 
tion is to be made. 


ECRAT 2 Contraction ratio at which none 0 
the ballistics calculation 
will assume the nozzle in- 
let to start. 


FSUBM = le..gth of submergence/length none 0.0 
of internal motor. (used in 
caiculating the empirical sub- 
mergence loss for the nozzle). 


AEAS = nozzle entrance area/nozzle none 1.0 
throat area (used in calcula- 
ting the emplrical submergence 
loss for the nozzle). 


AVELS = motor average L* in. 0.0 


KN@Z e nozzle type flag none 0 
=l1 stee] nozzle 
=0 all other types nozzle 


SEND 
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2.4 Tr ry Directive and In 


The TRAJECTQRY directive allows the user to input either the altitude or 
exhaust ambient pressure as a function of time during a motor firing in the $TRAJ 
namelist set. This data allows the SPP computer program to calculate the de- 
livered performance at the actual motor conditions. If this input is not speci- 
fied only vacuum performance can be calculated. 


The format of the TRAJECTORY directive and description of the input is 
as follows: 


The word TRAJECTORY must start in card column 1. (Only the first 4 char- 
acters are checked). The card following the TRAJECTORY directive should contain 
STRAJ in card columns 2 through 6. The cards following the STRAJ card may be 
in any order upto the SEND card which specifies the end of the trajectory input 
data, All cards in the namelist input set must start in card column 2 or greater. 
The input items are 


Item Description Units Assumed 
Nalue(s) 

STRAJ namelist name 

PAMBI(1) = amblent pressure table PSIA 0.0 

H(}) m altitude table FT, 

TTRJ(1) = times corresponding to the sec, 0.0 


ambient pressure or altitu- 
de tables. Must be mono- 
tonically increasing in tlme. 


NTJPT z number of entries to be used none 0 
in the Input tables, NTJPT= }. 
If NTJPT#1, then the following 
is done: H(2)=H(1); PAMBI(2)= 
PAMBI(1); TTRJ(2)=ix1l0, 


IPH = 0 - pressure table is to be none 0 
used 
1 = altitude table is to be 
used 

SEND end of S$TRAJ namelist set 


Note: The PAMBI and H arrays have been equivalenced to each other in order to 
save core storage. Hence, if both array names are input, the last input 
values will be used in accordance with the IPH flag. 


2S Problem Lirective and Input 


The PR@BI EM directive and $PR@B namelist input allows the user to 
specify which calculation modules are to be executed and which losses are to be 
considered. The occurence of the PROBLEM card and associated data in the input 
stream causes the selected loss modules to be executed in the sequence Shown 
in Table 2.1. 


The format of the PR@BLEM directive and description of the input is as 
follows: 

The word FROBLEM must start in card column 1. (All charecters are 
checked). The carc following the PROBLEM card should contain SPR@B in card 
columns 2 through 6. The cards following the $PR@B card may be in any order 
up to the SEND card which specifies the end of the problei: input data, Ali cards 
in the namelist input set must start in card columns 2 or greater. The input items 
are 


Iter: Description Units Assumed 
Value(s) 

SPROB namelist name 

QDE = these input items call out the cal-  — none 0, 


culation nodules which are to ex- 
ecutead or specify how the dita wick 
would have genersted by these ray On 


CAL * modules is tc be obtained. ‘the 3 
| allowable values for these vartables 
< are: : 
BOK a= ; yone e 
Jik 0.0 elther not concivered ur ; 
} not exec itad. 
TDP ‘ | 1.0 module to t2 be executcd, none 0, 
| 2.0 data sreviuusly generated by 
wOtule ss to be read in or 
es iogtcal aniz INF, one 0 
baads = os e 
3.0 valucs which would pe gen- 
crated Lv toduies are to be 
read in, 
sREQ = 9 kinetic sti.ct conditions for none 0 


JDK are to be read in via 
the SPECIES gairective. 


{ kinetic start conditions for 
QDK are to be computed 
irom a'restricted equilibrium 
QDE calculation’. (see 
ECRAT in SODE NAMELIS7) 


Item Description Units 


Assumed 
Value(s) 
ETAI (1) rs input values of the loss mechanism none 1.0 
efficiencies fractions. (Used only 
if module option =3, is selected) 
I Mechanism 
1 finite rate kinetics 
2 2 dimensional-2 phase 
flow 
3 combustion efficiency 
submergence 
5 2 dimensional divergence 
loss 
INP = Unit number on which previously none 3 
generated data is to be read in on, 
(See Section 2.11) 
SEND = End of namelist set SPR@B 


2.6 @DE INPUT DATA (ALL PROBLEMS SPECIFYING @DE=1,0 OR IREQ=1) 


The @DE input data described here its as defined in NASA SP-273, Refer- 
ence 1 , except namelists input SINPT2 and SRKINP have been combined into a 
single list named $@DE, Only the rocket option (RKT in the namelist) is described 
here since all the other options have been removed from the ODE calculation moduie. 


The @DE input data consists of the following input groups: 


1, REACTANTS directive card, followed ny up to 15 data cards, 
followed by a blank card, specifying reactants. 

2, @MIT and INSENT directives to omit or insert species for equilibrium/ 
frozen calculations. 

3. NAMELISTS directive card followed by input namelist S@DE 


specifving input case data. 


2.6.1 REACTANTS CARDS 


This set of cards is required for all @DE problems. The first card in the 
set contains the word REACTANTS punched in card columns 1 to 9. The last card 
fu the set {s blank. ‘n between the first and last cards may be any number of cards 
up to a maximum of 15, one for each reactant species being considered. The cards 
for each reactant must give the chemica! formu!a and the relative amount of the 
reectant. For some problems, enthalpy values are recutred, The format and con- 
tents of the cards are summarized ir: Table 2-6, A list of some REACTANTS cards 
is given in Table 2-7, 


Relative amount. of reactents. = The relativ< amounts of reactants may he 
specified In severa] ways. They may be specified in tenis of moles, mole fraction, 
or mole percent (by «eypunchirg M in cer? column 93) or In terms of weisht, wetaht 
fraction, or weight percent (blank in co: cna $4), 


Reojative amourts of total Drei to totul o¢iccints can also ve input. For this 
Situation, oach reactaat must Le specifies as 4 cuel of cn oxicizer by xeypunch- 
ing on For O, respectively, in columi °2 of the REACTANTS card. The amounts 
Given on the REACTANTS carcis are relauve -o tutal fuel or total oxidant rather 


than tatal reactant, 


TABLE 2-6, REACTANTS CARDS 


| Order | Contents 


First | REACCANTS 


Any | One car: for each reactant specius 
(maximum 15). Each card contains: 


(1) Atomic s«mbols and formula num- 
bers (maximum 5 seis)* 


(2) Relative weight? or number of 
moles 


(3) Blank if (2) is relative weight or 
W if (2) is number of moles 


(4) Enthalpy or imernal energy®. 
cal mole 


 {8) tule. S LL. or G for solid, 
hieuid or cas. respectively 


(6) Tempera'ure associated with 
enthalpy in (4) 


(7) F iffuelor O if -midant 
(8) Density ing em® (optional) 
Last | Blank 


Sa ne eS EE eee 


[Format _| 
3A4 


84 to 71 


72 
73 10 8) 


*Praerom with caleulate the enthalpy or imterart energ: (4 for sorcies ta 
the THERAR) data at the temperature (6) of zeras ire punched tn card 
colunins 37 and 38. (See section Reactani enthalpy fur additions! in- 


lurwition.} 


bRetative weight of {vel in total fuels ur cxidunt i total axidams. All 
reacties must be given either all in relative wer-hts ur all in rember 


wf moles. 
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TABLE 2-7 LISTING CF SAMPLE REACTANTS CARDS 


REACTANTS 

Ww 2, 100, 0. 6296.15 F 

N .7808810 .209795A4R, 004662 100, ~7,2021646298,15 0 
REACTANTS 

N ae H 4, Chie 0 4. 72,06 70730, $298,185 F 

C p a a] 1.869550 0031256S 2008415 18,56 2999 ,082L298,15 F 

ALS, 9,00 +C,0 $298,415 F 

MGi, 0 i. »20 143700, §296,15 F 

Wh 2s O ie 016 6831744 1296.15 F 
REACTANTS 

H 2, 00 100, 0, 6298,15 F 

0 2, 00 100.0 6.0 6298,15 9 
REACTANTS 

N 2, nw 6, C 2; 50,0 42734,8 (298,15 F ,78&6 
N 2, M4, 50,9 12050, .296,15 © 1.003 
F ?, 100, @3630,892L 85,24 Q 1,54 
REACTANTS 

Lad, 100, Ce $298,15 * 

F 2. 100, ©3930,892L aS ,24 0 1,54 
REACTANTS 

N 2, H 4, 60, 42100, L.296,15 F 1,003 
BES, 20. 0.0 $298.15 F 1.85 
4 2, 0 2, £00, ©44860, 1296.35 0 1,497 


*Listed above are six examples. Each example must end with a blank card. 
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There are four options in the $@DE namelist for indicating relative 
amounts of total fue’ to total oxidant as follows: 
1. Oxidant to fuel weight ratio (@F is true). 
2. Equivalence rutio (ERATI@ is true). 
3. Fuel percent by weight (FPCT is true). 
4, Fuel to air or fuel to oxidant weight ratio (FA is true). 
For each option, except @DE with @DK#1, and TD2P41, the values are given in 
the @FSKED array of S@DE (described in Section 2.6.4). For @DE=1, @DK¥1, and 
TD2P41, the MIX array is used as described in Reference 1. 


Reactant enthalpy. Assigned values for the total reactant are calculated 
automatically by the program from the enthalpies of the individual reactants. 
Values for the individual reactants are either keypunched on the REACTANTS cards 
or calculated from the THERM@ data as follows: 


Enthalpies are taken from the REACTANTS cards unless zeros are punched 
in card columns 37 and 38. For each REACTANTS card with the "00" code, an en- 
thalpy will be calcuiated for the species from the THERM@ data for the tempera- 
ture given in card columns 64 to 71. 


When the program is calculating the individual reactant enthalpy for values 
from the THERM@ data, the following two conditions are required: 
1, The reactant must also be one of the species in the set of THERM@ data. For 
example, NH, (9) is in the set of THERM@ data but NH,() is not. Therefore, 
if NH, (g) is used as a reactant its enthalpy could be calculated automatically, 
but that of NE, (°) could not be. 


2. The temperature T must be in the range 7, /1,.2s Ts tare 1.2 where 


“low Tow 


to Ty, igh Is the cvemperatare range of the THERMG@ data. 

For many cises it may be desirable to modify the enthalpy entered on the 
REACTANT cards. ‘his cax be done by using the DELH input array. The DELH entry 
will be added to the system enthalpy as co.iputed by @DE from the reactants cards 
(see above). For example, overall system enthalpy of the propellant can be input 
through the reactants cards and the energy aided or extracted due to conditioning 
the propellant can be input by the DELH entry. An alternate method would be to 
input zero enthalpy on the Reactants cards and input enthalpy by the DELH entry. 
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The second and less important reason is that if one knows that one or 
several particular condensed species will be present among the final equilibrium 
compositions for the first assigned point, then a small amount of computer 
time can be saved by using an INSERT card. Those condensed species whose 
chemical formulas are included on an INSERT card will be considered by the 
program during the initial iterations for the first assigned point. If the INSERT 
card were not,used, only gaseous species would be considered during the initial 
iterations. However, after convergence, the program would automatically insert 
the appropriate concensed species and reconverge. Therefore, it usually is 
immaterial whether or not INSERT cards are used. For all other assigned points 
the inclusion of condensed species is handied automatically by the program. 


NOTE: For all metalized propellants it is strongly suggested that the appropriate 
condensed phases be included in the insert cards to insure convergence 


in the equilibrium calculation. 


2.6.2 OMIT and INSERT Cards 


@MIT ard INSERT cards are optional. They contain the names of particular 
species in the library of Thermodynamic Data for the specific purposes discussed 
below. Each card contains the word @MIT (in card columns 1-4) or INSERT (in 
card columns 1-6) and the names of from 1 to 4 species starting in columns 16, 
31, 46, and 61. The names must be exactly the same as they appear in the 
THERM@ data. 


2.6.2.1 OMIT Cards 


These cards Jist species to be omitted from the THIRM® data. If OMI 
cards are not used, the program will consider as possible species all those 
species in the THERM®@ data which are consistent with the chemical system 
being considered. Occasionally it may be desired to specifically omit one or 
more species from considerations as possible species. This may be accomplished 
by means of OMIT cards. . 


2.6.2.2 | INSERT Cards 


These cards contain the names of condensed species only. They have 
been included as options for two reasons. 


The first and more important reason for including the INSERT card option {s 
that, in rare instances, it is impossible to obtain conver jence for assigned enthalpy 
problems (HP or RKT) without the use of an JNSERT card. This occurs when, by 
considering gases orly, the temperature becomes extremely low. In these cases, 
the use of an INSERT card containing the name of the required condensed species 
can eliminate this kind of convergence difficulty. When this difficulty occurs, 
the following messaje is printed by the program: "LOW TEMPERATURE IMPLIES 
CONDENSED SPECIES SHOULD HAVE BEEN INCLUDED ON AN INSERT CARD”. 
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21060 S@DE NAMELIST INPUT 
The @DE subprogram contains namelist input sections $@DE. The Name- 


list S@DE must be preceded by a card with NAMELISTS punched in card columns 
1-9, 


The S$@DE Namelist is required for inputs to the $PRQB namelist set if: 


@DE = 1.0, 

@DK = 1.0, 
or 

IREQ = 1. 


If both @DE=1, and IREQ=1, are specified in the SPR@B namelist, two sets of ODE 
data will be required. 


The variaijes input by the S@DE namelist are listed in Table 2-8. Addi- 
tional information apout some of these variables follows: 
Pressure units. - The program assumes the pressure in the P schedule to be in 
units of atmospheres unless PSIA =.TRUE. 
Relative amounts of fuel(s) and oxidizer(s). - These quantities may be specified 
by assigning | to 15 vlaues for elther o/f, %F, f/a, or r. If no value is assigned 
for any of these, the program assumes the relative amounts of fuel(s) and oxidizer (s) 
to be those specified on the REACTANTS cards. (See discussion in REACTANTS 
Cards, Section 2.6.1). 


RKT problem. - Orly one value for chamber pressure, P, is to be input. 


Print out wili be given for the chamber pressure condition (i.e. stagna- 
tion) and the throat condition. Print out may be renuested at other conditions by 
use of the PCP schedule and the SUBAR and SUPAR schedules. The values of ASUP 
ana ASUB from the $SGEOM namelist set will automatically be inserted into the SUPAR 
and SUBAR print arrays. Thus any input into these arrays will override the values 
input ir the SGE@M namelist set for the ODE calculation only. 


The progran wil] calculate both equilibrium and frozen performance un- 
less FR@Z=F or EQL=F are input. If FR@Z=F, only equilibrium performance will 
be calculated. If EQL=F, only frozen performance wil! be calculated. 


TABLE 2-8, VARIABLES IN $@DE NAMELIST 


RA yr TH 


Variable No, of Type Definition and comments 

entries 

! 

| 

RKT l L True Rocket problem 

P(I) 25 R 0 Assigned pressures: stagnation 
pressures for rocket problems: 

[ values in atm unless PSIA, or SI= 
.T., (see below) 

SI i ib False “Values in P array are in N/m? 

PSIA l i False “Values in P array are in psia units 

XP(I) 50 R Ll. Not currently used 

QF | L False Oxidant to fuel weight ratios are 
to be input9 

ERATIZ l L False Equivalence ratios are to be input® 

FPCT j i False Percent fuel by weight are to be 
input? 

FA ] Ne False Fuel to air weight ratios are to be 
input@ 

@FSKED() 50 R 0 For a Roc! et problem, and ODK= 
l or TDZ2P=1, @FSKED will be used 
rather than MIX (see Reference 1). 
Relative amounts of total oxidant 
to total fuel are input as defined 
by GF, ERATI@, FPCT, or FA. 

DELH(I!) 5¢ R 0 This value .vill be added to the 
system enthalpy input thru the re- 
actants cards. Units are BTU/ 

# if PSIA=.T., otherwise cal/gram. 

IQNS l L False | Consider ionic species” 

PCP(!) 50 R 0 Compute end print solutions at these 
ratios of chamber pressure to pres- 
sure (entries must be <1,) 

SUBAR(I) 50 R 0 Compute and print solutions at these 
values of subsonic area ratios (en- 
tries must # 1.) 

SUPAR(I) 50 R 0 Compute and print solutions at these 


values of supersonic area ratio (en- 
treis must # 1.) 


orf variable is set to be true. 


DSet variable false ii these calculotions are not desired. 
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Table 2-8 (cont'd) 


Variable No, of Type 


entries 


ECRAT ] 
EQL ] 
FROZ l 
LIS TSP } 
KASE } 
AQLPT ] 
NLPRNT ] 


a ; : 
“If variable is set to ve truc. 


Definition and comments 


Subsonic area ratio to start @DK 
calculations with computed equil- 
ibrium conditions. The SUBAR in- 
put table must include an entry 
equal to ECRAT, 


Calculate rocket performance as- 
suming equilibrium composition 
during expansion”, 


Calculate rocket performance as- 
suming frozen composition during 
expansion”, 


List names and dates of all species 
residing on thermodynamic data 
used, 


Optional assigned number associated 
with case. 


For cases with condensed phases, 
the equilibrium conditions at the 
onset of solidification will be cal- 
culated. 


Write out the S@DE namelist data 
set using a namelist write.? 


“Set variable false if these calculations are not desired. 


“Required only for the restricted equilibrium cption, i.e., IREQ=1 in SPR@B namelist. 


“Assumed falce for a restricted equilibrium calculation, 
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Dal Grain Design and Ballistics Input 

The input for this module is described in three parts. Section 2.7.1 de- 
scribes the $BAL NAMELIST Input, while Sections 2.7.2 and 2.7.3 describe the 
formated grain design and formated ballistics input, respectively. 


2.7.1 $BAL NAMELIST 


In the following list of input quantities an * preceding the variable name 
means that it is internally transmitted from the ODE module if ODE = 1 or 2 in the 
$PR@B NAMELIST ‘at present only values at one chamber pressure are transmitted 
to the BAL module from the @DE module). Two ** are used to denote variables that 
are transmitted from the $GE@M NAMELIST, if used. Values input directly in SBAL 


will over-ride internally transmitted values. 


Variable Description Units 
*PC Estimate of initial chamber pressure Ib/in* | 
*CSTART Theoretical C* ft/sec 
*TFT T.enretical fleme temperature °K 
*CPG Soecific heat of the gas only Btu/mol/”R 
*CPM Specific heat of combustion products Btu/mol/R 
*GAM Ratio of specific heats - 
*XMG Moles of gas per 100gm of products mol/100gm 
*XMM Moles of products per 100gm of products mol/100gm 
*FCP Moles of condensed phase per 100 gm of mol/100gm 
products 
*XMV Gas viscosity lb/in/sec 
*XKG Gas thermal conductivity Btu/in/sec/°R 
*XB Nozzle erosion narameter, R - 
**RN(1) Nozzle entranc: radius in 
* *RN (2) Initiai nozzle throat radius in 
#* RN (3) Nozzle exit radius in 
DTIME Not used 
KBAL Logical unit number for internal informa- ~ 


tion transfer (set equal to 8 in MC module) 


MODE Integer used for control of internal informa- - 
tice transfer. (Set equal to 0 in MC module) 


**NNOZ Number of nozzles - 
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2.7.2 Grain Geonnetry Inputs 


The data cards required to define the initial grain geometry are divided in- 
to two categories: (1) Initial data cards and (2) Figure cards, This input must be 
preceded by a title card, This title card is not used, however. The title card de- 
scribed in Section 2.1 is used instead, 
2.7.2.1 Initial Data Cards 


INITIAL is »unched on the cards starting in column 1, These cards con- 
tain data which define the region tn a coordinate system where the surfaces wil] 
be located plus data which control the limitations and X, Y and bum increments. 
The burn interval may be divided into a maximum of four sub-intervals; e.g., BO, 
Bl, B2, B3, B4, where the Kth sub-interval (the interval defined by By}. B,) 
would be divided into NB, parts. This may be desired where, for example, it is 
obvious that a complicated shape will bum into a simple shape, warranting a 
change iu mesh size to save computer time. 


The subscripting for X, NX and Y, NY is also identical*® ana is related to 
the burn intervals cs follows: for each burn sub-interval, a complete set of X and 
Y values must be prepared, thus the reason for the double subscript on these 
quantities, Tre procedure for the first subscript is the same as for the bum and 
burn increments and the second subscript corresponds to the burn interval, To 
illustrate, Ssupmnose the burn interval is divided into two sub-intervals BO to Bl to 
BY with NBl ard NB parts; thus, two sets each of % and Y data are required, as 
follows: 

Interval BO tc Bi: 

XO1, X11, NAIL, X2l, NX21, X31, NX31 up to X41, NX41* 

NOW leven Vl ae al ee ee ee Le 
Interval Bl to B2: 

ROD GA eer Ie pare) NMA Y ee ae, Nc IN 

NO Zi Nae ULe sate g tie pen re 
*EXCEPTION: When the only burn interval is from BO to B! it is possible to have ip 


to sixteen changes in the X Interval, ‘".e subscripting on the X input would be as 
follows: 


From one to nine in.erval changes: rrom ten to sixteen changes (continuinu) 
X01 X21] NX11  # X6] NX61 X03 NX23 X43 NX63 

X21 NX21 X71 NX71 X13 NX33 X53 NX73 

X31 NX31 X81 NX81 X23 NX43 X63 NX83 

X41 NX4i X91 NX91 X33 NX53 


XSL NX51 
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where X01 = beginning X 
X11 = ending X 
NX1]1 = number of Intervals in between X01 and X11 


X21 = next ending X 
NX21 = number of intervals in between X11 and X21 


etc. 
If too large a mesh size is specified a plot of surface area, or pressure, 


versus burn distance, or time, will reveal a saw tooth effect. 


P effect of mesh being 
too coarse 


actual motor 


t 


To avoid this situation the maximum X and Y increments should be based on the 
following criteria: the maximum X ana Y increments, in ccder that a step or saw 
tooth effect in the plot of surface area or pressure versus burn distance or time 
is not obtained, are: 

max AX = AB csc A 

max AY = AB sec A 


where @is the angle between the X axis ar.d the line tangent to the point of in- 
tersection of the actual void and the representation figure that constitutes most 


to the surface area in the region. 


actual surface 


fitting surface 


TR EERE = 
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It is not necessary for the number of increments of X and Y corresponding 
to the first burn inte:val to equal the number of increments of X and Y correspond- 
ing to the second burn interval and also for any more burn intervals. The restric- 
tions on X and Y are that the initial and final values for X are equal and the initial 
and final values for Y are equal. To illustrate, suppose the following case: 


For the interva! BO to Bl: 


X01 = 0.0 
X11 = 5.0 
NX11 = 10.0 
X21 = 10.0 
NX21 = 6.0 
For the interval Bl to B2: 
X02 = 0.0 
X12 = 10.0 
NX12 = 10.0 


In this case X91 = X02 and X21 = X12. The same idea fo:lows with Y. 


The number of bum interva!s must be a whole number, and the number of X and 
Y increments must be whole numbers, and the X and Y increments must be even whole 


numbers, 


a 7 ere we 


The burn interval, or sum of burn sub-intervals, constitutes the total pro- 
pellant web, The aumber of burns, N8, determines the frequency of geometry out- 
put; thus, NB 20 will yield answers ror each 5% of web burned. This {s not to he 
confused with NX or NY, wkich define the calculational resh at a given time or 
wep burned, of with separate conirol of ballistics output. 

The computational time will be p-oportional to the product of the number 1f 
X anid Younerements, In other words, if both are doubled, the time will incresse by 


roughly @ factor of 4. 


2.7.2.2 Figure Cards 

The first card for cach surface contains the names of the figure starting in 
colunn 1, The cards immediately following contain the data defining the figure 
and how it will burn, The first card contains Information 2f whether the figure is 
to consist of propellant or void. If the figure is to consist of propellant, a CR is 
punched in columns ‘9-20 and the figure is considered to burn inward; otherwise 
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the figure {s considered void and burns outward. 


The burning of the figure is controlled on a BURN card, A figure may burn 


in any combination of the 


following types of buming. (Tnis card ts used only if 


a figure {s to be burned in a special manner.) 


Vartables 
Type Name Value Burn 
Burn Columns Columns Considered Action 
Nonburning NB Burn = 0.0 No burning occurs in 
the figure. 
Delayed burn DB a Burn = B-a The figure is delayed 
in burning the speci- 
fied amount (a). 
Rated burn RB b Burn - (B)(b) The figure is burned 
b times as fast (or 
slow) as the normal 
burn. 
Corner rounding CR Cee nnne The figure (prism, 


The data defining 


cone or cylinder) is 
backed-up the distance 
c and then bumed out 
again the distance. 


the figure is contained on the DATA cards. Punched in 


columns 1-4 of these cards is the word DATA. These cards will contain the control 
characters Xl, X2, X3, Yl, YZ, Y¥3, Z1, 22, 23, R, Rl, R2, H. The X, Y, and Z 
characters define points in an X, Y and Z coordinate system. R is the radius and H 
is for prism height. In the case of a cone, where two radii are required, Rl and 

R2 are used. When R1#R2; e.g., In a cylinder, the Ris used, which will cause 
the same value to be stored in R} and R2, 


The figures the pr 
Gone (right) 


Cylinder 
(ight) 


ogran will handle and their inputs are: 


Two points (the center of two bases) X1, Yl, Zl 
X2, Y2, 22 


Two radii (the radius of the two bases) Rl, R2 
(Rl is radius at point 1 and R2 Is radius 
at point 2) 


Two points (the center of the two bases) Xl, YI), Zl 
X2, Y2, Z2 


One radius (the radius of the cylinder) R 
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Prism (right Three points (corners of bisecting X1, Yl, Zl 
triangular) plane) X2, Y2, 22 
X3, Y3, 23 
Height of prism H 
Sphere One point (center of sphere) X1, Yl, Zl 
Radius of sphere R 


2.7.2.3 Propellant Periphery and Exposed Insulation 


The complete description of the motor geometry requires a series of cards 
in the form of a table of outer propellant radius versus axial position. The X con- 
vention should be consistent with that describing the configurations above, but 
the number of entries need not be consistent with NX. Only the extremities in X 
need be consistent with those on the initial *INIT) cards. Up to 100 entries are 
allowed: however, far fewer will generally suffice. 


Insulation area exposed as the web burns {is not computed, but may be in- 
put in the form of a table. (See Section 2.7.3, Treble 2-10) 


2.7.3 Ballistics Inputs 


2.7.3.1 Input Scheme 


Following tt:e last card containing grain desian information, but preceding 
the ballistics input, 4 card containing the word COMPUTE, starting {n card column 
1, must appeer, 


The ballistics input is controlled from subroutine TABIN, but is actually 
read in subroutine DATAIN, in the following sequence: 


TABLE 
COMPLETION CODE 


DATA INPUT WORDS 
SET COMPLETION CODE 
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TRANSFER N-Z WORDS 
TO WORKING REGION 


DATA INPUT WORDS 
SET COMPLETION CODE 


[RETURN TO MAIN3 
All of the data tables (1) are read first, with a limit of 50 tables and 3000), words 
of storage, both about three times more than normal input for one case, The tab- 
.'¥r input ts assumed to continue until a READ IN completion code (loc. 4=-1) is 
given. The input in the remaining two sets (2) and (3) is then interpreted tn the 
program. Integer words (2) are logic controls for the calculation, and data words 
(3) include coefficients, constants, and some problem logic control words, to be 
described in detail] in subsection 2.7.3.2. 


The input is transferred te working storage in the main program (4) if non- 
zero (nz) values are given. As a result, only changes in the problem are needed 
in successive cases. Also, only input pertinent to the case and options to be 
used is required, s> unneeded input can be completely ignored. Since zero input 
is ignored, no “filling” is required, but values in the cards do no harm If desired 
for column alignmert, etc. One limitation in the method is that options can not 
be “tumed off" by inputting a zero in a later card. The node points are derived from 
the NX values in the grain design input. Points defining radial slots include one 
point at the end of each grain segment and one (only one) in the slot. The number 
of NX increments is given as -2.0 for the special case of a radial slot to key the 
logic. Once accomplished, the number is made positive and the grain design 
program operates as u3ual, 


2.7.3.2 Detaiis of Input 


In this subsection, the mechanics o! card format and table input are con- 
sidered; then the parameters and requirements are specified. 


Format 
Three card input formats are used by the program, 
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4 (12, 14, 64), 12, 28A Logic, first of set 
6 (12, 14, &X), 12 Logic, remaining cards In set 
6 (i2, F10), 12 Data 


For either logic (integer) words or data cards, the input is interpreted as (loc for 
location). 


1-2 3-12 13-14 aes 

loc, input loc wom 
value 

[2, 14, 6X or F10 12 =e 


The input cards are read In seis, and a set is assumed to continue until 
a completion sign is encountered. Two completion signs are recognized: 
(i) A blank card 


(ii) A negative location, i{.e., location = -1 in any of the [2 fields 
on the card 


Note that the completion sign for the group of tabular inputs is different in mean- 
ing and form ( 4b -, in l2, 14). 


Rules governing the interpretation of a card set include: 
Any number of cards may form a set 
Any order of locations is acceptable 
Zero values of input are ignored 
If locations are repeated, the last value of input is used (overwrites) 


If a location is blank but data is given, a sequenced list {is assumed, 
so the following two cards are cguivaient: 


] 10. 2 ll, 3 12, -1 
l 10, iif 12. -1 


The sequencing doer not follow from card to card, so the first location on the 
next card must be given, 


Tabular Data 


Deta tables are assumed of the form z=f (x, y). Tables are termed rec- 
tangular if the number of x arguments in the rows is the same for each of the col- 
umns of y arguments. A non-rectangular table, then, has a varying number of x 
arguments in the colurans. 


The card sets to irput a table are given in Table 2-9 with the location 


So Sn NIT ARRAN. sa MR IONND DCR NNCRS LAME ce MAI dnee <n © oe 
a2 ea whee 


Table 2-9 input of Tabular Data 


Input for Form Z = f (X, Y) 


Input Type Description 


Integer card set Required Control constants 
Data card set Required if NY >1 Omit if NY=1 Y values 


Data card sets 1 set if rectangular NY sets if X values 
nonrectangular 


Data card sets NY sets required Z values 


Integer Cards 4 (I2, 14, 6X), 12, 28A; 6(12, 14, 6X) 12 


Description 


Required Table number (see Table 2-10) 
Optional Key to LOG-LOG interpolation (+) 


Optional Table type non-rectangular = (+) (Note that (-) 
is key to end table inputs) 


Optional Number of Y (columns) values, required if NY >1 


Optional Number of X (rows) values, required if NX>1 


Optional Number of X values in 2nd ---- 19th Ycolumns, 
non-rectangnilar tables 


Data Card Sets 6 (12, Fl0), 12 


1 


Optional List of X, Y, or Z values se yuenced from Location 1 
Required Completion sign 


NOTES: 1, Atable can degenerate to one point in X, Y, or beth. 


2. Interpolation is linear with no extrapolation: values at boundary 
of table are returnec, 


3, A missing table causes Z=0 and one (only one) error message. 


of the words in the sets. The logic (integer) card set is always required and must 
include a table number and a completion sign. The number of y columns (NY) in 
location 5 and x rows (NX) in location 6 must be given if more than one. Location 
3 is used to key log-log interpolation if desired, and a (+) value in location 4 keys 
a nonrectangular table. 


Usually, one set of y values, one set of x values, and NY sets of z values 
compose a table. If only one column is specified, the set of y values is omitted. 
It a non-rectangular table is specified, the number of x values in each column 
must be given, starting in location 7, and an x argument set is required for each 


y argument, 


Assignment of Input Tables 


The numbers of input data tables are assigned in the subroutine MAIN3 due 
to compiled-in constants in the call statement: 

CALL TAB (X, Y, Z, NO) 
The correspondence of input data and the table numbers, NO is given in Table 
2-10 with the FORTRAN symbo)] and assumed arguments of the tables. In Table 
2-10, data that is designated as “optional” is needed only if the related option 
is specified. Data that is designated as requirec is almost always needed since 
zero values of the parameter (Z) are rarely meaningtul; the possible exception is 
the list of times of datailed output where no input results in printout at all time 
steps, 


Integer Inputs 


A series of integer words may be input to control program options, as shown 
in Table 2-11, Loc (4)-NUM(4)=-1 is the required input to indicate the end of data 
tables, ‘she none number of each radial slot is required input. Other input in Table 
2-11 is used only if the particular option is desired, 


Data Input Words 


Up to 99 data (real) word locations are provided to complete the input 
constants and contro: options in the program, as iisted in Table 2-12. Note 
that one data set is used, regardless of the number of cards, as discussed earlier. 
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Table 2-10 Input Table Assignments 


Propellant strand burn- Required 
ing rate at conditioning 
temp. 


Pressure exponent Required 


Motor geometry (grain Required X123 
outer radius) 


Insulation area Optional TIME 
Time increments Required TIME 
Time of detailed output Required TLIST TIME 


Throat area table Optional AP (NN); TIME 
KODE=5 


Slot burn rate correction Optional DR WEB (I) X(I) 
Slot diameter Optional DSLOT TIME AFI(L) 
Slot burn-pack correction | Optional DX WEB(I) | AFI(L+10) 


Table 2-11 Integer Input Parameters 


Description 


(=-1) indicator of end of tables — 
Key (+) for ignition dt option 
Node number of radial slots 


Propellant number for following 
segment 


Num (4) 
Num (5) 

Num (11-20) 
Num (21-30) 


Required 


Optional 
Optional 
Optional 


eee ee ee a 


SE ee eee RN ae ee ae 


Z(1) 


Z(2) 


Z(3) 
7K) 


217) 


Z(8) 


aa tht tin a in il i ak ale EG aa ll alata lala, lik atari aatar lalallala siete aatla nal eaiaaaiae 
~~ 


Z(3) 


(10) 
Z(11) 
Z(12) 
Z(13) 
Z(14) 
Z(15) 


Z(16) 
Z(17) 
Z(18) 
Z(19) 


Z(20) 


TSTOP 


TRISE 


CONV 


PAM 3 
KODE 


KODE 1 


KODE 2 


per cent 
gm/mole 


lb/in® 
fractional 


lb/in® 
Ib/in® 


tractional 


Required 


Optional 


Required 
Required 
Required 


Required 


Required 


Requires, 
Required 
Required 
Required 
Required 
Required 


Required 
Optional 
Optional 
Optional 


Required 


2-40 


Table 2-12 Data Input Parameters 


Calculation run (expected burn) 
time 


Time of rapid pressure rise 
(ignition rise), an ald to con- 
vergence 


Convergence ratio (.00005-.005) 
Ambient pressure 


Nozzle material 

» carbon phenolic 

» ATJ Graphite 

- pyrolytic graphite 

e non-eroding 

. @ table of erosion rate is to 
be input 


1. if fully insulated motor 
2. if uninsulated 


Binder type 

1. PBA Binder 

2. HTPB Binder 
3. CTPB Binder 
4. NC Binder 
5. FU Binder 
6. PGA/NC Binder 
7. NF Binder 


On kd GS DO im 


Aluminum concentration 
Molecular wt. metal oxide 
Ambient temperature 

Loop gain ratio 

Trial value of motor pressure 


Nozzle efficiency; 1.000 in 
the integrated mode 


Propellant density 
Insulation density 
Insulation ablation rate 


Projected nozzle surface area 
for reradiation to propellant 


Recovery efficiency for radial 
slots 


Z(21) 
Z(22) 
Z(23) 
Z(24) 
Z(25) 


Z(26) 
Z(41)-Z(50) 


2(71)-z(80) | 2 


Z(91)-Z(100) 


in/sec Required 
== Required 
lb/in® Required 
BTU/1b°R Required 
°K Required 


in*/sec Required 
Optional 


Optioral 


Optional 


— S TTY erereio - 


Termination of Ballisytcs Data Set. 


Following the last card containing ballistics input 3 card containing the 
words END @F PROGRAM, beginning in card column 1, must be inserted, 
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Takle 2~12 Data Input Parameters (Continued) 


Erosive burning constant 
Erosive burning exponent 
Ambient gas density 

Heat capacity of propellant 


Propellant surface temperature 
(850 for AP-composite, 520 for 
double-base) 


Propellant thermal diffusivity 


Key word for radial slot ig- 
nition, 10 forward faces, 10 
aft faces (value = ignite) 


Table argument value for up to 

10 radial slot ignition diameter 
inputs (key to input table number 
13). Ignited diameter becomes re- 
quired input if slot face is burning. 


Key word to suppress pressure 
rise down the port. Requires 
definition of slots and segments 
to function. 


oo ae le le Pla oS 


2.8 @DK INPUT DATA 


@DK input data is required only in @DK = 1 as specified in the $PR@B 


namelist set. 
2.8.1 SPECIES 


Species used by the computer program are determined in several possi- 
ble ways, depend.ng upon the problem type. Methods used to determine 


chemical species for each problem type are discussed below. 
DK = 1, IREQ = 0 


For @DK problems, species names and concentrations must be input, 


see Section 2.8.1.1 


@DK = 1, IREQ=1 


For @DE frestricted) - @DK problems, the initial start conditions for 
the kinetic expansion are obtained from a restricted equilibrium calculation. 
The species list generated by the equilibrium calculation generally contains 
many more species than the 40 species for which the @DK subprogram is dimen- 
sioned. Therefore a selection process is required to interface the @DE calcu- 
lated oquilibrium start conditions with the ODK kinetic expansion calculations. 
This selection is performed using the following rules: 

hule 1 If a species appears in a reaction, it is selected for the 

kinetic calculation, 


Rule 2 If a species is specified using INERTS directive it is 
selected for the kinetic calculations. 


Kule 3 If any species has a mole fraction greater than an input 
criterion, it is selected for the kinetic calculation. 
Specics which are selected but whicr do not appear in a reaction are treated 


as inert and listed as such on the output list of selected species. 


2-42 


2.8.1.1 F L I NCENTRATIONS 


(APPLIES TO THE @DK ONLY PROBLEM) 


This input begins with a single card with SPECIES in columns 1 through 
7 and with either MASS FRACTI@NS or M@LE FRACTIQNS in columns 9 through 22. 
If the identifier for mass or mole fractions is omitted, mass fractions are assumed. 
Up to 40 species cards may be input. Only those species specified Ly input 


species cards will be considered for an @DK problem. The order of the input 


species cards is independent of the order in which the species appear on the 


master Thermodynamic Data file. 


A chemical species is identified symbolically by 12 alphanumeric char- 
acters and must correspond identically with the species name as it appears on 
the Thermodynamic Data file. A complete list of the current species names are 


listed in Table 2-4. The species symbol may not contain the characters * or =. 


Col Function 
1-10 Not used 
11-22 Species symbol! (left-justified) 
23-30 Not used 
31-60 Value of initial species concentration (if zero 
must be input as 0.0) free field F or E format 
61-80 User Identification if desired 


2.8.2 REACTIONS 


Chemica: reactions must be input for an @DK problem. 


Up to 50 reactions with an implied third body and a total 150 reactions 
may be input to the program. Only one card per reaction, and only one reac- 
tion per card is permitted. Cards specifving third body reactions must precede 
cards specifying all other reactions. Species names appearing in the symbolic 
reaction set mur correspond identically with the species names as they appear 
in the master Thermodynamic Data (see Table 2-4). A card listing for a sample 
reaction set is presented in Table 2-13, 


Table 2-13 


Sample Aluminized Propellant Reaction Set 


[ 
: 
| | 
F ’ 
i i 
k tek 
f oe { 
f REACTIONS 
; Cae 0 8 COP 6 ARE AREDE, Meh, 5 RED R, RAULCH (1OAADLY ' 
¢\? ® Ce © CL « ASP FTV 2M, 9 =O, OKEe SOLNMANI1O71) : 
| wel 2 Wee Ch © ATE RMEFA YT. Perr, © RELAD.37e JACMAS (18h) i 
4 MA? @ oh NP g ABP ANEDR, MEN » H80OG,O « NOe3A AAULOW (1060) 1 
“> e 1h Oo 1g ABPI,NCELP, Nee Khe RED? Ate NOe3 AROK AD CUOTLY, : 
ae Sn HP « ABIAHFIN, Al, CH, AS ALAM NOeS PREHN C1@AT; : 
0? 8 © We OEP, PEFIG, bet Ae ABLIA, Pe NOol JOHNSTON 11 OHA) ; 
or cLe C1 v AlCl? e ABUNFIK Ne Oe8 8 eHE ALA 3 : 
arcureek 8 etl o AEF ITA ewe Oe& Re Me 3? 
ale 9 ® aloe eo Me ®.NFih ene 8 os 0,8 e 6! 
AL ot ALe 4 Ae AAOETA . Ne AS 4 MO ALOe CST IMATE 
Me ALD © ALON , Ae TAF IA , Me AS g AB AL He FSTPMATE 
Od @ ALT © ALNPH « he (‘UAFIA g we ALK g AE Ayfe FSTPMaTE 
ALN ¢ © ALN Ae %, PC )H , Ne F,6 « AH AL He FSTPMATE 
ENN FAR hEAY 
CNM oCLN = CNP © Cl» aey.Fi) o ge, & , Rez, &, ESTIMATE Pd 
hr 2) cae o AZ, 7APIO, Neth, . ABP, 82, MAULCH ()@ARILY 
CH 6 ft B COP © og ABET) » Nem, 5 AML AR, BALCH C}O@KBADLY 
COD © 42 & CO CHD 6 ABILSFA 5 pee 6 ABTS TINMER () 047) ; 
thee 0 8 CM © 92 . Aed,9F93 © Pah, 5» ACHE, ESe RAULCH (1@ARILS 
CL © HAPS CLM © D4 6 ABELE LY . hwe Sh 4 ABz,K, ESTIMATF 
Ch oo HAPe MTL © PF 6 ABH TY 0g NBO 4 ABP, esvy4“arte 
Ci @ HS we MCL 8 4 8 ARE AETV eh, 9 ABbe%e WESTONAERR (AR 
Cr « fe w WL © 4 6 AFILME$D oe Se Le | e MAYFR (y@AT) 
Cie 4 Bre, © DF ¢ ASHLEY © Owe, G% 4g MEIKE, EST 4ate€ 
Cine O # CL © 972 0 ARS FDS rme.t g REY Rb, esSvrsare 
CL> eo MW 8 HEL © CL 6 AFTER LG, Wie, 9 ABA CHERAY (1QAT) 
Cir ¢ O BS CLA © CL o ARSLFYY 6 Atlee eh 9 8889, ESTIMATE 
wee 8 Mes PM 6 The AELME DD 6 bee. S 5 HERA CHERRY CIQATS 
Woe 4 woe MN 6 AMI, TEEL A, weN fh, HE GOLAN NO ROARAS FOTN, 
Wo ft Bo oD g ABPLFM1TI, Meh Ae, HS GL 9K NOo?G BALILCM FIQ@AA, 
Ne NY RB OND © AP 6 AB1,FE13 6 Ne, 6 = BPA Re PAULCH LEOAQIL 6 
Rue Be BE Mo HIM 6 ARK PEED, nen nn, AS H,PMe NOT MAULCHm eH OHA, 
nae 8 Me © Og ABD HOE LO BNL, MBA Me NOIR RELLES F107, 
arm oe ¢f we AL © CF GAB ASE LS eN Bee § oMa VPS « 219 
a, en e aL eno hee, FIT esse,S eAs 2B, F15 e Fan 
a99 6 CL e AL © ALAC L aha, CEP VV e¢ wo. § ehu 3,A271 6 7? 
ry) o (1s z oun eo eheh, 2 bhp eese.% aa Poh? ¢ 2708 
Af ee BLT L a Fam oe rt ehe FEF VY ehigenS eRe 1,79 6 22¢ 
irae & BATH 6 MAL AT AMERY} ot e000 = gke A, 5HK 6 23? 
a, ¢ Cb? es Art o@ @| CAFS, DV bV otweles ehe laa e 2x0 
As o mr after oe bh Ave. oe SO ry) eAvP ART 6 224 
aM oe EL we alt en she UP yet. seAeh eFur, A8V 4 rh] 
; Fh os Os ie ee Oa  ALNEL esa Cab Vem soley ehar RAR ¢ PY 
titan toa Gterr aari ye My Ay See ee oRha} ye e rd ¥ 
eee Cet eee, ee OATI LITE pet veh aha VGPr ¢ ri >) 
er le Oe er eo thd ee 2 ee ee) fo 2) eher?.ard o oe 
or oe ite (Cape ie eek ees Taiat soc oAnm LT Aneto Bohe® etetet e PaF, 
AL Cy ce 4 Bote cite ea tae oe ae ee Cure} eho StH 46 24a 
By 6 we Ce ee er hal PRS} Nees oBe7 Ar 6 26? 
Ts a e AL °F ow ef el POF oo pehes ea? ARI oan 
i ify e ne we Altr, a OAL AEE 1 tive teh eAat, PQA 6 769 
O Pore bLer > wo ALO. ee Or oe en ers 2 OY Par Shs ry Jha 
Ait ao ty? Creve bia) of Aces ohm TARE Sot wees o4elh e We 
RAC een ne re Lal | eho OE Lame ahey eher, KO, 3ya 
On 2 oe Be Veh E YL 9 we ol GS 4g He GC ArAS, ESTE MATIF 
apr ie Cll ew Airs eo OL BB YCARTY «© NO 90,8 5 Aa B,C 1KG, FSTIMATE 
BUCS ee a he fm Trt ht g om 90,6 , Me ERD , ESTIMATE 
ALM o Ie HB AENA © te g As bePS V3 9 he 00,8 » He BeanPh, ESTE MATE 
alle a we wh we ee e ag Vett: i 9 a of & , Ma BHU , FOtyuare 
Cr a es oe fe fete tt « ne of H , He RK ArAG, EST are 
AL TP © mh, F Ce @ Eye tee AG 1, AETT , Oe 0,8 , Aw Bi aVeE, FST Hare 
ALAM © 7 8 Vo FL tite Ae t,0b4% « es 20,8 » Re A, BAA . FSVIMATE 
a ee ee, oe | fe Fett 5 ow OF LR , Hae KL MPAR, ESTATE 
bp tte oe CLI hh © ALT eee 9 AM TA ETT 9 Be OHH 4g Me 1, 81KK, FRI mare 
BL 0 WAS we bln oH, fe «F698 5 owe oF, G , Re FLAS, ESTIMATE 
ALA 6 82 © Alte ow, ep Feth TT ¢ te 9O,% , Be KL TVE » ESTE Mars 
Aber 0 ©" Pp Bw Att 0 FL ef se Leh Ye 6 ow OAL , Me KR ATAS, ESV IMATE 
Alf eo fs we APH eA am J.fFT) 6 8a 00,8 5 he 8 APG. FST MAIC 
2 és QPF YY » oe oU,G , Me 1, 7e e EST ware 
ALw 0 %e @ al ow a ae fork yl 4» tw of, 8 » Me 1, Phe LAS BLT SG 3 
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BLw~ eo fs BA, we 8, og Awe Byte lt 9 te oO, 8 , Oe UL Phe CST mare 
bre REET MHETEL of f 0 HOB, FT ene eh SMe BGAN, Csrivar¢e 
ReproMvailable copy. 
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The symbolic reaction set contains directive cards and reaction/data 


cards in groups as outlined below: 


REACTIONS Directive for start of symbolic reaction 
input 

‘ Reactions with implied third body species 

END TBR RLDAX Directive for end of third body reactions 

* All other reactions 

LAST REAX Directive for end of reactions 

INERTS Specified inert species 

THIRD B@DY REAX RATE RATIGS Directive for start of third body reaction 

rate ratios 
; Third body reaction rate ratios 
LAST CARD Directive for end of REACTIGNS input 


The content and format of each type of card is defined as follows: 


2.8.2.1 The symbolic reaction set begins with a card containing the word 
REACTIG@NS in columns 1 through 9, Other columns on this card can be used 


for comments. 


2.9.2.2 Each card defining a reaction is divided into five fields, separated 


by commas. Each field contains: 


field 1 a reaction 

field 2 A = followed by a value of A 
: c . rate parameters 

field 3 N = followed by the value of N {oF tha reantions 

field 4 B = followed by the value of B, the 


activation energy (Kcal/mole) 
fleld 5 available for comments 
The general form of a reaction is: 
& * ct 83 
Ny Symbol, + N, Symbol, aes N, *Symbo. , + N, b 
where the left hand side represents reactants and the right hand side represents 


*Symbol, +... 


products. The reaction can be either endothermic or exothermic. 
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The multipliers, N, must be integers and represent stoichiometric 
coefficients. If no stoichiometric coefficient is given, the value 1 is assumed. 
The dimensioning currently in the program requires that: 


< 
N, +N, tose S$ 10 


and 


‘ < 
Nt tee $10 


The chemical species (denoted by the word "symbol" in the above gene- 
ral form, can contain up to 12 characters, each of which must match a species 


name contained in the thermodynamic data (see Table 2-2, card 3). Examples: 


Reaccion Interpretation 
NA++CL--NACL Na’ +Cl = NaCl 
R+2+@-2=Bd Bp’ + @ =s8¢d 
BE+2+2*@H-=BEGHGH Be +2@H = Be(@H), 
The value assigned to A, N, B defines the forward (i.e. left to right) reaction 
rate, k, as 
Vay aN : 97 (1000B/R T) 


‘e] 
in units of cc, K, mole, sec. 


All three reaction rate parameters must be input. The numeric value of each 
parameter may be specified in either I, F, or E format. If E format is used, 


the f must appear before the exponent. : 


2.8.2.3 The reac.ions with an implicu third body must precede other types of 


reactions, and must be followed by the directive (columns 1 through 12): 


Col 1 
END TBR REAX 


ull reactions prior to the above directive will have a third body term added to 


each side of the reaction, e.9., 


H2=H +H, eee 
END TBR REAX 


defines the chemical reaction 


H, +M =H+H+M 


where M is a generalized third body. Specific third body effects may be 
included by inputing specific third body reaction rate ratios as outlined in 
2.8.4. Cards encountered after the END TBR REAX directive card do not have 
a third body term added. 


2.8.2.4 All other reactions are input next, exactly as described under 2.8.2.2. 


2.8.2.5 After the last reaction has been defined, a card with LAST REAX in 
columns 1 through 9 is input. 


2.8.2.6 Reaction rate data for 14 dissociation-recombination (implied third 
body) reactions and 59 binary exchange reactions are listed in Table 2-13 for 
propulsion systems containing elements C, N, O, H, Cl ard Al. Cards can 
be abstracted from Table 2-13 for input to the computer orogram. For the 
implied third body reactions, the third bcdy for which the rate applies is indi- 
cated in parenthesis in the comment field (M represents a "generalized" third 
body) . 


2.8.3 INERT SPECIES OPTION 


Inert species (i.c. species not appearing iu reactions) can be included 
in the input by input of a card with INERTS in columns 1 through 6 followed by 


a list of inert species names. The species names must each be followed by a 


comma and each name must be written exactly as in the master thermodynamic 
data, The last comma must be followed by the word END. See Table 2-14 for 
an example. The species list can continue cn to the next card, but & species 
name cannot overl3p onto the next card. 
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2.8.4 THIRD BODY REACTION RATE RATIOS 


As describ2d above in Section 2.8.2.2 for the ee reaction only one 


reaction rate, k,. where 


j 


can be input. For three body recombination reactions such as 
ab = 
H+@H+M ' H,@ +M, 


the rate of reaction is in general different for each species, M i! depending 
upon the efficiency of the species, M,. as a third body collision partner. The 
temperature dependence of a recombination rate {s approximately independent 
of the third body; i.e., for the tee third body and ee reaction: 


° -N, -B,/AT 
7 = Ae je j : 


The third body efficiency of the ris species for the Phas reaction is then defined as 


m. =e 


i /s : 


j 
Thus m is the ratio of the reaction rate with species M, as the third body to 


the reaction rate input on the reaction card described ir Section 2.8.2.2. 


If reaction rate ratios, mie are to be input to. the dissociation- 
recombination reactions, a card with TH 8D B@DY REAX RATE RATIGS in columns 
1 through 27 must be input next. If this card is deleted from the input, the 
program assuves al! Roe lf this card is included in the input, it must be 
followed -ither by 3 card with ALi “QUAL 1.0 in celumns 1] through 13 (which 


sets all my, = 1) or by SPECIES cards as described below. 
The m ij can be input using a card with the word SPECIES in colusns 1 


through 7. Tais vord is followed by the name of the 1°" species followed by u 
comma, followed py the values m, in F format, each followed by a comma. 


These m j values can be continued onto succeeding cards. Note that the m, j 


| 
| 
| 
{ 
| 
| 
| 


TABLE 2-14 LISTING OF SAMPLE REACTIONS CARDS FOR AN 0,/ He PROPELLANT 


REACTIONS O-H MAY 3-4 1972 JANNAF 


PSwG 

He 1H = H20 0 AZT.5F23 « N&Qe6 0 REOee (AR) NO, 

0 ¢# 2 OH @0 AmS,HE{R » NEle .¢ RXOee (AR) NO. «¢ 

nent 2 02 eo Awl e2Ei 7? «© N&Lle 9 BOee (AR) NOe 

HeH # H2 eo Ax6-4E17 » NS&]e 2 B20ee (AR) NO. ' 
ENO TUR REAX Os 

HP? © 1H 2«& He H20 oo ASPeIDEISe Nahe ¢« BP5,18.s RAULCH NO, ? 

OH ¢ OH = 0 © HAN 9 AZSeTSEQPe NRO, 9 BRL TAHS BAULCH NO, ¢ 

HoH 2 0 ¢ HA 9 A&te33E1Z2e NEQe 9 RE7,300» BAULCH NO, ¢ 

0 + 04 2 H-¢ 02 e ABLOSEL3 « N&Qe © BOece . BAULCH NO. ¢ 
LAST REAX 


INFRTS NSOeARe END 

THIRD BODY REAX RATE RATIOS 
SPECIFS ARelooleoleoles 
SPECTES H205¢05e0Senbes 
SPECTES H20—20c eSe0Sen les 
SPECTES N205eeDeebeSodede 
SPECIES Node eheedoolede 
SPECIES MelreSeleeSelleierSes 
SPECTES 0ol%eSei2eSel2cSe2des 
SPECTES NHelZetsieemeleeseeves 
LAST CARO 
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values depend on the order of input of the reaction cards, i.e. the je reaction 
is defined by the en card input after the REACTIONS card, 


Table 2-11 gives a sample input for a hydrogen/oxygen system using 
third body reaction rate ratios. In this example the three body recombination 
rates are input with argon as the third body. The rate with Ho as a third body 
is a factor of five larger than with Ar as a third body for the first three reactions 
and a factor of four larger for the fourth (hydrogen recombination) reaction. 


At this point in the data input deck, a card with LAST CARD in columns 
1 through 9 must be input. 


2.8.5 $Q@DK NAMELIST INPUT 


S@DK Namelist input specifies the conditions for the kinetic expansion 
calculation. The input is read in subroutine @DKINP and consists of the follow- 
ing groups of data as outlined below: 


2.8.5.1 Inlet, Throat, and Expansion Nozzle Geometry 
2.8.5.2 ‘Integration Control 

2.8.5.3 Print Contrel 

~.8.5.4 Species Selection and Mass/Mole Fraction Check 
2.8.5.5 QDK Problem Input 


2.8.5.1 INLET, THROAT, AND EXPANSION NOZZLE GEOMETRY 


The nozzle geometry is defined in Figure 2-2. At the tangent point 
where the nozzle is attached to the throat section either a cone, parabola, 
circular arc, nozzle defined by a well table may be input. Two radi! of curva- 
ture for the throat wall (upstream and downstream) are allowed in ODK (TD2P 
allows only a single radius of curvature). 


Certain of the S@DM Namelist items are communicated from the SGI@M 


Namelist data set. These items will be prefixed with a plus (+) sign. Input 
of any of these data items in the SODK Namelist set will override the valucs 
input in the SGEQ@M Namelist in this calculation module only. 


ANewoey a1ZZON YO 2-2 eNbyy 


(2'o) ~—- (0° 0) 


Sdd 


aucg 


ALMY 
Imo\wuoy 


£5 
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Item Name Input Quantity Units SI Units 


*RSTA? = throat radius inches meters 
tRWTU upstream normalized none none 
wall throat radius of 
curvature” 
*RWTD = downstream normalized none none 
wall throat radius of 
curvature” * 
+e. 
THETAI = inlet angie degrees deyrees 
tRI = normalized inlet wall none none 


radius of curvature 


For a Conical Nozzle Option 


eomtRNg Ee Fee ee ne ~ 


item Name Input Quantit Units 
$ 
; 
+ WALL > ] requests cone option none 
*THETA = cone angle degrees 
TEP = nozzle expansion ratio 


($400). If EP is input, the 
@DK expansion will be 
through EP. 


EPS = set equal to EP none 


Note: Parabola (IWALL=2) and circular erc (IWALL=3) options are available 
only through use of the SGEOM Namelist. 


*The tlansonic analysis requires that a value of RWTU 2 .S be input. RWTU = 
RW if data comes from SGEGM. 


**If a corner expansion (i.:. Prandtl-Meyer fan) is desired, a value of RWTD = 
.U5 is recommendcd. !xperience has shown that values smaller than this give 
the same result but are computationally less efficient. 


For a Contoured Nozzle Option 


Item Name Input Quantity Units 

+IWALL =—4 requests contour option none 

THETA = wall angle at tangent point degrees 

*THE = wall angle at exit degrees 

*NRZS = number of points in none 
table, $20 : 

* PWRS(2) = normalized radial wa:] none 


coordinate table* 
* pwzs(2) = normalized axial wail none 
coordinate table* ; 


RZNYRM = optional normalization none 
factor fcr input wall 
coordinate tables 


2.8.5.2 INTEGRATION CONTROL 


The program begins its calculations using an implicit integration 
method to integrate the fluid dynamic and chemical relaxation equations. 
For near equilibrium flow an implicit method is used since it is inherently 


stable. However, once the flow has become sufficient frozen, explicit 


basta aetna aes Gaal eee Te 


methods become oumerically stable and can be used more efficiencly. The 
program uses temperature as a freeze cnterion in order to switch from an 
implicit to an explicit method. This criterion is T < TFEXPLI where TEXPLI : 
is assumed O°R but may be input, The program always prints the axial posi- 

tion at which the switch occurs, For high area ratio nozzles this procedure | 


can be especialiy useful, 


The integration routine controls the step size such that the relative 
error in the dependent variable increments are less than a prescribed fraction, 
DEL. Only doubling or halving of the step size is permitted, and on option, 
either all the variables may be considered (JF=0), or oniy the fluid dynamic 


variables (JF=1) may be considered, 


| 
* PWRS(1) and PW2S(1) are internally computed at the coordinates of the contour 
attachment point. 


When the flow becomes supersonic and the area defined fluid dynamic 


equations are used, an additional check on continuity is applied in the form 


(PVA) 41 7 PVA). 


(pvA) < C@NDEL 


N+l 
where C@NDEL is an input relative criterion. 


The step size is held between the two input bounds HMIN and HMAX. 


Fixed step cases may be run by setting input values for HI, HMAX, HMIN all 
equal. 


Assumed 

Item Name Input Quantity Units Value(s) 
HI = initial step size none 01 
HMAX = upper bound on step size none 0.10001 
HMIN = lower bound on step size none 005 
IBYSIE = fractional incremental none 001 

error 

° 

TEXPLI = temperature below which R - 

explicit integration will 

start 
C@N DEL = relative error criterion for none 1x107° 

continuity check for super- 

scenic flow 
pe + Q all variables considered none 0 

for step size control 

= ] only fluid dynamic variables none 
considered for step size 
control 


2.8.5.3 PRINT CONTROL 


Output from the Kinetic Expansion Calculation consists of complete 
output for each print station selected, The end point of the nozzle is always 


printed. Print stations are selected from one of the following options: 
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Assumed 


Item Name Function Value (s) 
JPRNT = -2 print throat and input area ratios = 
= -] print at internally set area ratios for conical 


nozzle.* Print at input wall contour points 
fer contoured nozzles 


= 0 print at every integration step 


+] print every ND3rd step up to the throat and 
then nominal area ratios 


= +2 print every ND3rd step over entire nozzle 


If JPRNT is +1 or +2, the following must be input: 


Item Name Function 

NDl = first integration step to be selected for ~ 
print 

ND2 = last integration step to be selected far - 
print 

ND3 = print every ND3rd step between NDI and - 
ND2. 


If JPRNT is -2, the following must be input: 


Item Name Function 
*ARPRNT(1) = requested area ratics for print, must be = 


monctonic increasing and greater then 1.6 
(usually entries are the same as those 
used in SUPAR of S@DE). 


number of area ratios requested for print - 
$100, 


*NJPRNT 


An extended print option may be selected as follows: 


Item Name Value Function 
IDYSCI = 0 no extended print requested 0 
= extended print option selected (aot suggested) 
*For JPRNT =-1 and a conical nozzle (i.e. IWALL = 1), the intemally set area 
ratios are: 
ARPRNT(1) = 2,364,0+-,39,40,42,...,58,50,64,...136,120,128,..., 200,216, 
220,..-.400,412,444. 


2.8.5.4 SPECiES SELECTION AND MOLE/MASS FRACTION CHECK 


In order to interface @DE equilibrium calculated start conditions with a 
the kinetic expansion calculations, special considera.ion must be made for 
inert species (those not appearing in the reaction set). Inerts may be selected 
explicitly by use of the INERTS directive or by use of a relative selection 
criterion. The INERTS directive is described in Section 2.8.2.5. 


The relative selection criterion is described below: 
item Name Function 


EPSEL = all species which do not appear explicitly 
in the reaction set but whose mole fractions 
are greater than the input value for EPSEL, 
will be retained for the kinetic expansion. 
Species selected under this criterion are 
treated as inert. The program assumes 
EPSEL = 1,0E-5, unless input. 


In some instances it may be desirable to use input species concentra- 
tions which do not sum to unity. Species concentrations, either input or from 
equilibrium start conditions, are summed and the sum checked as described 


helow: 
Item Name Function 
XMFTST Ze input species concentrations are summed 


and checked versus unity using this input 
criterion, If 


[1 - species concentrations| <XMFTST 


then the test is passed, The species con- 
centrations will then be normalized such that 


3B species concentrations = l. 
".he program assumes XMFTST = 1.0£-3, unless input. 


If the test is not passed, an error message wil be given and the run 


terminated. 
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2.8.5.5 @DK PROBLEM INPUT 


This input is required when an @DK calculation is requested without 
restricted equilibrium start conditions being requested. A kinetic expansion 
from input arbitrary start conditions is to be computed. In addition to the input 
items described in Section 2.8.5, an @DK problem requires input of those items 
described in Sactions 2.8.1 and 2.8.2. 


Item Name Input Quantity Units 
PC = chamber pressure PSIA 
T = initial temperature R 
V = initial gas velocity ft/sec 
JPFLAG = 0 pressure table calculatea none 
internally 

= ] pressure table input none 

ECRAT = initial contraction ratio none 


For JPFLAG = 0 option, the following must be input: 


Item Name Input Quantity Units 
PI = initial pressure PSIA 
PESTAR = throat pressure PSIA 


For JPFLAG = 1 option, the following must be input: 


Item Name Input Quantit Units 

PTB(1) = normalized pressure table none | 
entries* 

ZTB(1) S normalized pressure table none 
coordinates** | 

NTB = number of pressure table none | 
entries, £127 | 

Z = initial axial position none | 


*normalized to input chamber pressure, PC 
**normalized to input throat radius, RSTAR 


"CORA MARINER nel Bandits nt ses aye = 
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2.9 Two-D:mensional-Two Phase Flow (TD2P) Module Input 


Data required to execute the TD2P module is read in under the namelist 
name $TD2, Certain of the input items described below are either communicated to 
the TD2P module by the GE@M, BAL,and/or @DE modules or are preset in the TD2P 
module. Hence, these data items do not necessarily need to be input in the $TD2 
namelist set. However, any values that are input will override the communicated 


or preset values. 


The input data items to the TD2P module are as follows: 


Item Description Units cen 
$TD2 namelist name 
The Sites items are communicated from the @DE or BAL module 
PR = Sar ance none 0.0 
GMGY = chamber gas viscosity coefficient Ibm/ft.sec. - 
CAPN = viscosity temperature exponent none - 
CPL = liquid particle heat capacity (T>T,, ) ft? /sec® °R = 
TPM = pe particle solidification Superdiice °R - 
DHM = latent heat of melting for the particles ft°/sec” - 
EPM = expansion ratio at which solidification uf none = 
the particles begins in an @DE type calcu- 
tion 
UGM = gas velocity at the @DE solidification ft/sec 7 
point 
TPS = gas temperature at some point after soli-  °R 2 
dification has been completed. 
UGS = gas velocity corresponding to the value ft/sec - 
of TPS 
PC 3 chamber pressure psia - 
TGO = chamber temperature °R = 
WPWGT = ratio of particle to gas weight flow none > 


XMLW 2 molecular weight of the condensed phase 1b/lb-mole 101.96 


On the follcowing geometric input data, the items preceeded by 
an asterix (*) will have been communicated by the GE@M or BAL 
module if they were input in the SGE@M namelist set or calcula- 
ted by the ballistics (BAL) module. 


= te i 8 aan yet ata 


Inlet and Throat Parameters (see Figure 2-3) 


Item Description Units Assumed 
Value(s) 
DZI = Az, particle trajectory integration step none 002 
size 
DZMIN = AZ in! inlet step size parameter none 002 
NILP = Ni. number of initial line points none 20 
*RRT = R_, throat radius of curvature. A value none - 
Re 1 ts required 
*RT = r*, throat radius ft “ 
SAUR(1) = first estimates of x_, u_, a, 8, andy none -.15,1,.5, 
for the special throSt expansion. Re- 3,71 
quired only if 9, > 6, 
THFD = 9,, faring angle (6,> 6, = >No faring) degrees 5. 
*THID = A, inlet angle degrees 
THIW = 0, , intersection of initial line and wall degrees 12, 
WwW 
THJD = A,, angle defining the zone farthest dov'n- degrees 9, 
stream 
VAR (1) = first estimates of x_, u_, a, 8, andyfor none .3,0,0, 
the zone farthest upstream Be one | 
LAX = 2 avis? intersection of initial line and axis none - 
rA\ 3 Nie number of upstream zones none a6 
2] ee n,, number of downstream Zones none os 
NTBL = number of points in the subsonic-tran-~ none 1S 
sonic region to be output for Input to 
the TBL module 
XITBL z starting normalized axial position for the none = 
Output to TBL 
*XLSTAR = motor L*, if input in $GEQM, it will be in - 


used unless it is calculated In the Bal- 
listics modu'ec, In which case, the 
average L* from the ballistics module 
wil) be used. A value input here will 
over-ride previous values, 


Characteristics Mesh Contro! Data 


DL = A£, maximum LRC mesh width none v4 

DTWI = AA , maximum flow angle change along degrees a 
thé wall 

DR = Ar, maximum RRC mesh width none 2 
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Item Description Units Assumed 


Value(s) 

EW = ew end of nozzle wall criterion none .001 
IMAX = , maximum number of iterations per none 5 

“TRA point 
Nl z= Ny, select each mS LRC for print none 1000 
N2 = Nos print each ny’ point on selected none l 

Nozzle Wall Contour Data 

*TWALL = Option flag none = 

0 = >tabular input (not allowed in $GE@M) 

1 = >cone 


2 = >parabola 
3 = >circular arc 
4 = >wall contour (spline) option 
If the wall is to be input in tabular form (IWALL=0): 
PW(1) = (r,, z,), wall coordinates none « 
{21,2, .., n points, viz.: 
PW(i) = Pye 2ye Foe Zoe coos Cys Zoe 0,0 
Note: 


3) always mark the end of the table with 
two zeros 


b) ns74 !s required 
c) always setr ape rcee THIW) and 


ar alae itt Ta aii shalt abet enim tale tuned 


2. :BRRT stn 
®EPS ss E, Rosle expansion ratio none ° 
Ifa cone, parabola, or circular arc contour {s to be specified 
then: 
*THIW = ein attachment ongle for the contour: degrees = 
e.g. for a cone, the conical half angle 
EPS a c, nozzle expansion ratio (cone only). none - 
*RWMAX = , nozzle exit radius (parabola or arc = none ° 
NN. 
SZWMAX 2 nozzle length from throat to exit none “ 


(Btibota or arc only). 


Description Units Assumed 
Value(s) 


For a Contoured Nozzle 


* THE = wall angle at exit degrees 


[ eed 
et 
© 
=| 


*NRZS = number of points in table, < 20 none 
*PWRS(2) = normalized radial wall coordinate table” none 
*PWZS(2) = normalized axial wall coordinate table” none 


Particle Data 


iT@T = nuinber of particle groups to be considered none 5 
SMP = particle density Ibm/ft” 250. 
R(1) = - , the radius of each of n particles is ft 0.0 


to!be input so that oS Sf «eeSf 


set r =0,n< 10 ts required. 7 
Pratl 


WPWT(1) w Pave, , particle weight flow fractions none 
) 


oH 


corresponding to each of the above particle 

radii, r 
Py 

Note: The mean particle diameter end distribution will be calculated from the 


following data if R(1) = 0.0, The equation 


d= XK - poPEXP | (ey XTEXP (1 _pxp(xteXLSTAR)) 


*(1424°XD-r*) 
ts uscc to calculate the mean particle diameter in microns, 
XK = constant in the d,, calculation microns 0.454 
PSI 
PEXP s constant fu che q, calculation none 0339334 
AL = constant in the d., calculation in”! -.004 
XD = constant in the qd, cuiculation inv! 0045 
XI EXP = constént in the qd, calculation none 033334 


In the above formule: ¢ is the concentration of condensed phase (moles/100 gr products) | 
calculated from WPWGT and, PC is in psia, r* is the throat radius in feet. 


SIG = geometric standard deviation used in none 1.9 
the particle size distribution 
SEND end of $TD2 namelist set 


*PWRS(1) and PWZS(1) are internally computed at the coordinates of the contour 
atta hmecnt point. 
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2.10 Turbulent Boundary Layer Module Input 


Data required to execute the Turbulent Boundary Layer (TBL) Module is 
read in under the namelist name STBL. Certain of the input items described be- 
low are either communicated to the TBL module by the @DE and/or TD2P modules 
or are preset in the TBL module. Hence these data items do not necessarily need 
to be input to the TBL module. However, any value that is read in will override 
the assumed or ..cinmunicated value. 


The input data items to the TBL module are as follows, 


Item Description Units Assumed 
Value(s) 
STBL naemelist name 


The following items are communicated from the WDE module if 
@DE = 1 or 2. 


T0 2 Stagnation temperature (from @DE-chamb- °R* = 
er temperature) 

PO = Stagnation pressure (from @DE-chamber psf* - 
pressure) 

PR = Prandt! number none ~ 

2M U0 ce Value of viscosity, u, at the stagna- lbm/ft-sec,* - 
tion temperature 

ZMVIS = Exponent in the viscosity-temperature none - 
equation 


we 2Mu0*(T/T0)2MViS 


The following items are communicated from the TD2P module 
if TD2P = ) or 2. 


PO = Stagnation pressure (from TD2P-the value psf* = 
of stagnation pressure at the initial line) 

GAM0 = Stagnation value of specific heat ratio none = 

RBAR Xs Value of gas constant f/sec°R* - 

IXTAB = Number of points in the x, y and Mach none w 
number tables 

SCALE = If it is desired, non-dimensional values FT* 1.0 


of x and y may be input and then multi- 
plied by a single scale factor. 
TD2P, SCALE=r*) 


Stagnation temperature (from TD2P the ORe = 
value of the stagnation temperature at 
the initial line). 


= 
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Item Description Units Assumed 
Value(s) 
WDYT = Mass flow rate lbm/sec ~ 
XITAB(1) = IXTAB values of x, axial distance, in none* - 
monotonically increasing order (from 
TD2P-XITAB()=x,/r*) 
YITAB(1) =  IXTAB values of y, nozzle radius or con- none* - 
tour height in consecutive order corres- 
ponding to each x of XITAB. (from TD2P- 
YITABQ)=r /x*) 
ZM TAB(1) =  IXTAB valves of free stream Mach number none o 
in consecutive order corresponding tv the 
contour points Xie Yi 
The following items are preset or must be input to the TBL 
module, 
lAZETA = Velocity pewer law exponent, n none 7 
IPRINT = 0 Prints omy ot input intervals none 0 
1 Prints at al] calculated points 
ICTAB = 0 Indicates a constant specific heat none 0 
calculation 
= 23 Number of voints in C_ versus 
7 table P 
ITW TAR = <1 Le = adiebatric wall temperature none |} 
0 Ty = constant (see TWTAB) 
1 IXTAP values of Twill be input 
(see TWTAB) 
2NSTAN = Interaction exponent n, in Stanton number none od 
relation 
DXMAX = Maximum Jenath of step size, If DXMAX ft* 01x 
is not input DXMAX is computed to be 
1/300 the length of the contour 
EPSZ = Geometry indicator none 1.0 
0. Two-dimensional Planar flow 
1.0 Axisymmetric flow 
FJ = Conversion tactor betwcen thermal and ft-lbf* 778.2 
and work units BTU 
G Proportionality constant in equation F= ft-lbm* 32.174 
(m/g)a sec” bf 
TCTAB(1) =  ICTAB values of temperature in incres- *R* 


sing order, defining the C.-T table. 
The maximum T should exceed T, and 
the minimum T should be below fhe low- 
est temperature likely to be encountered 
in the system 


Ss f2 
= VE 


Item Description Units Assumed 


ams Value(s) 
*cPTAB(1) =  IGTAB values of C, in consecutive order  BTU*/ : 
corresponding to the values in the temp- _ lbm°R if 
erature tabie. [1 
TWTAB(1) = If ITWTAB = -1 omit this input °R - 
If ITWTAB = OQ input one value of wall i | 
temperature if 
If ITWTAB = 1 input LXTAB values of wall 
temperature corresponding | 
to the wall locations of the a 
x and y tables 
T@LCFA = Tolerance in Cg-C¢Rg interaction loop none .0001 
TOLZET = Tolerance in zeta interaction loop none .0003 
TOLZME = Tolerance used in gas property eval- none 107? if 
uation loops 
THETAI = Initial value of momentum thickness. If = 


THETAI is input as a negative value, the 
sonic point start procedure will be actuated. 


PHIT = Initial value of energy thickness, This - 
need not be input if the sonic point start 
procedure is used 


a ee ee ee 


Edge condition option flag 


IFEDGE = Q edge conditions are computed from the none 0 
Mach number table 


= i ‘ablee of P-/P., T_/T_., GC. |, V..0 
are to be eee, Feldw Pe © © 


> 2 tables of P_, T., GC. ; ve 0, are to 


NERA NE NPR ELENA PARR ED McG a RABE 


be input as belbw Pe 
PITAB(1) =  {XTAB values of P_/P. or P_ versus X (sce = 
IFEDGE) pe a ae 
TITAB(1) =  IXTAB values of T./T_ or T. versus X(see - t 
IFEDGE) Soe ae 
*CPITAB(1) = IXTAB values of C,_ versus X BrU/Abm* . 
VITAB(1) =  IXTAB values of V e versus X ft/sec* = 
R@QITAB(1) = IXTAB values of s 2 versus X lbm/fte* “ 
SEND 


TNOTE-Any consistant set of units in accordance with the values of FJ and G may be 
used. The units shown are those consistant with the values communicated from the 
@QDE and TD2P modules end by the TBL module back to the SPP control routines. 


*NOTE-If a C_ table is not input Into the program (see CPTAB and ICTAB), the pro- 
gram vill attempt to edit in a C_ table from the values of C e V5. Ta. However, 
the edited values mzy not (espécially in the cold wal! case) cover the entire range 
of values required fer the “un, The safest procedure is to generate by hand a CPTAB 
table which includes several values of Cy for both the upper and lower temperature 
extremes, 
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2.11 Linkege 


The linkage between the five basic computational modules in the SPP com- 
puter progr>m was designed to give the user of the code flexibility in selecting 
which of the modules need be used in any one computer run. While all of the 
modules may be used in ary one computer run, this is not usually the most efficient 
approac!:, since input mistakes or the use of the program in parametric studies 
could resuire the repeated successful execution of several of the modules. Thus, 
for the sake of efficiency, and overall flexibility, all of the linkage data between 
modules (except for the restricted equilibrium interface between the @DE and 
ZK modules) is written out on logical unit 8 as formated data. This data may then 
be manipulated in a number of fashions depending on the computer hardware and 
utility routines which are available. On 6000 series CDC equipment which the 
program was developed for, logical unit 8 was equivalenced to the formated punch 
file, Hence, at the completion of 2 run, the linkage data would be punched out 
anc made available for subsequent runs (described below), eliminating the need for 
repeatiag an already successful! solution. 


This type of linkage has two major advantages. The first is obviously a 
savings in computer run time on series of cases wherein the output from one or 
more motules is unchanged, The sccunc is the ability to use other analyses (if 
appropriate) to “fool” the SPP code into accepting them as an integral part of its 
own colputational procedure. This second approach allows the engineering user 
the flexibility to substitute the outnut from any other code or analysis Into the 
procecure defined in the SPP coae as long as it conforms ic the proper format. 


Data is read in on subsequent runs by the SPP program on logical unit 3 
if the salculation module cquals 2 option [s selected in the SPR@B namelist. 
However, this logical unit may be changed to any value by input in the SPROB 
namellst. This feature allows the puncked output linkage data to be read in as 
part ci the norma] card input stream as long as the linkage data is in the follow- 
ing order: 


SPROB 
module option = 2, 


SEND 
@DE data* 
BAL data* 
@DK data* 
TD2P = data* 
TBL data* 


the rest of the input data 

*this data is only needed if the module 

option =2 is selected in the SPR@B namelist. 
The data does not need to conform to the above order a: long as the data resides 
on logical unit 3 or any other logical unit which is not assumed _ to be the card 
input stream, 

Certain complications can arise when previously calculated (punched) 

data is used out of the ordinary calculation sequence presented in Table 2-1 
of Volume III or Figure 2-2 of Volume Ul, For example, if the @DE module is ex- 
ecuted while the ballistics module (BAL) data is read in from cards, the chamber 
pressure from the @DE module is used in the TD2P module instead of the average 
pressure calculated by the BAL calculation, While nominally small, the difference 
between these two pressures can cause slight differences ir the performance re-~ 
Sults caiculated iy these modules, In any event, this deficiency v-ill be cor- 
rected in any further versions of the SPP code, 


siege) 


The punchec (or formated output) from cach of the computational modules ’ 
conforms to the seme basic format. The first card image contains the name of | 
the moduic,while the last image cuntuins the word END, The data between these 
curds consists of sequenced “packets” of information which contain the Jinkage 
data. A packet may consist of one or moire data cards, cl! of which have the 
same sequence number, Table 2-15 contain. the basic format used in generating 
the linkage data. Tables 2-16 thraugh 2-20 uescribe th specific data transmit- 
ted between modules, 


2-67 


Table 2-15 Basic Format of Linkage Data 


first | Ad name of module 


al! others 5X, 5E14.6, i5 data and packet sequence number 
last A3 the word END 


Note that in the following five (5) tables the first and last cards have been 
omitted and that the five flelds refered to are the five floating point data fields 
described in Table 2-15. 


Table 2-16 @DE Linkage Data 


1 l Kreg ~ reference viscosity - Ibm/ft. sec. 
2 w ~ viscosity exponent in u=y og (T/T,)” 
3 T, - chamber temperature - °R 
4 P ae chamber pressure - atm, 
5 , ~ chamber viscosity - lbf > sec/ft? 
2 1 %, ~ chamber gas thermal conductivity - lbf/sec-°R 
2 Yea chamber equilibrium ratio of specific heats 
3 8 - erosion parameter 
4 Ch ~ theoretical c*=P_A*/m - ft/sec 
5 Cie chamber equilibrium sp2cific heat - BTU/lbm °R 
Cc 
ths ] €melt 2°68 ratio at solidification 
2 Unelt ~ Velocity at solidification - ft/sec 
3 ia elt ~ temperature at solidification - °R 
4 a - fraction of condensed phase to gas phase 
| 5 AH ee > heat of melting of condensed phase (ft?/sec*/lbm) 
| 4* ) “py - specific heat of the liquid - (ft?’sec?)/lbm °R 
2 nts specific heat of the solid - Uft?/sec*)Albm °K 
3 
3 M.. > molecular weight of the condensed phase 
CP 
4 Toit - temperature at the exit plane - °R 
5 este velocity at the exit plane - {t/sec 
5 ! M,. ~ molecular weight Jf the gas only 
Gas 
2 M.. - ino.e: thar weight of the mix 
Yn ix 
a M,, - effectiv: molecular wight of the mixture 
eff 
4 Pr. - chamber Prunctl number 
‘sb - theoretical ley at the exit plane 
@ fee 
6 l ‘spe frozen ep at the exit plane 


*Not used if no condensed phase !8 present. 
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Table 2-17 Ballistics Linkage Data 


Packet Card in Field Description 
"Packet" 


t - time - sec 
Po - chamber pressure - PSIA 

oe gia maximum pressure to this time - PSIA 
R* - throat r-dius - in 


L* - motor L*® - in 


na ~ fractional c* efficiency 
m - mass flux - slugs/sec 

m fn - ratio of insulation to propellant mass flow 
P_. - head end chamber rressure - PSIA 


Cl 
fmdt - mass expended to this time - slugs 


F - thrust - Ibf 


* not currently used 


Table 2-18 @DK Linkage Data 


Description 


| Packet | Field | 
laN* €pp - area ratio corresponding to I sp in field 2 
I Soar - "Restricted" equilibrium I 5p 
€pp ~ area ratio corresponding to 5p in field 4 
I. Pre - “Restricted” equilibrium Wee 
N+]+END €, ~ area ratio corresponding to I, p in fleld 2 
Isp - “kinetic” iF as calculated by @DK 


*The packets ]+N are terminated by a negative number in field 3 (if the number 
of restricted equilibrium points is odd) or, by the start of the @DK data (if the 
number of restricted equilibrium points is even). If no restricted equilibrium 
points are available, the kinetic eificiancy will be calculated a8 y,7=1.C. 
The restricted equilibrium pairs need not be tn order as this table is sorted. 
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Table 2-19 TD2P Linkage Data 


"Packet" Description 


1 1 Cy ~ discharge coefficient 
2 R = gas constant, ft?/sec*/R 
3 + - specific heat ratio as calculated by TD2P 
4 C. - specific heat, ft?/sec?/R 
5 m - mass flux, lbm/sec (for single nozzle in mul- 
tiple nozzle cases) 
2 l P, ~ stagnation pressure at the TD2P wall initial 
point, psi 
2 T ne stagnation temperature at the TD2P wall ini- 
tial line point, °R 
34END ] t/t* - r.ormalized nozzle wall radius 
2 20" ~ normalized axial] distance 
2 M, - Mach number on wall stream) ine 
4 i Bees 


Table 2-29 IBL Linkage Data 


Description 


r- wall radius - ft 


‘2+ TL decrement due to turbulent boundary iayer 


Pan axisymmetric nozzle hos been as. umed here since output vill be for A/A* 
ener ps ly where AZA®  (pcr4¥ 


~ 
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3. OUTPUT DESCRIPTION 


As a result of the many types of analyses performed by the computer pro- 
gram considerable output can be generated, Since the amount, and nature, of the 
output will vary over a wide range, as a function of the options selected, the out- 
put descriptions will be keyed to the individual overlays types of analyses). The 
sample cace output in section 5 is a valuable adjunct to the output description, 


and should be jiberally referred to. 


3,1 Overlay O ram control ar.d summary) Output 


Since all of the program Input cards do not appear in formated output 


statements, it is recommended that a card image listing of the input cards be 


obtained, via control cards, at the initiation of the ru. 


Pe eR RUMR Geren prep yAeriign ox Ae laM 


The first output from the program is the case title. This is followed by 
an acknowledgement of the program directives that have been considered (GEOM, 
TRAJECTORY, etc.). After acknowledging the TRAJECTORY directive the trajectory 
information is cutput in tabular form. 


Following the PROBLEM directive, the program cnters tnto execution of the 
various cverlays as directed by the SPROB input. The output produced during the 
exccution ot these overlays is discussed in subsequent tubsections, 


Upon compietiun af each of the overlays from -.01> 5,0, the following 


outputs ure produced. 


OVERLAY 2. If BaL = l or 2, the calculated values of P_ (psi), L* (in.), ©* (in.). 
de® “dt (in.’sec), Miclugs/cec), and». are ortutes os + function of time in see- 
i. 


ands, «‘:ble of the averajye, moximur. and mins cum ufe ack of the above quanti- 


tles is 216° vatput, The sveraye values oc cam teu tare. on the total caleu- 


lated bum tima urtess the following messaye is printed out below the table: 


VALVES AT TIMES GREAISR TRAN miAVE BEEN EXCLUDED 
PROM [ab AVERAGE CALCULATION 


The excluded pol:ts (snould there be any) represent points during tail-off that might 
produce spurious », valucs that unduly affect the calculate averege » ce Value. 


OVERLAY 4, If ODK= 1, or 2, a table of the one dimensional kinetic specific 
impulse, ISP Op x: and the kinetic ISP efficiency factor, Mine ae printed out a8 
a function of nozzle area ratio. 


OVERLAY 5. If TD2P= 1, or 2, a table of the calculated two dimensional, two- 
phase, perfect gas, specific impulse, ISPonop- and the <-D, two phase, I 


Sp 
efficiency factor, 7p2P' is printed out as a function of nozzle area ratio. 


If RSD@T#0 in $ the SGEOM namelist, or if Overlay 3 (Ballistics program) 
has been executed and indicates that the throat erodes, then the nozzle expansion 
ratio will vary as a function of time. 


Ir, such cases Mop op and ! SPeaop are time dependent. The program cal- 
TD2 
culates and prints out the time averaged values of these quantities. 


The program also calculates and cutputs a value of "rp2P to be used in 
the overall! efficiency calculation. When there is throat erosion the use of this 
modified efficiency corrects the predicted overall delivered specific impulse for 
the fact that We varies with expansion ratio (see subroutine SUMRY), 


QVERLAY 5. If TBL= 1, oc 2, a table of the incremental loss in dee Al 5p ; 
is output as ¢ function of nozzle area ratio, T6L 


After all of the individual loss calculsutions ‘uve oven conpleted a sum- 
mary page is printed out, Figure 3-i reproduces a tyr7ic.+!) output summary page. 
A series «f deseriptive notes are given bclow, Trese cin be correlated with the 
cucled numbers on Fieure 3-1, 


Note: 1. These values are caiculoted when the relevant parameter in the SPROB 


i 
w 


“a 


Ld 


Input equals io 


éo The calculoted 2-0, two phase efinciencty alsr tacludes the so-called 
erosion loss, The magnitude of the cr -ion joss {in terms of an efficiency) is 
ual to th ra | octet to, thimate an) to 9% ” 
equal to th 10010 Of Nery, (corr at erosion) to 'rpop (tine averaged: 
vee summary output from TD2P module). 


3, A “caleuleted" vaiue of combustion efficiency 1s obtained only if 
Bal. = !, or 2 and TD2P 2 1, of 2 in the SPROB namelist. The “caleulated" com- 
dustion efficiency is equal to the nozzle discherge coefficient (calculated by 
TDP) times the “empirical” combustion efficiency. 
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4. This efficiency cannot be theoretically calculated with the present 
program. 


5S. Tre calculated 2D 2 phase efficiency includes this effect. 


6. Empirical efficiency factors are always calculated, (Except for 
empirical combustion efficiency which is calculated only when BAL = 1, or 2, 
and TD2P © 1, or 2). 


7. Includes the two-phase efficiency only, If BAL #1, or 2 this efficiency 
will be calculated only if AVELS in SGE@M Namelist is #9. 


8. If the real nozzle s;eometry has been altered as a result of the restric- 
tions imposed by the TD2P program, the output value of this efficiency will be 
erroneous when the nozzle attachment angie (supersonic) implied by the input 
nozzle contour differs from the actual nozzle attachment angle. In such cases the 
correct value may easily be calculated by hand, 


9. These values are equal to unity unless the reievant SPROB parameter 
equals 3, and a non-unity velue has been inrut. 


10. For each of the Individual loss mechanisms, this column give: the 
efficiency factor selected for use in the overall efficiency calculation. The 
selection criteria is as follows, in order of precedence: 

1. Input values. 
2. Calculated vaiues. 
3, Empirical valucs. 


11. The product of the selected efficiencies is not, in general, equal to 
the product of the efficiencies calculated by any one of the methods, as it may 
eentain efficiencies calculated by cifferent methods, 


12, This cozumn is 1 key to the ethod selested for use in the overall 
efficiency calculation, The meanings of the ietters arc as follows: 
! - Input 
C - Caleulated 
E - Empirical 
N - Not considered (a value of unity 
is used In such cases) 
13. The values of delivered vacuum specific impulse are equal to the 
theoretical specific impulre at the initial (t=0) nozzle expansion ratio multiplied 
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by the appropriate product of efficiencies, and then decremented by the corres- 
ponding turbulent boundary layer loss. 


14, The empirical delivered vacuum specific impulse does not include the 
so-called erosion loss-see note 2 for how to include this loss in the empirical 
result. 


15. The average ambient pressure is the average of all points in the tra- 
jectory table, The average chamber pressure is obtained from the summary table 
following the ballistics calculation (or is input {f a ballistics calculation is not 
called for). The average expansion ratio is an average based on the whole calcu- 
lated burn time. The discharge coefficient is obtained from the TD2F program, 


If BAL = 1, or 2, and TD2P=i or 2, in the $PROB Namelist, = page entitled 
“Thrust - Time History” is output at the end of a run. The table is basically self- 
explanatory, however the following items should be noted: 

1. The mass flow (MDOT) and the integrated mass flow (MDOT-INT), 

include the corresponding total mass flows (i.e. both propellant and in- 

sulation). 


2. Due to inaccuracies in the input grain geometry the total integrated 
mass flow may not equal the known propellart (and burned insulation) 
weight. In such cases the calculated total impulse (F-INT) should be ad- 
justed by the ratio of calculated to known propellant weights. 


3. The average I. ts the average delivered Isp ai.d will in general differ 
from the "I. delivered to the average ambient pressure" which appears on 
the aforementioned summary page. 


4, The average MDOT may d'!“fer from the corresponding value in the Bal- 
listics Ave,~Max.-Min.. T2ble as a result of different time bases for the 
two averages. 


5. Negative I-DEL values appearing during tailoff do not unduly affect the 
calculated average, AVE. ISP = F-INT/MDOT-INT. 


2 Overlay z (ODE) O t 


The first Jine of output says that ote zone is being calculated (multi-zone 
calculations are not permitted in the present application). The input reactants 


cards are then reproduced, followed by the INSERT/OMIT cards (when applicable). 
The NAMELIST input is not printed out in the program's present form. The word 
NAMELISTS and the print out "NO SODE VALUE GIVEN. . ." should be neglected, 


A list of the species being considered in the calculation is then printed, 
Each species in the list is preceded by some identification such as J12/65. The 
J refers to JANNAr data! 8) The letter L refers to unpublished data calculated at 
the NASA Lewis Research Center. The number refers to the month and the year the 
data was published or calculated. 


Following the list of species is the current value of O/F. Thisis follow- 
ed by a listing of the enthalpies or internal energies of the total fuel and oxidant, 
and of the total reactant. Following this ts a list of the kilogram-atom per kilo- 
gram of each element in the total fuel and oxidant, and in the total reactant, 


The next output yields information regarding the iteration procedure used 
in obtaining equilibrium solutions. User's interested in interpreting this information 
should refer to Reference 1 . 


The primary output of this program is basically self-explanatory. If 
EQL = . TRUE, in the SODE trvut the results of an equilibrium rocket performance 
calculation will be printed out. This output contains the following: 
1, Heading 
. WChamoper pressure 


. Proportion of oxidant to fuel 


2 
3 
4, Thermodynamic mixture properties anc curivatives 
5. Rocket performance informetion 

6 


- Composition of the combustion products 


Items 4, 5, ani are cutput ut conditions corresponding to the chamber, 
throat and requested area ratios (SUPAR and SUBAR arrays). They are also output 
(first exit column) at the arca ratlo cz which particle solidification begins (if 
applicable). 


The aforementioned tables are followed by a listing of the trace species 
Which were considered In the ~alculation, but which do ot appear in the mole 
fraction table, 


A table of the mass fraction of cordensibles is then printed out. These 
values correspond to the same relative locations given in the main table, 
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The molecular weights of the condensibles and gas only, in the chamber , 
are printed. The value of the erosion parameter, 8, used in the Ballistics Module 
throat erosion correlations, is then indicated. 


The next page of output gives the results of the transport property calcula- 
tions. Viscosity, MU, frozen thermal conductivity, K, and frozen Prandtl number, 
PR are given for the chamber and at the throat and exit plane. The power law vis- 
cosity parameters calculated for use in the TD2P program are then printed, followed 
by a list of the species considered in the transport property calculations. 


If FROZ = . TRUE. in the $ODE input a second series of output is generated 
for the results of a frozen rocket performance calculation. The output is identical 
to that described for the equilibrium solution, with the following exceptions: 

1, There cre no table entries at the area ratio corresponding to particle 

solidification. 

2. Mole fractions are frozen and hence do nct va. y with location, 

3. The information concerning the condensibles end transport properties 

ig not calculated or printed out. 


3,3 Overlay 3 (Grain Desig- and Ballistics) Outpyt 

The first output from this overlay consists of the current values of the 
variables in the SBAL NAMELIST. Values input directly via the $BAL NAMELIST 
over-ride internclly generated values, 


The next ser:es of output is a repeat of the informatior. contained in the 
formated input card:. 


The main outpu’ from the Motor Performance Module is of two distinct types: 
(see sample case) 

1. Results of the grain design calculations 

2. Results of the bailistics celculations 


Grain design results are printed att = 0, and then each ¢ime the number of 
bums becomes equal to NB. Thus, it is possible that more than one ballistic dis- 
play can appear between the grain design displays. The first column of the grain 
design output is the station number, The number of stations is NX + 1; station NX 
appears on a succeeding page, and station NX+1 is the plenum aft of the grain, 


or the nozzle entrance. For a submerged nozzle, the nozzle entrai.ce is artificially 
located at the aft end of the grain. The second column Is the axial location of the 
station tn inches. The next two columns are web positions between successive NB 
increments; the first is the beginning web, the second is the ending web, in inches, 
and the difference is equal to B10/NB10. Variations down the port stem from un- 
even burning. The fellowing two columns are port areas, in sq. inches, corres- 
ponding to these two web positions, The next column Is the average local grain 
perimete; {n inches. The next column is the local volume of propellant in the in- 
terval in cubic inches. The next column is ihe local bum area in sq. inches, 

and following that is a running sum of burn area down the port. The last column is 
a running sum of local void volumes, 


In the case of end-burning grains, where the buming follows the X-direction, 
results which seemingly vary with Y are artificial and should not be taken literally 
without interpretation. Port areas will be zero, opening to the motor area across 
the end-buming face; this interface will propagate upstream with time. Burn 
volume and local surface area also will appear .o propagate upstream, being zero 
both upstream and downstream of the interfacial region. (Web position is really 
longitudinal even though the output may make it appear to de radial). 


In the set of ballistics oriented output the first column is the station num- 
r, followed by a loop counter, This {s followed by local values of gas density 
(lb/in’), mass flux contribution (Ib,/in?~-sec), gas velocity (in/sec), pressure 
(lb/tn®), port arca (in*), weight flow contribution (lb/sec), burning rate (in/sec), 
cumulative web (in.), insulation area (in), burn area (in?) and the axial position 
of the station (in.). (Again, in the case of end-burner, values appearing upstream 
of the burning surface are artificia!). 


The suquence of output appesring after the columns are time (sec,), vacuum 
thrust (lb.), ambient thrust (Ib.), delivered specific impulse (sec.), motor stag- 
nation pressure (Ib/in*), nozzle exit pressure (1b/in®), nozzle exit velocity (in/sec.), 
insulation area exposed (in*), burn area (in”), throat area (in*), head-end pres- 
sure integral Ob-sec./in®; ignore for end-bumers), aft stagnation pressure integral 
(lb-sec/in*), insulation weight expended (1b.), propellant weight cxpended (lb.) 
und total impulse (Ib-sec.). Ail integrals are up to che particular tir-« in question, 
This summary output is Identified by the throat nede number, NX + 2. 

If this module is run in a stand-alone mode the above motor performance values 
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will be the only ones calculated, hence, nozzle efficiency, ambient pressure, 
etc. must be input. If a grain design and ballistics calculation {s combined with 
a nozzle performance calculation these thrust and impulse values should be ig- 
nored {n favor of those appearing in the THRUST-TIME HISTORY (see section 3.1). 


The unlabled output (6 columns) which may appear from time to time is 
debug print which was never reinoved (see subroutine MAIN3). The statement “Ref- 
erence To Uncorrelated Table Number" also is debug print and states that the de- 
signated table would have been used had the associatec option been exercised. It 


is not an indication that something is wrong. 


3.4 Overlay 4 (ODK) Output 

If IREQ = 1 in the SPROB NAMELIST (and ODK=)) the first output related to 
the kinetics solution is not from Overlay 4, but is a special equilibrium solution 
referred to as the "RESTRICTED EQUILIBRIUM OPTION." The format of the equill- 
brium output has been previously described in section 3.2. 


The equilibrium output is followed by the values cf temperature, pressure 
and velocity at the equilibrium (starting) contraction ratio (ECRAT). The species 
mole fractions at the Initial station are then listed, 


If IREQ = 0 the above output Is not obtained. The output described below 
is obtained in all cases, 


The first output from Overlay 4 is a listing of the input reaction and third 
body reaction rate ratio cards, as read. This Is followed by a table of all the 
species contained In the input reaction set. A reaction tavic, containing infor- 
mation from the input reaction cards is then printed, 


The next output is a table of all of the species that will be considered in 
the kinetic expansion folution. The numocr of the species, and its initial mole 
fraction is also given in the table. 


A formated table of reaction rate ratios (or the message - ALL REACTION 
RATE RATIOS INPUT AS 1,0) is then printed out. This is followed by a me-sage 
relating the type of nozzle geometry optinn which was selected and a table of the 
relevant wall geometry. The two values alpha snd alpha-bar are the particle to 
gas weight flow ratio and the mass fraction of the condensed phase, respectively. 


The next output describes the flow properties, nozzle geometry and chemical 
composition at the initial expansion ratio. This output is self-explanatory. The 
same quantities arc then printed out et a series of axial stations. The number ond 
location of these outputs being controlled by the print flag (J7PRNT). 
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When the temperature becomes equal to, or less than, the solidification 
temperature of the metal oxide being considered (if any) a message is printed out 
to indicate the beginning of solidification. Likewise, when solidification is com- 
plete, a message so indicating is printed. 


The note EP not reached means that the final requested expansion ratio was 
not reached before the solution terminated, In some cases this is due to an input 
error which causes abnormal termination. In many cases, however, the message is 
printed somewhat erroneously. In such cases the program actually has effectively 
reached the final expansion ratio but due to round-off error the message ts stil! 
printed, 


3.5 Overlay 5 (TD2P) Output 


The case title card is printed followed ty a list of quantities which were 
either internally transrsitted or input in the $TD2P NAMELIST, 


The following quantities arc calcylated in subroutine AGP, Only the math- 
ematical symbol anc units for each quantity is given here, See subroutine AGP for 
definitions, 


Variable Math, Symbol Units 
Gl) y - 
EqI) om - 
CPG 6 ft*/sec*/R 
9 ¥ = 
RCAP R {t’/sec?/R 
CPS Cs ft?/sec*/R 
HPS ne it?/sec*/R 
HPL ho it®/sec* 
: 2 
HPO ho a ft?/sec 
PR Pr “ 
RH@GO “90 lbm/fe 
5 B = 
RS B. ° 
GSB Ys - 
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The specific impulse, ISP(1DOL), and exit velocity, UE(1DOL), calculated 
from a one dimensional zero lag analysis are output at the nozzle exit expansion 
ratio. The mean particle size and particle size distribution (iny:ut, or internally 
generated) are then printed out, 


The next series of output comes from the one-dimensional inlei and transonic 
throat portion of the solution (see subroutine PARTIL). The output appears in the 
following sequence: 

1. The values k and m_ for gas-particle equilibrium. 

2. Initial and final conditions for the one-dimens:onal inlet integration 

(see subroutine @NED). 

3. Converged values for Kye Us Oe 8, yand f\. f.. f., £4. f. for each 

transon:c flow zone (see subroutine JAMES). 

4. The corrected estimate fork, currected values for Xye Uge a, @, and 

y in the transonic zone containing the initie) supersonic data line, cor- 
rected values form_andm.. 
g p 

5. Items 2,3, and 4 are iterated twice, 

6. The gas-particie flow properties Poe Sg: Uae Vg" r, 2, bye oe “p, 

and Vo along the initial supersonic data line. 

7. The nozzle weight flow, WDOT. 

8. The nozzle discharge coefficient, CD. 


A message indicating the selected wall geometry opticn is then printed, to- 
gether with the input (or internally transmitted) nozzle geometry. A finely spaced 
internally calculated, table of the wal! coordinates to be used in the solution is 
also output. 

The results of the ‘wo dimensional, two phase, supersonic solution are 
printed out along left running cherecteristics., A header is printed, for identifi- 
cation purposes, edove the output for each left running characteristic. The out- 
put may be identified as follows: 


Kow One: 
Item Header Meaning Units 
L.R.C. number LRC left running characteristic number none 


Ident. number ID type of point (see below) none 


emt 
ct 
st) 
a 


iteration no, 


VG 
THETA-G 
TG/TGO 


PG /PGO 


DG/DGO 


SDK/DG 


CF 
ISP 


17 


Rov s two througo hal 


Meaning 
r position coordinate 
z position coordinate 
Mach numb” 
gas temperature 
gas velocity (scalar) 
streamline angle 


ratio of ges temperature to chamber 
temperature 


ratio of gas pressure to chamter 
pressure 


ratio of gas density to chamber 
density 


ratio of tutai particle density to gas 
density 


thrust coefficient 
specific impulse 


number of iterations required 


ft/sec 
degrees 


none 


none 


none 


A&A 10ow is printed ter each particle size present at the point In question, 


item 
k, 
Re, 
V 
hy 
‘ 
@) 
" K 
T 
Py 
0. 9%. 
~h 4 


Header 
K 


TPK 
DPK DG 


Meaning 


narticle cize number 

particie Reynolds number 
particle: velocity (scalar) 
particie streamline angle 


particle temperature 


ratio of particle density to gas 
density 


none 


sa ak Na A ic ER al ll al 


Item Header Meaning Units 
0. /o DPK/DPO ratio of particle density to chamber — none 
Py Po particle density 
c p RPK particle radius ft. 
k 


The table below relates the printed identification number to the type of 


point calculated: 


Type of Point Calculated 


initial line point 

interior point 

axis point 

Kth particle boundary point 

wall point 

point inserted on the previous L.R.C. 
point inserted on the R.R.C. 


| ed 
~woaurswnre [5 


The quantity IP printed out at the end of a TD2P solution is the number of 
mesh points on the last left running characteristic. 


3.6 Overlay 6 (TBL) Output 


The case title is printed, followed by a Iist of the current values of the STBL 
NAMELIST variables, Values input in STBL over-ride internsa‘ly transmitted values, 


A table of the input (or Internally transmitted trom 77 2P) wali geometry and 
wall streamline mach numbe: is then printed out. 


At every pri:.cout station the program then prints six groups of quantities. 
(Note: DELF and DELI were added to the TBL output and do nct conform to their 
respective group headings. These two quantities ers not output until the wall 
slope becomes positive.) 


Contour Properties 

x Axia! Distance 

XL Are length of contour up to X 
Y Radius or height of contour 
DY/DX Slope of contour 
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Flow Properties 


M Mach number | 
TE Static temperature 1 
TW Wall temperature I. 
TAW Adiabatic wall temperature i 
DM/DX Mach number gradient | 
UE Velocity of fluid 

PE Static pressure 


Boundary Layer 


DELTA Velocity thickness 

DELTB Temperseture thickness 

DELT* Displacement thickness 

THETA Momentum thickness 

PHI Enorgy thickness 

H Shape factor 5*/6 

DEL F Thrust decrement due to boundary layer effects 
Heat Transfer 

HG Heat transfer coeffictert 

Qw Local rate of heat transfer to wall 

SUMQ Integrated heat transfer rate to point X 

FORCE Drag force in axial or X direction 

LAT.F Force normal to X direction for two-dimensional planar flov, 
DEL | I sp decrement due to bounc ary layer effects 


Internal Integrals 
ZETA, 11, 12, 13, 14, 18, 16, 17, UP, T2P, 3P 


Coefficients 

CF Skin friction coefficient 

CH Stanton number 

RETH Reynolds number based on mementum thickness 
REXL Reynolds number based on arc length 

REPH Reynolds number based on energy thickness 


Reynolds number based on displacement thickness 


In addition, if the sonic point start procedure is used, the initial values 
of momentum thickness, THETAI, and energy thickness, PHII, are printed out. 


After the last printout station the value of the throat radius corrected for 
displacement thickness, and a table of the normalized wall contour points cor- 
rected for displacement thickness are output. 


3,6.1 Dimensions of Variables 


There are six physical dimensions used in this program: 
M = Mass 
F # Force 
T = Temperature 
S = Time 
L = Length 
h = Heat 


Any consistent set of units may be used with this program. If input quan- 
tities are in units consistent with those below, the output quantities will have 
the units indicated pelow. 


RBAR= Lx F/MxT C, = Dimensionless 

C = H/M xT Ci = Dimensionless 

FJ =Lx F/H f) zL 

G #eMxl/FxS £ =]. 

TE #«T | = 1 

PE =F,'‘1° = 1 

rm weM’SxL 5* el 

F) zM/L o = HL? x 5 

UE «L/S SUMQ »* B/S axisymmetric flow 
Re = Dimensicniess SUMQ = H/S x L plener flow 
xX -#L hG HL xSxt 

yY -t FORCE = F(FAL for 2D flow) 
XL o=L LAT.F © F/L 


The built in values of FJ and G require use of the Foot, Pound,,. Second, 
BTU, Pounds, Degree Rankine, system of units unless FJ aid G are input. 
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For example, using the above conventions, and using the built in values 
of FJ and G 


Mo= Ib-m 
F = lb-f 
1 =eR 
S = gec 
LD = ft 
H = 3TU 
R =ft Ib, /b °R C, = dimensionless 
C, = BTU/Ib °R Cy = dimensionless 
FJ = ft Ib,/BTU 6 = ft, 
G = (Ib, /lb,) ft/sec* g = ft. 
TE =°R A = ft. 
PE = lb,/ft? " = ft, 
use: Ib /sec it 6* = ft, 
0 Ib it q = BTU/ft® sec 
UE = ft/se- Q = BTU/sec (exisymmetric flow) 
Re = dimensionless Q = BTU/sec ft (planar flow) 
X =it h, = BTUAt? sec OR 
YY “= 2% FORCE = Lb, 
Ae LAT.P = Lb, 


Note: When this program is used as an integrated module :n the SPP Program 
certain constraints are placed on the units by data commuricated from other modules. 


[ 
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4, RECOMMENDED PRACTICES AND GUIDELINES 


This scction has been included to provide a place for some of those comments 
that do not seem to belong anywhere else, but should be recorded. Thus, the follow- 
ing is a potpourri containing recom: endations on using the SPP program, miscel- 
laneous notes, and guidelines for selecting values for certain of the input quantities. 


The input oriented comments are presented in the order in which the vartables 
appear in Section 2. Amplifying comments, of the type contained herein, for the 
grain design and ballistics inputs have been integrated into the text of Sections 
2.7 and S. 


The best recommendation that can be made is to encourage prospective user's 
of this program to carefully study this manuel in order to learn the capabilities, 
intricacies, and limitations of the program. Time spent in this way will be repaid 
several fold when the user tries to run cases. 


Until the user has achieved a high confidence level in running the program, 
it is suggested that the modules be executed one or two at a time, The output can 
then be examined to see if it appears to be reasonable, and that the case that has 
been run corresponds to the one that the user desired. Th.s approach can save the 
user considerable computer time during his initial attempts to run the program. 


OM NAMELIST 


Many of the inputs from this Namelist are autom.:'cclly transferred to tne 
rel-vant calculation modules. These geometry inputs cen be over-ridden by direct 
input to a given module, but remember that such a geometry change is known only 
to that module. 


ASUP It is suggested that some extra values of ASUP be packed 
around the exit plane, average expansion ratio, and any 
regions of particular interest, since interpolation is per- 
formad in this table. Otherwise, e reasonable number of 
fough!ly equally spaced values will suffice, 


ASUB These need to be input only if subsonic information ts 
desired (unless an ODX solution is called for, then 
ASUB(1) must be input equa) to ECRAT). 
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TD2P requires that the throat be circular with a single radius 
of curvature. Different upstream and downstream radii of 
curvature are allowed in ODK. If the throat is not circular, 
RCURV must be obtained by a best fit to the specified throat 
geometry. Values of RCURV <1,5 should be avoided; the 
TD2P transonic analysis will fail for lower v- ‘ues. 


These values are over-ridden by the results of the ballistics 
module if BAL = 1 or 2 in the $PRO3 Namelist. 


The transonic method in TD2P will fail if a, > 45°, 


If the real nozzle shape cannot be input due to program re- 
strictions i. is better not to try to match THETA with its 
actual value. Instead, one should try to select THETA to 
minimize the differences between the real and input nozzle 
shapes (see comment on THJW under $TD2 Namelist). It 
should be understood, however, that if THETA is not equal 
to its actua! value the empirical 2-D efficiency calculated 
by the program wil] be erroneous. 


Remember that this table starts at 2, so the total number of 
points, NWS, is equal to the number of input points +1, 
Check the wal! derivatives tables printed out in ODK and 
TD2P for wiggies. Don't try to input too many points, or 
get spacing to close. If table is bad (not due to input 
ee it is usually better to throw out points rather than 
add them, 


If BAL#1 or 2, AVELS must be input In order to calculate an 
empirical two-phase efficiency. 


Setting ODE, BAL, etc. equal to 3 «llows the relevant loss 
efficiency tc be input in the ETAI(!) table. This allows the 
user to bypass both the calculated und empirical options 
of the program in favor of a value obtained from a different 
analysis (if the user has reason to believe the present 
analysis ts inadequate for his problem), or for any other 
reason, 


When this fiag is input as 1 the initiel data for the kinetics 
calculation is automatically obtained from a “restricted 
equilibrium” solution. Unies the user has a strong reason 
for wanting to input his own initial concentrations, etc. the 
use of IREO#0 is not recommended. Also, if IREQ=0, the 
srogram cannot caiculate unless it is “fooled” by user 
prepared punched cards containing tne Igy values that 


are normally punched following execution of a restricted 
equilitrium solution, 


F 
gf 
a4 
z 


SODE NAMELIST For @DE=1, ODKy#1 and TD2P#1, the MIX array can be 
used as described in Reference 1. 


PQ) For the present application the ODE program has been 
modified so that it can execute only rocket problems. 
Thus, only P(1) will be considered. 


OFSKED (I) In its current application the ODE program has been modi- 
DELH(I) fied so that only one zone can be considered, 


ODK USAGE INPUT - ‘The user is cautioned that ODK execution times become 
very large as pressure increases be ond about 600 psi. 
Fortunately, at high pressure, the kinetics loss tends to 
be quite small (a few tenths of a peicent of loss) and the 
empirical relation for the kinetics efficiency can be utilized 
without qualms, 
The ODK running times are not insignificant even at moderats pressures 
(P < 600 psi). The size of the present reaction set (Table 2-13) and the speec of 
the aluminum reactions are such that execution times on ‘he order of 2-5 minutes 
(CDC 6600) can be expected at moderate pressures. Thus, ODK should not be ex- 
ecuted indescriminately, Its use is warranted only when the utmos* accuracy 
(that the program is capable of) is desired, or when it is suspected that the kinetics 


loss may be relatively large. 


A completely satisfactory set of reactions and reaction rate data for alumin- 
ized solid rocket propellants does not exist at this point intime. The set given 
in Table 2-13 represents the results of the limited study that was possible within 
the constraints of the current effort. Should the user be interested in modifying 
this reaction set it should be pointed out that execution time is linearly propor- 
tlonal to the number of reactions and proportional to the number of species to a 
power between 2 4nd 3, 


With the previous comment in mind the user should use OMIT cards (section 
2.6.2.1) to eliminate as many of the unimportant (from a performance standpoint) 
trace species as possible in the ODE restricted equilibrium solution. 


RWITD The kinetics program allows for different upstream and 
RWIU downstream throat radi! of curvature (the TD2P program does 
not). If a combined solution including both TD2P and ODK 
modules is to be executed (and in reality RWTD#RWTU) it 
is probably not worth changing to a different nozzle shape 
for the ODK solution unless the downstream expansion is 
so rapid that the single radius (less extreme expansion) 
zpproximation would lead to serious underprediction of 
the kinetics loss. 
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PWRS and PV.ZS are already normalized if they have come from 
$GEOM. Thus, RZNORM is not relcvant unless ODK is run 
in a stand alone mode. 


The nominal assumed values should work well in most cases. 
User's should not change these values without first becoming 
quite familiar with the ODK program; they can have a very 
large impact on exectuion time. 


It is recommended that these values always be input as 0 
because of the speed of the aluminum reactions. 


It is recommended that the -2 option be used. The other 
options (except -1) can yield voluminous output. 


The user is warned not to change the assumed value. Tum- 
ing on this option uses alot of trees. 


To eliminate species for the ODK solution the use of OMIT 
cards is somewhat preferable to EPSEL. With OMIT cards 
the user has direct control over which species are to be 
eliminated, 


The nominal assumed values for these variables should 
suffice in almost all cases. The user should be quite 
familiur with the TD2P program (and be unable to run a 
case otherwise) before venturing to change these values. 


More than 20 points should not be :2quired, but one should 
probably not use less than 15. 


Should be greater than, or equal tc, 1.5 to avoid problems 
in transonic solution, 


THIW should be selected such that the Mach number at the 
wall initial line ooint is spout 1.1. THJD should be set 

to about 3/4 of THIW. THJW must be greater than THIW 

(the initial line must be on the circular portion of the throat). 
See comment on TEETA under SGEOM Namelist. 
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An estimate for this quantity may be made as follows: 


1 eat A ) where a attachment 


4 
om Ww 
ZAX ~ R. sin te + 1 a exit angle 
Ww 


This estimate is based or. a first order source flow approxi 
mation. 
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The nominal value should suffice. 


The results of the TBL program should not be too sensitive 
to this value (for further information see Reference 4). 


A value of 3 is recommended for IT@T. A larger value, 
while possibly a little more accurate, increases the risk 
of impingement. 


Since the particle size correlation is empirical the user 
is proceding at his own risk if thece values are changed. 


Should be left = 0. IPRINT = 1} 'eads to voluminous output. 


The solution stiould not be too sensitive to the value input. 
These thicknesses decrease as unit Reynolds number in- 
creases, 


The user should be familiar with the TBL program and/or 
boundary layer theory before trying to change the nominal 
values. 


S. SAMPLE CASES 
The Extended Delta Motor! !0) is representative of modern finocy] grain 
designs whose complexity warrants the use of a three dimensional grain design 
program. A sketch of this motor is shown in Figure 5-1, 


The input and portions of the output from two sample cases based on this 
motor are cresented in this section. The sample cases serve to facilitate program 
checkout, while serving as adjuncts to the input and output descriptions. 


The first sample case utilizes all of the SSP program modules except for 
ODK, Thus, ODE, Grain Design end Ballistics, TD2P, and TBL solutions were 
obtained, as were complete summaries and a thrust-time history. The second 
sample case represents an ODK (and restricted equilibrium) solution for the Ex- 
tended Delta, A complete summary performance page based on the combined first 
and second sample case was not available for inclusion in this manual, 


The sample cases are preceded by card image listings of the input cards 
used to run them. The following discussion oc the Grain Design and Ballistics 
input for the first sample case is pr2sented to further enhance the user's under- 
standing of this module's relatively difficult input requirements. 


The first group of cards present the grain geometry input. B10-11.92 denotes 
the web of 11.92 inches. NB10=20 denotes the number of burn increments; its 
significance is the number of new geqmetries called by the ballistics calculations 
(returns to HERCULES) without interpolation. SX00 denotes taking advantage of 
the symmetry of the motor, and 16 represents the fact that the diametric cross- 
section may be divided {nto 16 ple-segments. Such may be done with an 8-point 
star, each star segment being divisible in half, The next card defines the longi- 
tudina] grain mesh; the grain goes from X01#0 to X112*45.72 inches, and is herein 
divided into a mesh of NX11#100 increments. Thus, ballistics answers will be 
ees at 100 a stations, N $ Vv 
jon (0.4572 


inches inthis case) Js Dt (burn rate times input ! 
| 
| 
| 
| 
| 


calculational time cena The next card defines the erred grain mesh: the 
grain goes from a lateral distance (radius in this case) of YI1=0 to Y11=20.00 


inches, and is herein divided {nto a mesh of NY11=40 increments. Note that Yl) 
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is selected to be larger than the motor radius; this merely assures that no dimension 
will be inside the motor in case of error. Note also that the computer time will 

be proportional to the product NX] 1 times NY11. A100 x 60 matrix was the lar- 
gest found to be necessary to avoid saw-tooth pressure-time trace outputs in 
complex grains. 


The significance of the first group of cards is the definition of a geometry 
which is completely filled with propellant and divided into a calculational mesh. 
The following groups of cards define the burn surface by subtraction of volumes 
and codes. 


The first cylinder group (CYLI) subtracts out a cylinder which fills the 
central port of the motor, and is tangent to the star points as soon as they begin 
at about 1/4 of the grain length. The centers of the circles bounding this cylinder 
are located at (0, 0, 0) and (46.00, 0, 0): (10, Y10, Z10) and (X20, Y20, 220), 
respectively. The cylinder may extend beyond the boundary, 48 done here where- 
in the 46.00 exceeds the 45.72. The last card sets the radius ROO=1.88 in. 


The second cylinder group (CYL!) subtracts out the cylinder in which the 
subr 2rged nozzle throat region is snuggled. she coordinates of the detining 
circles ere centered at (34.22, 0, 0) and (46.00, 0, 0). It may overlap the pre- 
vious cylinder, as is the case here. A prior card inputs the rounded comer seen 
in the drawing, This is the BURN card with CROO and 0.75 in. as the radius of 
the comer round, The radtus of the cylinder is ROO=4,55 In. 


The next group labeled CONE subtracts out the conical volume that paral- 
lels the submerged nozzle exit cone. The surface of this cone, as the previous 
cylinder, is tangent tc. star points. The coordinates of the defining circles are 
centered at (45.72, 0, 0) and (39.40, 0, 0). The radii, respectively, are R10=7.30 
in, and R20=4.55 in. Note that a cylinder is really a special case of a cone. 


The next cylinder group (CYL!) subtracts out the aft volume and allows the 
end-fice to burn. The coordinates of the defining circles are centered at (45.72, 
0, 0) and (46.00, 0, 0). The radius of the cylinder is RO0=#12.00. (It might be 
noted that, if the aft end were inhibited, the second card would have to be a 
BURN card, with NB00 (no burn) as the proper analogy to CRO). 


(As 8 further note, tf the aft end were, instead,a division of segments of 


5-3 


PERU Onde eee pen Cree ne ee eee a ee ee BE oe te er oe pra S Seether ee Se ee 


a segmented motor, then in addition to this cylinder with a no-bum card,there 

must also be a division of mesh in the INIT cards. The first segment would go 
from X01 to X11 as here, with NX1l mesh. The first slot would go from X11 to X21 
with NX21 set to -2 to key HERCULES for a segmented grain, The next segment 
would go from X21 to X31 with NX31 mesh, etc. The reason for the no-bum cylinder 
fitting the slot is that the HERCULES scheme needed to be modified to properly 
compute the burn-back in the slot region. Other inputs required for a slot are 
associated with Table 2-12. The new schen.e was verified for the Lockheed Pro- 
pulsion Company 156-5 motor). 


The next two groups Subtract out the star volumes, The first of these, 
labeled CYLI, deals with the circular cut at the front end of the star cavity. The 
second, labeled PRIS (prism), deals with the remainder, 


The cylinder group, CYLI, attempts to lay a pancake cylinder inside a star 
so that the curved surface of the cylinder approximates the curvature of the grain. 
The corner round card takes into consideration the 0.25 radius of the star valley. 
Thus, the fitting cylinder does not have a sharp edge, but a rounded edge fitting 
into the star valley. Because the axis of the cylinder is toward the perpendicular 
to the paper, the X coordinates are constant. However, the Y and Z coordinates 
are offset from the plane of the paper and the centerline of the motor as appropriate 
to lay the cylinder inside the star as closely as possible, 


The prism group, PRIS, begins at the end of the curved surface and extends 
down the grain for a height H00=27.80 inches, Three coordinates 6"e required 
to define the prism, which lays inside the star, Again, the X-values are constant; 
the Yand Z values define the triangle at the constant X position, and the H00 ex- 
tends this shape al) the way down the grain where it overlaps the cylinders and 
cone at the aft end. A corner roind card ac~-unts for the curvature of the star 
valley. This completes the grain <iefinition except for the outer periphery. 


The next group of cards tabulates a series of pressures followed by cor- 
responding stra:d burning rates. Five entrics are used in this case; the maximum 
allowed ts 99, #9 this presents no limitations. The number 5 appearing at location 
6 on the first card below the logic cerd must equal the number of entries. The 
burn rate table is the first table at location 1. Location 3 specifies one rate for 
each pressure. It is good practice to input burning rates over a wide pressure 
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range to account for exponent variabilities, and to include a pressure higher than 
the expected maximum. The following group of cards is e corresponding group 
for pressure exponent, or slope, which is the second table at location 1. 


The next group of cards tabulates the axial posit‘ons of the propellant 
and the corresponding propellant radii. This is the third table at location 1. 
Again, 3S entries are allowed; this case used 36 entries, and this number of 
entires must appear at location 6 on the first card below the logic card. The 
inittal and final axial positions, zero and 45.72, must correspond to the zero and 
45.72 appearing in the INIT cards 10=0; X11=45.72). No other consistency is 
required. The fourth table at location 1 is insulation input. The input is values 
of time (first card) and corresponding insulation areas exposed (second card), ; 


Time increments are the fifth table at location 1. In selecting the cal- 
culational time increment, the product of this time increment and the number of 
burns (NB in the grain design input) should approximate the expected bum time, 
This increment will impact computer time, but is of secondary importance as com- 
pared to NX and NY. The next entry is the desired inte:vals for output. This 
latter time interval can be larger except where detailed study or plotting of results 
are desired; it has 2 trivial effect on computer time. If is is desired to have 
variable time increments in the course of the calculations, the tables may be ex- 
tended by inputting additional times and time increments associated with those 
times (¢.q., finer increments for Doost-sustain transiiion). 


The last group of wards are Table2-]2inputs, in the order of and with the 
units given therein. The 4, -1 ig a required key to end table loading. 
The first entry, TSTOP, should bs sclected to be somewhat greater than the ex- 
pected burn thc, If TSTOP is icichec before the computer time reaches DTIME, 
the computations wlll terminata, so there is the possibility that the tail-off wil! 
not be computes {f TSTOP is too short. The input value of nozzle efficiency in 


(15) should be 1,000 for an integrated made calculation, and can be something else 


for a stand-alone mode calculation. The projected nozzle srea in (19) is defined 
by the annulus formed by the outer nozzie material shell and the open throat. 
Valuc:: of ambient density (23) and thermal diffusivity (26) are, respectively, 
4.43a107 Wh /tn? and 3,0x107 “in?/sec. The user is cautioned to remember a}! 
-1 completion entries. 


Other tubies do not appear because they are not used, 
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SUMMARY OF RESULTS FROM SPO COMPUTER PROGRAM VERSION 38 OATED 3) JANTS 
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ALN 6 MEL a 3LAm 6 CL oe Az FeAFTV © sm 20,5 4 AB SKATING. ESTIMATE 
ALN « Ow & sie 0 9 Ae 1,OFTY 6 ne 90,8 , Ra S&S, 4768. FST vate 
fiw ¢ 8 B OL ef SF As &,9F91 » Ke 00,6 « BB AV «9 ESTIMATE 
ALw 0 Tw & Bl 6 4195 6 Az 1.0711 © Nw 90.5 6 Re U7ee «ESTIMATE 
Lv e M 8 aL oe Wa Ae EAE 1t 6 we 90,5 » FS F,760 ESTIMATE 
aLu e e. B® A oe <EL és ROFIL «© V8 20,8 « BE 3, 7be ESTIMATE 
Br © ACL PMALFL CAL Che BBS, EE} 1 0NB—0, 5.286. 15A4, estieare 
pac?’ whee 


retQS AXfiv 2Fae Garr Baveras 
ALE EAWSL 100 
wes? Caan 
ane 
Feed Selies 
en) 


JOAN  aede FOOSE 0) Foa, 


on - canabie Ate dents rate il niki hamster akan dal abaséc 


ee ee ae EF 


TreLE EVTENDEN HELwA VALIVATICNW Cage 
AS ne 
aatum 
ACialrppae, 1. vaSutale 
Ved LL Bae 
Eraarar, WN, 
BEUPLL a2 o heheh Ahe We 
NAS 1OBKRe 
TueF ah be 
TruoF FAD V0 
Ful ale dde 
DAijpayvare 
Beta r tale 
720798 WELILY S2ReR ict iat arth 
© £1 He OENTT RSME VAL MAIKO. TOKAI» 
Deep aL SPL PIPH oD AR VOVe! AMHEAR OD, WIFIVAT OZo TTPO SOS 
2 PROS. Do ANDVI 0. VSI od, WAIHT ~.5S1Z4H1 OF 6 
we. Salt. 
een) 
yoo ef t9a9 
erng) 
@aual t 14 Be 48 e OPK NIM g! eole 
THD FIP IOM, OL M, OW cHbr eT AA Gs 
Ut IST BKe 
end 
pera ee 
aga 
aawate Tw*Aele 
eR 


af eCranes 
mY, * 4, ce 169) =f 48 Pd | oveand, S 278,95 
au 1.* th a." % 299,15 


C %39% = 10,509 9 %, 002 Par } 1,794, $ 297,35 


ee ee ECE ON Cw NOR IY TIO DI LAMENT rE OR TNS 


MQvEL TSTS 
aang 
Feaates, a. 

Octt efale 

seraat. 

eatate 

% Baul ef « 

FarTat. 

weep el OM eOH Gali tp ears 
ee,3 

@F LETIAeg 
faa f 
f. 2 
ere 
oD 
a® 


REFS. Ha Foe FE SRRBIL 

ouF Fy Ve we, LAL ba) $2, SOAK epevhs 

neh, G@aFP}, aX0?, acyso.aF, er a4 344) 

har, SAF1G, teh Ao SOdSLF 8 er ee eee Oe SD 

Ba MOHD, © en Sq HET Ne re ee | 419 7G6 
. 


ens 
of a it 
oe ft 
a &@ 8? 
=u es 
431 
aA e@e * 


aeG 46816, wef, 


“a 6 7m 20%, 7st. eal h4. 44 Cee deere aueeeS) all -% at30?, 
as Aer OFT a, bai hh, SPITS, Pe lb ASTRA £900 R, 
Ye 7 id SL ik oe Bet, 4B te ote MR Be TU vt 
ee e@ a.*_t Ae. es Oe ee efie .*% ) Ww? 
@foe 9 east. e bat, th ate Fh #e *,% ‘ «1 
ate me Ae eo ee ee eS a wart ware | 
mo OW & ALAN ae «0, 7818 6 we 4e8 we She ECT eg tl : 
fae AL t © ALTPe e fe UTE le ee Se ee he 
oe ee) do U,9FIH , we Meh 9 BH 46% Egtt oaté 
qx% ts wfae 


oeoeovev# ees @ 


Ce tee ee ee ast,Pih « &ee,% 4 99,6. ASAE TEL 
ome ® 8 ead 5 ABt Putin. Beh, 9 BOD, S Ve Bein fe tyeeangs | 
en a Me O KAP ee Bah AEE 6 SOM, 6 t09, 98s agai fa Ut eearye 

C49 6 oP @ FA Oud 6 a2t,S6a « wer.h 6 FOS Le Timer he te98 48 


Ngee NII SO 


5-7? 


‘ 
ae 
faze it 2 29 © 57 6 Ot, WF « AwO, « EHO, Re MaucCee (POARILY 
7. @ 4VPR ELA © Smeg ORK, FEY wee,S » AOD, ESTy ave 
“9, @ GNPE a7 © WF 5 Bnh, FEI we-.§ . FFFLH, ESTI are 
Co 8 wo 8 ACE Oe lw BOT OF EY, Nah, «¢ BBAc%e e€STONRESR (AA 
Pr oe Fe 8 OL 8 tg MEAL ERE. rPee,4?, a8.) e wave § 11047) 
Po Me oe 8 LL 8 Fg ABET 4g eh 4 199.8 be FSTIMATE 
fife © @ Cu 2 32 » O86,F11 6 %ee,S « HED be FSTTVerte 
CL? © * #& MCL *& *: 9 BBD, AELO, Sa, ¢ B8Ve% eee eey 11047) 
Zur e HM © ELM © Fh + BOKER «NOLS 8g TEV EM EST WATE 
a9, © THB APR 6 7) Og ABE LAF IES Nae, 6 FBAA « Copeay 113067) 
> 0 % 8 me © Sg BOF, TAF 13, NEh,Ofe “8 BO, Rhe NIRS SHRGRG 67997, 
wo 0 Ne 8 wo 43% 9 BBP,39F19, 4B9,0%, 42 $180 NDeFh Mater 6FOGM, 
wae 89 6 HF O97 6 OBLFID 6 Met, 6 NOTA Ay Jain Ce CYQASd 4 
Nae Me © Jo to » BOR, FEFIP, wee 8%, 4B ALPRe Nde2} BQuLCw 619A%, 
427 © = 28 fee © M6 401, 60ETae wet. th, MFTA ARS NOo® MELLFES #1 OP)) 
a9 ¢ CO e 6. © FSP ehad ACEI Need eGo 3.715 6 at? 
eH 06 eo 4. o a? o82S,% Ertevee.d Aa 3,798 © 2en 
3X e CL os 4) rs OL obet CGF Tee. Be 471 © ne | 
® ° fa e &9 o obe®,% "11ers. Ge 2.8%? » 776 
eA e acy e $39 o Fy of PeMGF 1% ewyeeeS Ge Pe Fr? © ree ) 
BCL Pe « 9 6.7L 8 mt, ehet  ABctPeuechs§ 802,348 6 <3? 
a: © CLF = ALCL FL BAST Fizemere§ § 9Mat.a ° 236 
a o oL eo if © eh sh,d EVV onweieds Aad, O5% © 234 
aA eC€. s ae. e # ° e%,3 FY} esee%e§ oa? 962 ° 23? 
OA © LCL? © OCR © BLA Leet MSERD veh eh ofer, S50 0 20An 
4: ae Cd Cv? e a tie a Gs oben A Aneta ore} Bel ,& ry 74? 
Bye CNP @ alcPL eo * 060! AEFI} Ne rhe ds eBad.San 4 74% 
a tL ae fa SALVA 8 OM eB ae LIKE R os weheS «= 9Be27. 68 0 ray} ® 
afi eo €.2 @ Arey 2 6 oAal AEE guv]eh.§ oPel.4 ¢ 206 
aC Pe & @ alte, 8 el obeT MEK YL os oh.§ Pay oY Pea 
eA 6 why @ One, ow OS eT MEFS Per soteG ALY? 1k ey Fa? 
mC, e 18 Oe nl ee ePe lL FTP eneete§ oad, 88% 6 ae 
ar f_ 0 9? @ ney eo FF ohal  PGEt Lo 3-8. Ff AY Re 200 
OF Po BLOF & ALL 8 A ML Ved eI AEE Les eete § har, SHR 6 vee 
OP ve CL? eaectde * of rE MEF en wet. § o*et.& ° VI} 
6 Cy ,ae =» @ 8C1 2 © 7) ead AFI] os Be 1e§ o9aP, 80) » 344 
6% © Bw © 4 AP 0 wt 4 as Po.cKHY » Ae ALS , Me 8 ,496S TCT Emre 
arr 6 THe 47h o OL oe As PeAFUP og se - 0.8 . Be FR imh, CSTE ware 
Oc 0 MP 2 ete 6 Do Se PeAP7: «© wt 7.9 5 BO HEE? , OST I ware 
Bum © 297 © OL49 © Me Re LT .ABND 9 20 -0,8 » Be Bi arFd, FOV PmatE 
BLAP © oF © 4 6 I MPH be oFPy 6 Me 84S , Me GL mtD . KEV PwatE 
O,9P @ ome Fo 8 TDM de FeBFyt « te =tL8 5 BB & OFHG, FOP UQeE 
0497 © of, © rE 0 619d, fe oAhFZY 9 ~ 9 -7,8 , Bp F,arTeQ, FIP MATE 
61% © 93 6 Do 4198 S@ 3.0799 6 2 -9.5 , Be FOr ©» LPG wAIE 
Ay so 3H to AL ADH QU 1.02 5E © B= $4.8 , Se KH, aPar, FET earE 
Owe OLD w@ Me 01 NP 6 88 TA AFET 2 ae 6,5 . Me BLM ING, TOPLmATE 
6. « 285 © awe ww, Bo GAGE 9 2 eo 07, 4 Be FLHATSe FETT eg 
O14 © we © ALM ow Be L)9%HFs « Owe 9,8 , Be BPE . FLPo iy 
OL% 2 mF, o OLA ec CL oe Om FoF Tg te ede , Be GOR, F 
O6,* © Be 8 ate 6 + ee FARE ew te wes , he fete, F 
b.e eo @ eb, 0 Mh 2 Bo Be 8FHs «ee 29,8 , wm BL he og OSES 
6.s eo YW © aL wo WEY Of AS YTethE® «2 £ @=',% ge VL Pee 6ST Su. 
aoe % F€ Bo 4e Oe 6,471 29 ea 95,8, Be TL F%e6 tS°yu .€ 
ayo @ €: © a, wry ee QoPrst » ae of © , Be Q, Pune <P rwmaty 
6 © DT CEPOL CK CAL EL oe 80%, LEP et aod Ssuee. hae, Cc*ivare 
aes’ OF os 


“oe 'SD a%1e OF He BATH BEI 5¢E 


Orr 2950, 948 

ae? Gear 
aqe7e 
Fr eOy Faroe sPSuteoda 
@ean9 


FVAE 08 Foe, 


5-78 


ee ee a ee ee ee 


SCECEELECE C.sQmFar FeeCUTTAN THE Tete Cage © oTES MED SEREECAEEE 

QE AC WIP as 

fae 8 8B faite 2 BBB FEF IQ, meh, ¢« APPR, Ran Cw (pPORMILY 
of ? B CL & SL « O89,F1% Heh, « HRA, ON, SM NEN C1 AFG) 
rae r) oo PL fg AKL OED}. Pet, 0 FEIAD.ILP. JONAS THIKHAD 
wry s + 2 57 5 aeF,4fF 1G, wei, fh , BEO4,@ » N9e3A MSIL 61966) 
rr) a we © Hg ABILVCEILA. Kee, S4, ABQ? Afe 3,3 #Pheis 49197 u5 
seh e MIN g BET, 0AF IN, wal. Ah, As AL Bhs NS.5 PEE MN iTQAP, 
ao s Ao 1M . AuP, HRRELO, Hel, HA, ABEYR, Pe NOe) JOMNSTAN €1998.¥5 
Cre tL 2 ality? o Fat NFA we Meh Pe ¢,8 ® 39 
Or CL Perr 2 aLey3 o Aart. AFIiA awe AL Aa 8,8 r 32 
Os fo % @ are. o aN OEIG one 4 ofe A," e al 

AL ce ww bE Be UF LA . me ALS gf BB ALM e FST watE 

ao LMF @ OL APM , BE UL AELE 4g we 94h g MB Me ELT PMATE 

Od o AL 8 BLM . fe TARTAR . we 944 4 AE MLAs EST I wATE 

Or 2 = & ate ec as 3,96 2h we AS 5» ME Ae LSVPMaTE 

€vn 7OR UPas 
Ch oft% 8 C32 © CL + 809.F)) ~ne,$ , M09,%. ESTTaTE ; 
fn eo 7 & ro? oe 207, 98F 18, wef, . Med, 44, MAUL TC (PSCORILD 
Pye fw @ CNP 0 wt 5 BOK.AFET , Neh, 5 88 AA, Ray cee 1 OABILE 

COP © 42 BS EA ouPN « O00, EFA uer,G& , &048,. Tim FR CL OAPs 
Cade M 29 CD © 1D 5 Bet, OFT} , Her, 5 2850, 15s Bari Cm CVOMAMIEG 
o, @ TAPS CLA eo Tg AOE.FI} , Mee, 5 483. Le ESTI Mave 
Po @ e598 aft © 9F » anG, FE? , wer,§ , 87,4, Csvruare 
fy o@ 972 BS ate © 4 « HOT. MHI, wet, 5 Aea.® ef STOVAERR CAF 
Cr 0 Mm ow weL 8 Fg BBP, AERY, wee, kt, AB,t e waytD (QAP) 
PFi%o © @ WEL © S  ¢ ABEL FET 5 wel. 59, Kbe CSvivere 

(i, fie 9 @ TL © 37 « 68E,FH} , Men, S » MIRE, ESrieate 

CL> © 4 B Wh 8 Oo FOUL MELO, meh, 9 ABV Me CaugGav (yQary 

C9 e A SELF 8 Mm . OBE LFTS mee,% 5 BRL, FST hare 

mr_ 8 Tee MPM oe ML 4g BBE AF II wer,& . HEALA 8 CmB HA OT) 
wpe & eo Me 6 BoP, CAF, i, AeA Mh, OE 8A, ~%.1% B2,R8RG 41 OPHs 
ad a Bae BS tg MPT gg BHF, 19TR, NRALAM, 438 &,160 Sde70 MAIL Cm £FVVAR, 
vA © &Y SF YP © 99 5 20, F EY wef, « BSIA,A, “hor la itaegile 
te @ Me BS Fo MPA, FOE, PEERS, Beh, AN, AB A, THe fied’) @4 4.08 1t Jab, 
Ape ue 8 fee © 9 5 B83. 6bF lay FEAL AA, = MBYLS  AAe NeTA BE LFS FE BPliy 
af © C3 s 4, oP oper FEET LNA 5 oMe W7'S 2 air 
eo ef a a. ar J oder, t EV1.u20.% oMa 1,716 6 Pos 
be PN eo Fy a Ay o AML eset. 7 Feb eewe.s fe 3,878 . “34 
bs ere e @% obe%,% FUL erse.” eRe PLHP 4 Pt 
aM ¢ Ley 2 699 8 , eee} 26Ft 1 wee,S Pe 4,76 o 27% 
ad SL e {fi ot Se Sie eee eS ehs7, S86 6 rir 
a oie eof. 6 Oy shah, 0 FV ouseles ohoi a a exe 
& e oft 8 aft eo & 69h, % Et eweehat » 07, ag a aa 
eA ef, we &ly «0 o9eG, Bt Filemectes Mor, O68 4 PV 
eR o OLPLe w At © Pe FF gmat, Me UL eweet& ee nL ee | Pen 
re © @ ALP pd 0 FF ob ed POF Et angates ~*et le € Pas 
e& Cec e (7:39 ee Ce, e #4 obat, GF} .eea7% sh otB tar 6 Par 
Olas Aa i nt ed eS ee Oe. ee coe De POs? ead o Pos e 
a7 © CLP oo ALNPL 8 eRe lL PEF te s9e 78S + tte ® 46 
Pee er @ 8 PL 8 My e8a gAF Sev eahe& » >, 446 6 Fea 
6, e¢ wey @ 8 fh 2 eho d ASEtt ar goeheh Pe OhD Pat 
Pf, © fe e army ow re eee Cr S| ete? 848 Pay | 
Qe: fF, « AD e €. 7, 0 & od EePEF EL ewan. etat Fm. 746 
Oo Py Po OLOLF @ Flr © 8 PL Vebe .FEFUY ge Oe Mek ota hoeba ¢ Be 
Ory se (LP ett Ye ehot  A8F rt aegehed eet 6 e Lan 
Ce e rae © o why 8G AEE EP ar Geheh eMe?,41: ® yea 

O.% © We Oo ALPS 0 mg Be T,OFTE 9 we of,% . Be SF, KPO, FET PATE 

aye ie Ce @ aypP? o ME oe 2G teMFtT « nw of,8 , Re %,6 246, EST fware 

ee AFP wD aNP eo M , Pe 148838 6 ww 00.4% , Be 6,0GF , ECTE gre 

0% © B19 © PEMD a Men y Oe 8.8891 « wo 8,8 , Bo 3.0575, CETL war"e 


TITLE ExTENIFN DELTA VALTOATION Case 


GFA 


cS TAWECTrRYV 


5-80 ° 


TIME SEC) 


100,A6A 
PRAALES 


abil Ne rE ee CRE TT I Tey eT aS EP TT 


tTeajrcrnay TNFARWATION 


PoguA (PST a) 
ebonnoneenl 
2 IENHHAS OHI 
eVOnNAKEoeh} 
220N004 
2 tAndod 


9-8] 


48 


RTE ES ATTRA EE NR 


7% Fa 
RESCTANtS 
N 1, fag 
a 9,066 
¢c 7.3998 
NavELISTS 


" &, 0040 
af NhAn 
a «696.8000 «6A 


Cy feanna A 
of 4000 


20990 


@.%000 
of, GH6N 
=A, 8006 


=0 008A - 
of 0An4 
#8, 6A 


wl GONE VALUE RIVEN Cha SF, FraAr, Fee 02 FOCT 


706008 


2166080 
0340006 


SOCOIES MELWG CONSINF AEN Th Teg Sys tem 


392765 
t &s7A 
J 306? 
JT) sH8 
J 3764 
J #770 
J 5465 
J 4H 
Ji2s¢e8 
1 306) 
J127H8 
J 9708 
aged 
. Jd26t 
vw sant 
i &¢e3 
312765 
Jase 


4,153 
ari sess 
aL 
anew 
BL INR 
cc. 3 


J$9768 
J a47h 
J &ITA 
J A776 
J AAL 
J17768 
J Wnt 
J1 7/66 
J WA? 
J AHA 
J WAS 
Jp F7KE 
J13770 
L 13/66 
J1277a 
d19786 
J1 3766 
J 9466 


eL Ics 
ALCL DEL) 
ato 

aL Ir7_s 
at 3) 
ere 
rms 
e|ney 
eu 
ren 

fL 
fLan 
“tA 
woh Ss 
ne 
wht. 
a0 

"2 


J990465 aL 
J arte 
J 0764 
Jd 9765 
J wet Cc 

J37/67 Cw 
4769 «CN 
J hl 
J WAT 
J19/62 €3 
J 6766 
J as45 4 
ShPSTO 
32705 
399765 =H 
J @rhe 
J GFh4 
J #46} «403 


af a 2808049 
Enea, Fv 
(eHoV LS (EG Ceres 


Gq At VEG hE 
vt 
*L 
a 


al 
¢€ 


CRFECTIVE #LZt 


of 
i] 


> 


ME sOTa 6. PUREE 


Of Os 1, Fah iB 


PP PRPAAP CP wD Hb A 


+] 
o}3,338 
ot +, a%e 
ols, uF 


#/) 6, +&f 


eti,e} 
oe} 3, +66 
ot i kee 
eo! 5, mee 
ote, tan 
o>, 106 
ota, Was 
o}e, +06 
ele, hh 
ate,689 
ote, ay 
ote, 687 
oe, 878 


ag 
#8, 20A 
a4, 308 
oS, 436 


o%, +S 1) 


04, 19? 
0%, 198 
eo, 198 
ow, tev 
eo, 20! 
o%, 37m 
o¥%, 356 
oth, 391 
oth, ore 
offi, 214 
otf 318 
ot, eh4 
#10, e6 


¥e 


ve 


“pati 
eo 20 3 7ARKhAF OAT 
CLL eB ey 2) 
25987 99REF oN? 
o 125 21 AQEE on} 
oe RPS 7OSRRE HAD 
eo? VET 4OEPE CR) 
oS879Q9AG 1 CaN? 
eo LPZUP AAG OA} 
eu tf 
222, 77A of 9, O06 
Py eS a o1@, 26° 
oF, THR eA, T Tt 
0963,46 
oF2 FP? af8,7 46 
S101, 36 
22h, edt ath, 976 
oPA, SYR of HR, DB. 
@FA, ATA 21h, 9Ge 
oPA, Th wih, DAE 
aPY, 048 oPH, 458 
«FV, 39m FA  OPe 
oFt  G8h eA, 472 
oF2,99} 37,94? 
oFP, 7Eu oF, HO 
oAF FHF oA at? 
oF2, 192 oP), a5? 
OFF, OTE 094,459 
ePU Aah sPP APE 


a 
wht, tas 
ot? O78 


2! 3,489 


oi, 39 
#7 re 99 
ot %, 144 
OS mae | 
of 7, 846 
ot F, See 
ot 7,849 
oth, aot 
ott gt 
ot FF a4 
oP, 196 
ef%,*18 
we he eed | 


PFFECTIVE OXIHANT 


nO EPS 


Oe 


BOP CTs 


¢ 
olf ,%e? 
w39, 989 
o}3, 76} 


of?, 348 
oi, dhe 
oi?, 93 
717,008 
012, 6646 
oft, A%e 
of I, HFe 
ot 3, 8% 
ol, tte 
oft it 
ot ayy 
ot3,e78 
o13.149 


13,068 
a, 666 
olf, 798 


o3}7,80) 


#9P, 135 
ol, 379 
oi)? 196 
“17,130 
one) 
eof*, 184 
een, yee 
oF7, FOR 
of) 94? 
oF), dia 
of}, ¥Ge 
771. 9% 
#22 466 


§-82 


3,000 
2,668 


&,a8¢ 
2,008 
7, GAN 
1.666 
@.866 
9.008 
$404 
@,998 
1,008 
3,0ea 
2,090 
3,000 
3,668 


#©7n600,0% § 
6.66 4 
11706,68 §$ 


my RTURE 
uSuee 
+o 270970796%E 083 


aged) 
o S96 799586 of? 
e tha 38 G96F 2g) 
«S84 7995452 wa? 
e207 HAS Mosh 
PLLEL Lilie ty y 
elites} Tee cel 


208.158 
20A.250 
208,10 


Berets 
=6,000606 


Sean ws 


ot6, 838 
19,0346 
33.00? 
#38, 206 
o44,2716 
14,354 
@{5. 308 


10.604 
@96,SA5 
o9&,879 
o76, 78% 
714,709 
3h, 8% 
@18,%99 


223.402 
©23,,.%75 
23,816 
oF6, ahs 
@26, 926 
o26 693 
©2$,6?7 


oP? .90a 
022,460 
oP9,4 34% 
26, 4% 
* 24, S49 
«27, Aba 
Od Pink a) 


of 7,873 
of a,2a? 
oh, 266 
o1%.043 
of@,318 
o36,682 
19.05? 


#13215? 
o3P 77h 
of? AAT 
ot2.167 
o32,h} 
ol f 448 
o1 1 69 


@22. 7% 
oF, 3% 
=26.107 
o7n.570 
©2808, 756 
@31.392 
#31.365 


2.0006 
4,000 
2.666 
%.00¢ 
2.008 
&,006 
3,908 


ercrasen ttt 


DIME ELSE 20 AS PA AO et 46 TSA SPI RT NO ai 2 ; 
sas st tie 


ae eS cs Pn pre Seaeety oF 


eee th ee ee ee ee ee ee ee: ee ae ee a en oe ee 


THEAMETICAL QNEKET CERFHIUANCE AGCIMING FQOUILISATUY CO<PNSITION OUAIWG EROANSTION 
©0690 000006000050000606000500060608 060060860000000500000008006060008 000000000006 00556008050000000000H0080OS SEE 


GESTATCTEN ENUILTaRtUS DPT ION 


009008 0000858 06 0008 0690860006866 006600089000 HSCOE 06500000000F00 060660206 000560008 66600008060006000606000565500 


CF se 66,8 ecta ; 
wt eMactinn ENtMaLey STATE Temp = ss OES TTY 
Cetutca: FNawiLa «S€¢ NOTE) Cal sud, 0€6 « 6te 
Pet, eu ee ee Ye ee Pe 79608 3=— + 704 98, 008 $ 208,15 8, n0en 
Frey, &L 1.40000 214060 8.060 § 206,55 60,0008 
Fey C 7,249KH 86 1H, 80089 6 688088 14000 11768, 0600 $ 208,14 =%,0000 
seer, SFAFENT Fwy Se, FOULVOL ENCE Barz %s ,1ACSECG) STOIC “ieTURe Mavidsee, OF WIT¥sae. 
Cqau~agta = FRAP exty Evry Frtt Exty ext7 exit exty 
or ss LLL ey ee CS St ee, fe | 
#,. SETS 441," 221,23 $20,7 70,99 77,9) 16.3? 8.155 3.206 1,969 
T, ACA 9 4nS$i Shee aN2Y 4446 0347 aan Jang 3428 ame 
Me RTA Lh o79R,.6 == GHK,% = FOSLF «099579, He 1MP2ALS © 1960,0 9F150,6 9 2969,9 @9ERG, 8 
So MTHFR ad CID =>, UUMP = «2, NFAD = «9 AAD =—«s-_, 192MF —-_s-,92KF = - 7, 39MP —s_s«4>2IKD_—ss 2, 93387 ~=—s_s«_—, 89D 
NEw MIvsFTZz, 2296 6 t4S PY oo o FTF LAY ,169F HM g10GE0} .693f on? ~PPOERO? , 1 TRE COD 
. wa at 27,243 27,608 397 PRA 27.930 898,083 = 2ALNTK 8M. 1SO BBL POR OR MS 
th: vshieyt ST AMENGE SE LATTAS Of, ‘1HPS $e IHOTE SF NMIGA ©3.G6H79H =], 481KY 01 HON2?F of ON1GE 
™ vontee® 1, tf T.7faR = 9 EKG ULL, eee TL] Se Ce Se ey ke Pe bs 
CO BT 14 By CO) o@Tt? overs erst? «S687 6,090 SANG 26799 065A 24a 
Gbuseg 1¢) ay 2 ee eS eS 2 TS eee LG a Te Pa) Pa S| 
SM VEL ST She y88a,0 YY Pe ee) Sr, ee oe Sh ee 2 Le ee. LY re ie.) ee | 
mera syeRED A nana 1.0664 oh ee eee ee eee S 2.) 7o7ASH =, 5307 = 9, 0BDY 
VEi oF THE %, 3643, O7e, &IA2, vV99, rage, Hes’, tOA2, 0702, 
at at 1,6f0n 8 =6—69,904n) =0 2 fiane fi hane = 64,8002) «618.998 = Of6 88 8 Ge 
cirae, erste Sina @149 @10¢ &1ee e190 &16Q g1¢9 g148e 
ce ehad i7a ¥,708 4.626 1.519 1.474 1,73? 1,406 
TVS oh BO oSE Cry 9 30,2 397,0 23a? Paa.a 772,9 24s,* 09,4 307,48 
Te LatocsC rw ber,e da? 19,9 9aa,? 200.9 74>,5 aayt 200.8 
uM atluoys) BS, VPA 78 473 94, 309 245,240 76,948 26.086 24.18) 24,2728 26,270 
wo F FRaCh IVS 
At. oFFOEVQ ARAB BQ AMALIE ATNNAEPR CAH AREY ENBANZ = AER SFA CL ANHESA A, ONaROe 
Yel Ut eh ee ok ee ek ee CT ee ee Lob Mee Te 
o,7? eAA TERA =, AMOTTD = ,AATAVE AMER ZAHETED ,AAPENE J ANPEMe CON PANH ORK R RY 
Ye a) eMfAUGO FOF IAF AMAT OR 448HOT AFANTAPA CARPE AAT ETT LC ONASTR ANAM RY 
au~ eAAAF IH = FAADED = ARAAT A = AUAASAT AAR AN = AAA JAF E AAA OT RAAAS A BOROOR 
9 eAAALA® =, AOAMED ,ROATNE .GAAPER .BAARSZ «4 NfAERE AA AAAN B ORAAOA A, nARege 
ayer, o8THOCe = AMIKKD = AALEH® = ,CHTKAD ,AATUPE LAPT IAw HAP SOP LBARere 088188 
yf eAPAAFP =, AAAIAT ANMEPE EAA IPR ,FAALSE {FAAP SH 4 AAeATY CA AAARR FARRAR 
aye? oMAAP PE = AAHATR ,ABAAIY AQ DANAAP AAAAAA LFA AR AN AFP SAR 0 ARAAAN A ganage 
Yet eAPPEIA =, AAARAA {AAARPA = ,OAATAH AAAAZER 4 70N22R AAMT Te CORA HOL J AnAale 
by Bs ePAAAPT = AMARZY ,AHARPP OF AAD ,AARAAKL , FRARG JAAP ARY AAR RAR {ARERR 
a, 996, TeAAAHOH GF, AOADAA A, AQRELA A,GAKAHH HVPE} _HFTHA® ,API*AH 726M 4071609 
O, PAL) oARDEPE ==, AKIRA _AOPTR CFA PRE ,83RODO O,ONEHOD O,ONAAAD 2.089008 Hn, BBRneA 


Bl Se A tn NT A tre He ERR ee SEEN Mi MO he 


AREER eRAU RENE ARNE 


® 


oAhhQne 
0 364614) 
ohh hag 
of} 2Ga2 
089 968} 
afhAPAAns 
eAANHAlA 
of Meee 
Pe 4S | 
PR 
Ahanre 
ehOhF A} 
o FRY) 20 
08} 7689 
a AhANne 
oA@AhQae 
eo hherata 
onanete 
oA Ohern 
of FS WA 
mY Ty tt) 
o®AAP GA 
oO hHA TS 


@.48F o%8 
P.PSFonR 
1, JaFor? 


O, Woh? 
&, SAF AAP 
9, PAF oR 


ofannny 
« 7680 3¢ 
400004 
of 137ea 
07a 
eAQHORD 
ofhOaAia 
oG26%—R 


Sageesn 


1, *IFege 


a, 
T.1aFona 


9, 83 oar? 
9,1 a6? 
@, 


ehhaha? 
9 PRONG 
oh8°One 
et }tAya 
A by do 
ehOanns 
ehhAnys 
oh 34795 
oT PRRTE 
eAAAABA 
otAnnte 
eMAARKY 
o PRI4TA 
et h7A2 
oNGANae 
ongagna 
2 AAANA® 
oMOAQRE 
ohBansis 
aAPEGLT 
ehhneen 
ehtAAaa 
40RO?a 


2eh2Feta 
a, 
3.a6F one 


Oa AFa AR 
} ASE ath 
1. 75T oe? 


eonnnan 
o 249ga8 
2oranad 
eo DTASTE 
9 OHFA4? 
Onenan 
e8nana? 
2 8ARE3A 
2233108 
edGhnat 
s8nAnel 
adnraan 
o 39 %aRe 
ot Peag? 
Oenane 
2 8oaange 
Canary 
oGnann? 
o@agnte 
otv7etea 
o@anais 
e@earen 
9eena? 


eu | 
&,9}€en0 
fe 
a, 


ehhnane 
o FeR7EN 
AP 
ohh ange 
o HA9 206 
eannanas 
eOnran} 
9 606186 
PS YD 
afO nang 
ohOHANG 
@,cannad 
eo FARES 
el Aeat? 
GC, fHARAKD 
8anane 
eAARhana 
ehOranny 
Ahr naa 
e887 7639 
-annnn? 
o8HNAD76 
2 OGARAA 


e. 
9. FAP AR 


o9800n8 
eo FOR208 
.088an8 
eS I RTS3 
0 6O9H45 
@,cnnanh 
0860690 
e@h3eeC¢ 
0133383 
efORNNA 


0060006: 


@.600088 
2 30O9N? 
o 129996 

0.%e0nAn 
e6eeana 
2906880 
oe SOGAARY 
of90a05 
eO77ShRA 
MTY TS 
20083 a9 
2880800 


f. 


1, 06fo@A 
Oe 


5-85 


é , 
aAOnnen 
9 246605 
a OOnnnn 
of?} 48? 
o hGAWMs 

6, akOnAN 
oS8ANNG 
o8GL 92 
o B B9A646 
o8hPeaneo 
o89nnat 

6, 00NAKA 
9 186868 
Partasa | 

0, CaAAhA 

@,00nan0 
ALT | 
PSL TD | 
eOhKAh} 
sA77Ohe 
ehfnane 
ehOnne? 
eoananand 


098008 
9 26 29h6 
©, 6nante 
62585? 
ohnals 
0.980006 
040066 
CAO 262 
1278S) 
.0nanes 
-0annee 
§.0n0068@ 
371498 
o 189636 
@, s0nace 
0, 0080066 
0, caneee 
080000 
2000006 
ora 
0, 600906 
.08006S 
0, 000000 


e@ 
o 


@®@eeeese¢€e¢ee¢ea 


OO ASD D DOSS ES DHDBOSSODD® 


2880066 
o2e08i? 
0.090000 
oO PAnnze 
0000024 
@.¢e000n 
0.008068 
2000685 
03 FORD 
2600006 
&.0000fn 
0. 8ONOO6 
© 326726 
ol ARSE 
&,8esane 
6.008080 
@.409nGa 
0880660 
0.09n08@ 
oOPAZSA 
Ho OROGRA 
0800061 
0.600006 


wAnti 8. fa, fi, fi, 6. 6, 8. 6, ®e 
47? 2 HPFOAT APPT us 1 OIFon? 4, By 6. 8. ®. Be 
urn 2.45% oHK 1, AGECHA 4, 29F OHA 6, Oe 6. 6, 4, Oo 
Silat & SP FFA LFF OHA 9, AAFOGR fi, Se 6. a, C, Oe 
who 2. a AA é, q 3 GeFofa he Ge Os a, i Oe 
isk of M e. a, a, Oe t, 3 a, e, VU. Se 
Noe 6. &. @, Ae 6. Ce o. ®, Oe 
ya 4, 9We AD 1, PTF CGR 02, 1GF-OR Oe Be @. 6, a, Oe 
NINe e. a, Se Be Cs GO. @. , be 
my? 8. 6, ae fe Py Oe Ce ®, Os 
04 e. e, Je : Oe Ge Oe e. 0, Oe 


“Set FROCTIDN OF THs, CHWMENS BLES 

COMERS TELES ePPVOVTE «=| PPVGsy = DHVONG «=o F77ONY A PP7O11 «L277O2) «6 PTVIOTY «6, 27PORE «= LRP IORY 
m5: FC se wETGue W TOTAL. CONDE +S TOL ES tw THE “~auBERs 101.9612 

hee VALET I 49 wT IN THE CHAMRERS 8O,AN67 

FANSIAN HANAUETER AFVb0 o 13084 


WATE, EFTGHY FRACTION OF FICE, Ty TATA, FUFLS abn OF ONTDANT IN TOTAL Ont hans 


EDIE IAS 10 COmTMACT TAM crn rans ZONE ff 


TEwrRATURE OF RED 
PaTSGURE (esta 
veLnetity (et s%) 

SOFCIES MHC FOACTIWS 


acid 


Ct. 


.6n29Ger3E obs 
eSUOAGH TCE CHA 


4.976907 72E00% 


@. 


e 
© pAsTPAISEERAY 


BIRD ASIE SOD 
9 9096981 3E 08? 
@, 


e 
© 106) B7S2EOS 
® 109aagaiEots 


¥+ npeareeoe? 
9 1A07 SINE -03 
9 1S00AIPOL HOP 
9026 2QA81E HHS 
9 PAAER PORE HOO 
9 SPR0S 90ST OOS 


a 
* 9091 PAR IE SMO 
 3ani2701E 066 


-aarnanreeeoat 


* 1SAROrTAE SHS 
a2 TATE TE CHT 
9 arvarese oo? 
. 198an} e0F ets 
a, 


6 
<PER@AE VE CAA 
9 PA WATE ZAE CNA 


a, 
JV UngarprEoas 


Stapra rend ors 


oP 2?) PENA ROP 
a 


e 

ol eG neeoe obs 
a, 
0 34798) PIE WO) 
049 Sa2nkeE ong 
oF P4R7E6 FE cee 
oh 156 HO60E 00g 
9 199 7A1 OE oGe 
oF WMEKOTEE SRE 
of HP24R6OEE CHA 
o 3° 46004 7E oss 
PPT E TL Ry | 
a, 
a, 
ry a0 FRULGSE 0h 
2 1 ONT) 01 6E OR 
o 46 7RAQPIE WHE 
o 37 T81 FI2E OY 
o® 27 Gn TPKE OG 
oA O72977E SRE 
oP WOPnA STE OE 
miFATY ber 
o 34844 2OPC CAP 
0 IQA ARSE CRP 


6, 

2 7801) 48RGE og) 
r 

WETT TC TT Sey 


6 
@. 
246 h2a2 16F 09 
o TPHRVIEE Ge 
¢. 


steps iDd * 402% 2 ~e 


ee rin ec, 


33$93S983S (nin 


5 wet? 


- 4$¥9 Dirta 


Sell NULANI“NS LK 3dGCTD 69994949908 


$ Swe -e eesw 6 2¢ ree. » ce © agen &<- 2@ 
JaUmiaS2 *GZOV°E sy * H°Ge an * Lisuei ey * wy ¢ guile s dun ¢ vit 
24¥mlaS2 © doe®vy sy * G°Q- en * Lisutt = ay *vuecculg se aueuty 
2aVmLiS2D *S9en°E 2Q@ * SO" BK * Lisu’s sey * VW © 2ute 2 O19 ¢ uil¥ 
2a¥nlisd *S9e¥s BY * G*ge SN * Lige*t sy om © aevty os my ¢ uil¢ 
eit © UeS*Zeye Seucsvetiseerisve toe ¢ giply 2 vet oh ee | 
ite ® v* lane S* Gesr* 1 i seu® Lave Wee EIz%e @ 233 e123 19 
V0€E © Vow *Zag® = SPYCAM ETL ssU Levee TITy © = 61d5v 8 elute edtyte 
ove * Yod*esu® Seyosnsiianu® tage ve Wu @ cue Igig 6 
wee * thuetdsy® S*yome*tigve® i sy: se ¢ Toute @ wy e Ire 
av2 * thu¢esn® S*urBr*Lisvu isee we Yyute a Ign ¢ yly 
¥e2 & V¥n*esy® S*yoer*lisweu’ go % Ie © Wou"y a QD efty ty ¢ 
ee2 6 v° isu = Sfuranetbasu’l © «ly 6 lyulw 2) €W © lt 
bd 92 © Cuereday® SPuesrtisswal’ sve ton e Iyyte @ ny ef 43 i¢ 
eee ¢ vest 2Q* S* yore issue sye uJ © Iou7e 8 €03 « Ty 19 ¢ 
eve e v*isye e*Omosntt isos Ll sve Tg e 2@ta"v 8 eV ¢ Yoie 
vee ¢ YY¥y*euy? S*yese*hLisav*teveiguty «© late 8 eJatw ee 6 te 
d€2 © ton*esy® heyosselig wc °yszve ve Wty 2 J ¢ ut ¢ 
¥e2d ¢ ~¢90°dsyt S*yrsn iia v®°ysve ee ‘Tote 8 4 e rT) 
st2 bd 9° Lag? S°Prenelid &*yave 32 iyte s é13 « ‘9 
€€2 © 99g°%esH?e S*OSBVEEL su °dse® WG 2 ‘Igte «# mM ef ty te ¢ 
e22 © ove*t ay® G*-sn*Lignu* jaye tye vente @ Ag e ye 
Vé2 * @cv°2 suf qPoset tig u*yzrs nm © vie @ mi « ct 
tee * 1¢e9°€ sy? ie elie eo cme ee ae Ty a “J «© vec '# ‘ 
Uéed * ote°E sue B°osr Lbs Lessee cs e@ ‘ys vu ¢ ui 
412 * ote°e ay? SON *Lidsv isve guy e “8 vJe u'? 
{oz61: $37°92@ vleon 8US°vlah SUUcuEN “Vi sue *isy © Ce my awe ¢u j 
(e96t) mdIne@ t2°ON U4°O au 8UUT VEY °S'9Ga Sev © Ber © GC 8S my ¢ mm 
PIL6YOL) HO ING *O°ddew 8 fUEN Ff Ei sthse 8 2u 6 eR Oe UN 8 OW 
Cyyol) = M™INGVE O2°ON wis ey FUUTUBN *Eidoe cee © Yer 8 me os hy? ce U 
{02643 Sudvow 6f°ON SUM’n an “Ul oURY CCL aveO°C ae fore km z= Ge om 
(2960) Apoagnds © U*%ey © y®esy “Ti ae tse 8 7, oe Cee ce, 8 Tun 
a4Vnl is] *e°ben ° wPote § Uls'esP * jn 2 TS st yect> : Z 
; (249610) ApysHd *ur°eew * “yan “vidu*ese © We 13> ow we etd ”) 
24¥n13S3 *9e°Ean © w®mae °& ilatesy * 26 2 03 8 © oy lg baled 
340m 3453 *ee°Eay 8 BPmBd § bso uee 8 be Ve we om ey ly 
ee (4¥61)) Bw d2AVE bd bouy 49% ony *“Ulau°cse © | © Bon ws mye 13 
@9) Vesunci$ 36 *.°esu *§ “°yas *CR9e° bse © & 6 Toe we ew ee ly 
aL0n18$3 *vecew © SPoae * Lt steer 8 gue Ve aden © iy 
3i0nlas3 *w°faw © S°mBN § Listsee * we 6 UTS Ogun e fy 
UiWevel) ad Wee *oi°eSan * “ume © Cisectsey * dee ud 2 y eéug 
(24960) opus **Staw ° S°-aN * Yos*ese © Ucre U3 8 ene 603 
wn. ve . EUV00900) MIG ee Suv'leu ° °usm * Tldv’sse * mw 6 cud 2 hy © ud 
bTivvel) #3 Inve *En°cay ° “VEN “blac lay @ €cJ en oe | 
24V02 492 *w°Ceym ° G°9Be * The*tov * WD o €c3 2# vile og 
: aV4a"N oes und 
Ziv £483 *c°u at & wou Bh °F yiau’t sv ee "van ¢ uly 
BLUNT LSD PUP se 8 SoU 8H Visut2 sv: ful 2 Ute @ nO 
ah Id 2L9A14S2 *U°O ay * Boe BN gigucs fe Of mcute 8 éu'le © 
Si0aliS2 *0°v od * Gou Bh & O12u% ay * ply ane Ig 
te 6 e°o ay? S°Q@ BSNS Q9470°L ay ¢ uw s & ot3 le 
2¢ ¢ —©0°O Bye) = atg ame) play t ey * CW & ‘Weg ID ie 
ac ‘. ® v°o aye S*°o BN vigd*esy © 21, 8 WP ot2 19 
(y96l) NOAGNHOS 1°ON *2°wllew “uRri#ts Soljnu’cs¥ * ge v s eu 
tavour NKMIGd =6§°ON *u0°u su fUUlIBe “fal guo ase © vex a wy ¢ p 
(hee) “vnOue C°ON *Uv°doau fuscosy ei avy cay & we s cn 
(090i) HOWiey 9E°ON * ©°btan * uusn “el sueresy ° du eo @ Cube 
(9¥0l) Seuour CZ2U°euisw °° 880 fzsee’vey *§ WG ne Dare 
(14012 NGnO WSs -*939°9ee@ * “oar *RLte*ecse * Ine 13 ra?) 
(Tivevel) money °s*cen * 908M OM eeu ise © gtd 2® 0 ¢ 
Shu iD* ou 


‘ ‘ mt é a s, FERRIER eee SERPS Re rer errr bE SY SYR AE He AL AREY PPT PEPE ET? SUR TPL 


sareotGuten cnet Momeen IRE aE SARC ER 


ie SecusngE erpgnntnnecaryes nein tarot OM SRR TER EEN 


Jev> ase 
vel teres Ve 
Glave Ja0h S¥de 4th Leber) 
tv3e gse* 
2a0etasa SUNS 1 OMe*S ursNO Ll? “yes "TDive WlesZis we ‘Ff 
24001483 oL°E yf S*Oe ON *e Ilsves se ® Wee Te s W3 0 te 
24001453 o¥.°E 8 8 SNe ow fe fiasts vy omy © te 8 Le we 
9itwiliSs mde te F&F te whe fF FO oe See 8 Ss mm fF me 
BieniiS39 & ect sy 8 S*He Be Fe Thats sa Ge 8 OS a Oe eg 
JaveliSa “G9CR"S 4H © BP Qe A F Lear ay i a ee i J 
Bilas *Sves’s ce * wtue we © Tesues 2 ee es Ce ee 
BQivilas2 * tea°s 89 * SPD] BH * Fh gOrt BY ‘-w © ~. be ee © we 
BaUels,S2 *Sluvr2 ey © ogee are & Fada hr 3? a ee ee 7 
BavwliSZs *GHts*c eH © SIO AH * Liaues eae Re a a ee a 
Bliss d Psves*G tH ® yi ge Sw * Chara ty a ee i 3 
2iveabisa © wes’y =a © e®Ge By © Lasers ty ee a 
Bard idsd “Gueves zy % B®ve Be o Lidar SY fo weve Ft 8 Toe ¢ cv 
ZiWwlasa *evcre®s seo ® Bye su ¢ Tiaset 28 wees & GS ow & Cw 
Jivelisa © leees st * S° pw ae * [havi as eC we,te * » wo ge e@ cute 
ee ee ee ie ~ » @ + Gx - : @ eon + @- “ 


CPeCtesS tame 


{AATAINEN FROM OFACTIONS) 


i co 
a ar 
3 rh? 
& c.2 
s 
6 HCL 
i 7 re] 
] e ne 
j ] na 
4 109? 
4 . 13000 fae 
: 12 PO 
q 13 aur, 
’ we at 
: 1S ots 
3¢@ To 4 © 
wo 
$e aun 
ys 44 59 
ee aya 
zoos 
‘ at at 
pr eit 
20 low 
a) "> 
7 


5-91 


| 


2 
4 
= 
= 


Teron KADY 
Furap Aeny 
Teoen GAny 
Teen PANY 
Furead anny 
tute) 6°O" 
yures AaLr 
wureed we 
ween atv 
acer aaf¥ 
feee® ARDY 
qeeah Sane 
toras 80 


e085 ante 


eeact ts ’ ae 
AE ACT ONTS 
earspur ts 


asacttiw ? ae 
SEACT ANTS 


PRNTYUE TS 


acactiw 2 ae 
a¢ actauts 
eennur rs 


agactie: 8 66 ae 
a¢actayts 


pauls 


agactinn %& fs 
a¢actauts 
SRRHuUL TS 


agartiw 4 se 
wg af PANTS 


sunnyr t% 


egar lt 9 ae 
efectavts 
Sanne t¢ 


ecactem 8 ff ee 
or actar’s 
eugnve rs 
ag ag t?s4 ® ke 
eg ast aut 
Beryl rs 
agactrwm 98 oe 
e¢ 6c T4yTs 
oa wl t¢ 


avacte%, 3) 48 
were tarts 
sesrue tS 


avartitu UP ae 

at ac tevte 

pex.cur r< 
aga¢tt*s ty se 
ot eT ar%e 
a tipi FE 
aeaett% te a0 
ecec touts 
ante *¢ 


erar ety 14 ae 
Oc ac te.*s 
aeteur t€ 

qa 


eg act 3% ae 


.Pactants 


wgectran reac€ 


0% :SOnOk os 
¢i 
tC? 


o Pf AAAGHE 016 


ce 
eu 


eA PANAORE 0392 
wrL 
” 


ePOQPARAE OFA 
wre 
” 


oP ONAROE 1D 
bet 
od 


9 WAANRGNE 030 
) 
HP 


ePAEPAHOF Pt 
ap 
@ 


o TOHAHOFE 01? 
ae, 
aLPL? 


9 BARAAA TE OVP 
iL? 
ou°Ld 


by ald 2 


Tes 5 
e ary 


eV BRAM Fat? 
a 
rT ad 


9 WeraaHaA ot? 


& 
TS? ae) 


g BANANA el ® 
fie 
ay, Adee 


9 WAH MRPAE ON? 
ao 
The 


a LAPRAAHE 00D 
a 
ene 


ot FER ERE oft 
ea 


Os 
a We 
te 


ee 
ba] 


oPeneannt edt 


o2beSQnnd oO? 


oleai ranged? 


eoSennaresd? 


»PPA000NE 097 


@ 


ehgaraangeed 


A, 


9. 


9 WMRESA AM OMY 


a PRWHBARAE OAS 


we 


he 


2 2A0ONNCE 005 


we @,Seeenestos 


oh hesaeer® 009 


9 1 0GORORE OO) 


Seqenang ets 


-deeennng +00 


“e 


o*egeanss ota 


9 S0C0RR0E o6e 


SO HSC RAs 0OQ 


he 


9 5408008 o6@ 


98aneenns off 


eo 0, Sagqanrtoea 


we @, 


en wur. 1¢ 


@€arti3w 17 fhe 
OF ACTANTS 
PRMNUFTS 


eCarlisn a as 
Stactartg 


Ba Wue»rs 


@ da 
ogacTauts 
eanQuers< 


Qari we as 
SF acTanuts 
OBIE TS 


agactinn 93 te 
ecaclavtg 
sernur rs 


2¢artiss 27 te 
Og aclants 
Samu TS 


Ogartiw 23 Ae 
wr ac Pauls 
BEM U9 TE 


Otarlim Re fe 
ar ar layts 
O6 Qhur rE 


Ofarti¢w ee 
OF al Tayts 
oosnue rt 


ecerfts 36 as 
Or ac rau*s 
Ber me 76 


etartise FF ae 
Par fa 1%@ 
Ba Wear 7% 


or act iw 


aeeetiwm vA fe 
OF SC¥ OutE 
pewmer re 


oeartow Fe tes 
egacta> 14 
SB*ALIF FE 


orecta-ve 
DE Wyer Fe 


orertiw 3 ae 
SF af Tantse 


Sawjur es 


ww ae 
er ertagute 
de Wer re 


arartt ww 


cae 


e®4PRRONE O12 
¢n 
gare 


oB48GRONF 0nd 
Ci? 
ca 


ot FQAGGAF fe 
cn 
en F 


oe SRARQERE ot P 
ey 


ethngatafc,? 
Cu 
art, 


ot PuOROns ore 


1 
wr. 


o P*Qaager ot? 


a HANREAF Ett 
at 
om 


eS *AQngne ey? 
e =P 
tt hal 


a RHE REO OIF 


ooPt, 
mer 


PY te eo a 
rt) 
ca 


PY sh bees TE | 
ap 


TeneennE oo) 
eS SBQRODE oN? 


05818 enF 002 


o SSA AN ANE CO} 


252 3°3 g62 422 


ePRAAAAAE 08} 


oA 2O0RANE CC) 


2 thRgannl cee 


o Seonans es) 


o Wensnnger) 


« 9BARGAE of 5 


2 V7Rb ARE CR 


ghBI9F HAF 2G) 


owOe: 


“30R ah) 


Sh Lee. | 


Pe LS Lor ey } 


o PRaaL A 8E a%g 


Ge 885 See cos ARE Prk BES Bre > 


Ne 


Ne 


Ne 


Ne 


Ne 


“eo 


@,S5A300 HNC COR 


a. 


=oSh000NCE O89 


@eSONGANAF COR 


Ge 


we*?7IhE DFE COR 


©, 30HOA ANE OCS 


+, TeRSASAL eA, 


© 64 SARA ae Pa 


«. 51 G88. 0e4 


3e 


&. 


i 
q 
q 

2 
# 


8 

dPACi TIN Cg as 
REACTANTS 
ORAMUCTS 


QFarltw 36 As 
SE ACTANTS 
Pannut Ts 


araetyan 4s As 
HEACTAUTS 
Pe Mere 


QEACTIAN 36 Ag 
BFACTANTS 
PONnucTs 


REACTION 27 As 
REACTANTS 
Pe»ADUC TS 


QasacTIoN Yo As 
PF aCTANTS 
pasnucrs 


araeTrys = §639 as 
OF ACTANTS 
PUSHUCTS 


aFacTI®WN 00 8 as 
ee actanTs 
PRINUCTS 


aeEactTiaw 4} Ae 


. eae Tarts 


PLINVUETS 


agarTiayn 8? As 
SE ACTANTS 
ELAINE TS 


REACTTAN 69 as 
er actTauts 
ewsniie %4% 


areeti3n 6 As 
weactauts§ 
eanmyCcrs 


asactTtan 66 aa 
BFACTAYTE 
paaALc 7s 


ararTysn | 6.6 ae 
Meactlat® 
Benue 7& 


aFaeftas 67 As 
CER ACTANTS 
eesnye t6 


aract1an Of Ae 
gear tavts 
BQrnierys 


aFartian es Ag 
TAPP Ew ts 


a TeOHAONE O15 
0? 
t.t4 


a 10SAH0NE ot? 
a4 
aL 


ayo 
aL 


2 VOSHRONE +19 
auan 
at 


eR INAANGE CT? 
AL 
iO 


cA NAONE OV? 
xe) 
as rn 


a 11SHHOKFE oO? 
ALCL? 
r Yee & 


eS INAAONE 019 
aL 
acc. 


eSBNNNNAF O12 
aL 
aves 


oe BIAHAOKE +12 
4La 
oye. 


eT GAN M RF oO? 
a0 

le 

eV AGONISE LP 
acer 

ALeL? 


a VORANNEE SY? 


AFL 
ay OF. 


at GARAAF+Y1 
&L 2 
Bal st B 


yASAAerF el? 
Ally 
aye 


a ORNALAR CT? 
ari? 
ay WP 


a1 TRAAMF Oy? 
a4 


Aas ,PHANOAAR OO? 
4 
fn 


Re 29150006 ¢61 


ta 
ths 
Qa B7YSKANESOY 
0 * 
i? 


Aa ,3h71NOKESO} 
cL 
aLcsee 


Qe PRT ANNE OG) 
nn 
o¢ 


Se, tch900NES0} 
aLci 


* ¢kl 


Re , ASHKnnnked} 


az LEMON OAE SO} 


ae ,28S300HEOO1 


ae PAS INE KE OOY 


at ,ThRATAAE OA} 


TSA AAat ed} 
ant J 


nie 


Qe SPF HARE eA} 
ene 


s (Sh 


te |. FRM ANG eG! 


An 


oe ef 


Az LTA TAINAN CG) 


esos 
ue - 


wh 
Be. PSAANAAF ee HY 
“ad 


AG a EH UEAME CH] 


é 


Ne 


Ns 


Nes 


tee 


Ne 


Ne 


Ne 


Ne 


Ns 


Ne 


tie 


Ne 


Ne 


We 


§-94 


~.5006NNNE S04 


=. 50N0ANNEs ON 


= S0nonnneeta 


= ,S000NHNhE SCO 


= .500000NE+ 00 


o  S0CONHHE eA 


@~,5SHOCShO0Esbn 


~  S000NAAE SCS 


@  SAONDANE SON 


= eSONOHANE SNA 


< S509AN Eee 


@, SAK 0NNNECNA 


= 5HCONANE CNG 


=e SHONOHAE CHD 


e.SHHgnnneeta 


oe Snnnnangete 


sl 


ee es na FTE 


‘Parnucrs 


MEACTIAN Sn jae 
REACTANTS 
PRNDIUCTS 


aracti9N 6): Ae 
AFACYANTS 
PeaNuE r< 


Qgacliow Se as 
OFACTANTS 
ppnnuc ts 


REACT TIN S 3 As 
ME aCTaNTs 
PROnUCYS 


aracl isn a se 


DESC TANTS 
Panouc rs 


acactwww $¢ As 
CF actants 
ePRnnuCcTs 


REACTION 6A = Ae 
CFACTANTS 
Pa *WUC TS 


&7 és 
OF ACTAYTS 
PEANUCTS 


afactiw SA as 
OFAC TENTS 
earners 


aeaettiw 


CEartiw Se as 
erectaurs 


PEAnYyr TS 


Sractian 44 os 
REACTANTS 
Sa IHue 7¢@ 


eparctiToN 4) As 
Braclauts 
PAVIUT IS 


arariqsay 4, fe 
OEACT ANTS 
PL NUTS 


aractiw ay fs 
OF ALT Ayre 
eoqnie F< 
eFar¥ts, 46 Ae 
ar arTayts 
Pasnuc rs 


efarlisy 4» Ae 
Se actautse 
SusniE TE 


ee 


acete 


acecy 
Tmuga 


oVASCHONF 0}? 
ouee 
cOCe 


eo VOSHAONE 032 
acct? 
am @ & 


ol ASAANHAF 012 
aLCu? 
Aacce3 


oT 8S0R08E 01? 


cL 
aLCL? 


eo LONFHOHE 1? 
aye 
4.02 


eVINANOAE CT? 
aca 
AO? 


eo VAONQAFE 04? 
a, 
Ay Ap 


eVONAANANE 1D 
a4 
aro? 


eo V9ONNAAF et > 
ayte 
‘4 


ThAAAArE ey? 
ard? 
ry 


oe ONAANAT et > 
a nie 
e% 


ot OAAAREF ot? 
3( On 
* 


e*f*ARACHF et D 
Ay fom 
& 


eo VAC AMCEe 2? 
6, fu 
tr 


eo SORRAAL E05? 
a, 
ayru 


oe 

As .2AS3000F 60} 
e fx 
o 4 

8S ,3P9ANHNE CO} 
o "2 
te | 

9e =. 2S6HNNNE OO) 
® ALL? 
¢ ALCLy 

Re ,1ADOHONE oO] 
¢ tui? 
an = % 

Re ,2S9) arated) 
® «4 
¢ aL 

Re ,SAAARNAE CO} 
e (i 
om 

Re, IS BARNAE CC) 
¢ ¢(_4 
* tL 

Qe AR 20nAE CH] 
@ Ne 
en 

Re ,BMIAGAAT oH} 
@ 8D 
oe a 

Qe EUAN ASN 
« 4? 
+ Al Ved 

ae SF werk eD) 
e Ate 
e aL" ae 

Is VSN PNR MY 
o ae 
@ at 5269 

Re {AS Rane Are 0} 
* 
@ bye tPRe 

CT RS Oe Yi 
» At 
¢ &, A908 

Re PS 7 waFed} 
nt Wa) 
© AL APH 

Qe LS FARPGAE Oy 
@ APA 
o ww 


Ya 


he 


Ne 


We 


Ws 


Ne 


his 


ve 


§-95 


-.550000NE e606 


-.S00ONNOF SG 


=~. 50NOnNAECOA 


SOKO NNESOH 


~.5tO0NHhEeha 


~oSAN0NNHE COR 


=o 00ONNNE S66 


~.506 19AneGA 


=. S000NNOF ofA 


~  SCOdNNHE COR 


=“, SONF**906 ot 


we Sir iAnseta 


oe SANMAANE ofA 


e,SOA ADAP ON) 


we StAIAANE O14 


=, SARAAAAF COE 


weaceton fA Ae 
OF ACTANTS 
earnue ls 


agactiww 47 As 
Pg acrant’s 
PR MUL TS 


agartiqa «A as 
o¢actants 
sannucTs 


ae as 
AF AC TANTS 
eannuc rs 


agactiw 


tt ée 
AF ae TaNTs 
Panoucts 


aganly ww 


arscTiw ?} Ae 
agactayts 
eaanucrs 


agacticn 72 fe 
oracteuts 
pennucts 


aracttse 72 as 
ar ecTants 
eanpoucTs 


e TOHONOAE C1? 
ao 
Ann 


.1Conacere)? 
ayo 
aye 


2 1ORDAQOE CI? 
a. 
ay mm. 


eRhOWNONE CTA 
at 
aL 


a 1OGNAOHECT? 
Ath 
AL 


~SOQRRONE O12 
6. 

aL 
SVICRANEOL2 
Ayo 

aL 


eSNHAHONE C12 


at 
ALCe 


oe ¢ 


Rs 
0"? 
| 


a731000E 60) 


Qa 
CL 


eSATOSHNESAL 
SA2ERNNECOL 
Tene RSAEar 
2 37A000NEC Ot 


ga 
t 
nv 


2 374000NES 0) 


as 


cL 
“CL 


o37SCOONEOM) 


au, 415P0NNEON 
aLeLe 
@LPL 


Ne 


Ne 


Ns 


é a 
©  SONGENOE 060 


@  Sh00080F 006 


~  5000NARE Cbg 


@, 500000 S86 


w  SONOGHHESRE 


@,$085 RARE SOSH 


~~, 50n0nnneede 


o. S0000R0E 906 


5-96 


Sececren SPECIES FOR KTNETIC FRPANSTON 


AL her 

ay fe 

a. 0? 

OL O90 
aL on 

OL 22258) 
ar anmatyd 
ef 


0136703954603 


+ 482031 95°E 02 


9 VAASPP LIE HH2 
oV6G197S2F 63 
919260081 E=04 
91 NA 7IS2RESNY 
23540 3270E 202 
2320 19953E 263 
9 PAE ONZORE 08 
»82K05905E-03 
PhO) 7653F 008 
Se 
ghOP NGO 7GE att 
-BEAIWNG IVE OHH 
ePUOTAATIE MAL 
Verte Vrreogy 
eS IATHISSE HOS 
ePSOSRGSNE Oa 
mYSLLIRS, Cus! 


2) PEBPSOTE ore 


sUNAenlarledn 


ePATETOISE 669 


eT PTA VL SEE 000 


ear nnnreertd 


PLY S10) dole) id | 


.£bhan? 3HF eG} 
.o 1WPTAT LE OG? 


TRIP VIELE ote 


§-97 


TVERT SPECIES 
reat SPECIES 


3 
3 
| 
3 
i 
a 


q 


~~ 2 SC BPuPRCOwEN~ 


ae oe 


a, tate 


OuwZ5 et 


o TALOAE CHA 
oVPalSeon) 
eh SaAneeny 
o FRROPE OC) 
4g MANE ONY 
stSaGreen} 
eSSPUNT OAS 
eAAGSAF 06} 
Ahearn} 
eT AK26E 0%? 
ol PVOSECH? 


o 98693 


NISSACTATI AN RECHERINATION GEACTION RATE RATIAS 


Owe tt) 


a PS9NF CR) 
UAT I 2E009 
2 1ONRGAE CARY 
2 ARREE ONY 
o F321 POOH) 
oPI374F On 
o WPAOGF OFT 
0 ABD 9Ce 1 
oS BOSte: @y 
9 AOAPGS 0} 
oSSS32F 01 


6. 04s Bae 


alt @eacTinn mate MATTOS INPUT AS 1.6 


CO TAVQEN WOZZLE WTTIW SHLECTEMs wate TeA.€ 1S 
- WAL: SLE 


PYS 3A a7) | 
PY ali 
eh PROVE CNG 
94 POOOF ong 
eA hhG2E One® 
948) FOE one 
26h263E ong 
2 19) SOE Ona 
e WPAGIE one 
9 PIRGAE ORG 
SITTER £11) 


yyy XT YY ic tiisy ici yi rye ier iii 00690000060680 


Py thd D> | 


SLESTCESSSE CeMMENT EXFCUTIAN THuf Tats Cage & 


> 2 = Gf Gu 


FLfa Pnreratres 


~- + ee - 


tNTTiaL COWwTTTIONS 


KTHETIC EXPANSION 


WOZZLE 6, NME TRY 


2 


eaCe SUCRE e o PSASSGR0E 00 THROAT RANTua CFT) 28 7861667C 008 
VELOCITY CET EFF) 099434 979F 00% THROAT WALL @aNlUC OnSTMFae ,20R6aSEEE oO} 
SOCSK INE (OSTA) eT ISHSHPAE COD CONE ANGLE (N€6) 9 2300n0eed e062 
WaSilvy (LmsFT3) oPPTITATIVE COO FXPANSION @aTtO 8. 
TENDFEATYURE (NE AHR) eARPPSOP UW 0 O64 CONTRACTION natin o2 SAG 00C0E 85 
Fulscaley (QT107,8) @_ TK) LARA COS INLET ONALE (NEG) 0 34eseoeet +e? 
GAS WOLECIG 8 UF TENT o VOASRG1 IF 062 IWLET wary @anrus O- 
: mgay Cavacety (PVus_AehE Gen o4APISSINE CHR TeRoat #0, @antus veSsTAcau _2eeCeCCCE oO) 
3 dates 2 11075290Ge6) 
3 WHET O8 WELINT ICTR} oPSNOIATNE CER 
z euEMica, COMPOST TTON 
wo, SPPeCIES MASS FRACTION une FRACTION wn, SPECIES mass Feacting “OLE FROCTION 
| at PERS ek | ol 347TI ST PE HAY 2 ALCL 0129014600E og) 232238841 C 0? 
? CL? ol S1ARA1 OF 20) =, 98 28 8 09S K WNP & rye 9 & oA .672OF ong =. , BOK A FERECOS 
§ ays oP TWAGEPF OSG =—-_—_, NO NAENO LE OHS 6 ALO ol P2OR2EEons =. . 100 F020 2E 09 
? aunty eSOOSI ANGE oLP® = =—s_ 3% 14.1394 5E OA? a aL Om oTVSPO} 40 HNq =—«s_-. 6 800 9K LOD 
e a, h2 eS POREAGMFehd = =—«-_- . FF LOKI ST HOO if AL new elBOveZteeeep §«=—. , 820095 0583 
v1 6.70 PScU RAC SS 20! MEE et 1? €. 2031S} fe Ce 
13 AuPAney oPTVOMAMONE «=, HOPTEIOPEoe = 18 CO o27PTNOSLoNH =, BOnSRAEL SES | 
i 
YP st 2 ePPSAT TORR CH) =, DTANTILORCO) = 1H CL olederezetog, =. 16P72ESIE°08 
; 
7 «6&8 oF REEAIETH CAR §=—«- © 18 7902 97 08 a Ck? eAD17ICTOH Re =. 15060 004E 964 | 
48 ” et BAT WahWoePP =, 28.79 994GE OM 7A weL olARTALP2ZE Cag =, 1 ESSALASE ES 
2) “ap eo APR PAOAW RGA =, BOOT ZRF OK >? ne o728707TESE on} =. BBTV ONSISE COO 
2% ~on oAI40AG0CF WA} oT 1 7O19SPOF OHH Je nO oehFP2Q7OTEohy 0173 96h0E 89 
| uP oAWPAROTOG OH} =—«_ -, 754. 91 VOPK OL! 7% ry) of MeO iSE~ne = §=—. , 0400 399 292E O89 
2? fw AY IY CS) ); 25 pee ie eh fe 24 oF? sPAApeSsatoge «=. ,. 7200029 9E 9046 


nf, 


AT eae RR mS sre 


THERAT CONDITIONS 


Le PROOKOTIES 


‘ TR ERE OOO ANOS STINET Sa A Ee apm cei cme 
oi 


<TNETIC ExPaNslon 


PFOFORMANCE PAaRauETFas 


 vacuue THRUST COEFFICIENT 
varies SPECTEIC @PLSE ¢S€C) 
CHARACTERISTIC VELHETTY (FT/SEC) 


0 12363700€ 083 


wale SMAFE o*B4ORE2EE 008 
SBESQURE cHStad 2 37PARAN BAF 03 
vELVENTY CHV sere) 1964201 APF ONE 
AFNETTY (Lasers 0 FAQSUR VIF OOO 
TEMdG RAT OF (NERRO) oSA61 PRN F OOO 
ENTwALOY taTris_ad me 168001206003 
AAS SWECIR AD WHINY o HP ARACPATE CR? 
AFAT CAPACTTVIAT Us RoE) «=. . 0 FPR IS IE COO 
Gata 0} 10672678001 
UALECUL OR WEIGHT. UTES oe PRanaSCar eo? 


INTEGRATION Paha TERS 


gree SiZz€ 

axvTaL POSIT lOn 

prarent EntHaL Pv CHancr 
CEACENT “ASS FRACTING CHanée 


eueMsCay COMPNSTTION 


<orcles mace FRACTION wmF FRECTION ua, SPECIES mass Faact ion 
a eA®POPZ0l Fahd = =—«_ , BAAN TTPZOF OE ? arce 0 17301 39088) 
a1 Ci? oe LARGAESEF AF] —«_-. 6 3 91103 OF 20? ‘ ALCL 2 %6594620C 99 
aye oL2U66AGM elk =, BENGOTIAE SNS & ayo o 1AAQTOSNE ROS 
aLner eS TRIAIGOF CIP? 1960178992 ® atom oASAZ 251 OE ORS 
a, 02 oMPNPEMH Me «=, 1.301 0 3OAEW GO 19 a. Gan o 11786600Eon? 
41,20 OL PEBARAG CFE = =—«_ , P2E THROGE LNG 12 AL 2031S) A. 
& 20TtL) oPTIOVAIF SHG =, KDHE VIVE oO) 16 co 2271 7060SE Cog 
fa 3 ePIAWALTESA} «—«-_ , FUFVOL INES} 1A ti 0110e987SECn] 
¢.0 02070101 F oH =—«-_, 1 2HOHNZOF WGK ya 1.2 o2ThasrRetene 
“ eTCOKID2IG AD «=, PARSE PPOSESO! ae wey » }A9979OGE CRG 
wap eSSean7OMPn@A =. ,. | JRL AP PE aN 2? He o PIA PAL PAL CA} 
age eAON7POOEOR) «=—«_ , } FQKATORESAN re wd eth AP>TIESARS 
us eAIPIOARAF®:.) —, 79068 9OOF 0G} a6 0 o Le} r3R3ELCRg 
fs ePROREROIAG OAD «=—«-_ , EAR TPAOI IF O09 an oF 0 4h0e 2000 one 


5-100 


oSe0e0800E -62 


Oe 
©.1010208 3E 0) 
2 300 7ORGOE—16 


WOLE FRACTION 
0 40766609C 082 
01 9125796E°09 
9 6376921 9E-00 
o 3%750228C-03 
0011164 7ECO3 
0. 
0 24062300C 066 
» TOn8G000E 02 
010646 720Eo8s 
» 1 20S 2870C 000 
0209079008 008 
0169 9030Ee0) 
022471 600Lc03 
2 37600560 906 


E 
: 
; 
iz 
E 


« 
a. ~e Gee e- . @ eu 


ERPANSICN FANOTTIONS <TVETIC FrOeANSION 8 anfa RATIO 2.finea 


Sea SALTS IEEE 8 A ESTER LEP PFE PT: ROE TTA Se eas EET ED Cre 
es 


Fite PaQAFATT SS PFUFOWMANCE PARAME TERE 


wate NUIVRED 


AS ITELIEL ed | . VeCtnem THRUST COEFFICIENT 0 16GA07OZE 004 
PRECSURE (PSTai eo TAPNPATIFE ON2 VACUUM SPECIFIC 16M SE (SEC) o83602771 £009 
VFLACITY (FT s6EC) eh WAR2620F 00% 
ENC Ty iL. A/F T3) 2 2O7RGEPHF HO) 
TEVOFOAT UNE {NF Rol) eORRARPRAE ONO INTEGRATION PaRanE Teas {1 
Ful saLev Catuese_as oo 1 SPbaAAAF ONE | 
80S MAL ECIn AD WETAHT ePAIAIAI9E AZ Sree Size oS00edeen'e 09 ; 
AAT CaPanseTviBTUsLAeNEH) «=. & 380} FAVES OO ANTAL SOSTT ION o 398 26925E 004 | 
Raves 01 1eNeeacteel Pracrrny ENTHALPY CHanae © 0222736 34£-6} | 
Cumeee?e AnNeveos oSAHORI7)E 062 PFACENT “ASS FRACTING CHanee o0aSe21 E619 a 
MFC Oe WEISHT IME E) oPRRTAISIF COS i- 
rutmr ea, COMPOST TION ! 
0, swectes MASE FRACTION «=6—MMLE FRacTTON NN, SPECTES MaSS FRACTi¢n “OLE FRACTION i 
7 
t “tu oPOVUWIIF CHR «=—«s-_—-,. J PRANS BAF wd ? arey co LAP 7O2GPF ON} = =—«_- , OSNSNHRAT HOF i 
d t| 
a se an o DS TARORR AE =, LEN TANO DF AOS a = =ato cANOPGRORE CHA =—«_. 1 ZOGRARE 04 
3 v a MCL eSONOALTEG¢ 292 «=, 1 7A PASE TIF CRP a AL tn oASTAORIZECHD § «=—«_ , PUSS VAE~OZ 1 
i 
’ e are cOMANN TOR OR = -63N11090SE-as 1h ALOR 100725 TAECHD =, , 030801 BTEROD I 
at ayaa o F422IOGOF 0.38 PP eALREL oY 1? AL203¢S) he ee H 
3 OL PO3ZtL? oP PTOUAIE OGG §«=—«-_ , FNS NEP MEE WO) Ve co ePACUITTEECOR «=, PARTOTTIE COR : 
q 18 tat oPPUARTL IFO} = BANAT ANGE RO) 16 CL oSASeV920Eonp «=, 10D ERZOE—O2 | 
a? €.9 0201 1U9AAE HOG o17137CO2E 006 1&8 cL? 0519s Re Eons 1806S 320C°0S 
98 “ eOS7IPKIOE CHD «—«-_ HH TUS 2OGF oN 0 wcL olFTEOOAECNS =§=—-, 8B FRAHISOE OO | 
1 “ap eVORWAIRNE RAP «=—«_ LA NGQTEPF-OR >? “> oP? PERE ont =. FON PGCE 08 
92 noe eASTOMPSIF SNM) «=. 9.727792 9FonA ae no e4A799THNEONG § «=—«_ , 0 DAHORASL COS 
Pt ND oADPPORCOF..€* =—_ , PACK 29OF oD) 4 C oTASAlAD2Eong =. 26 7899S OF HOO 
v fa oARaPAQOFE oH9 o 1081 ASA cap rd) Oo? oStIALFIIEoNG =, 0b POR OTE 08 


0988 0000008 800009040 00000008! 6960000600000 08000 S000 08O8 0H 0O0O00 0 0O0 00000 00000600000 0000000000 0000000005 0000 


AEKIN SOL IATFICATION AF aL ereiL) tos 


Bo TAOTGI IPE OD OFLTING PFuve 


5-101 


©. 18 7NERROE CES 


eee ee ae a a a ran = eee ee oe 


® ; « 


CEPANSTON FHUNTTINNS BTNETIC EXOANSTON an€a ®arto 0.006 


FL fe Snrecet tes PFAFOVUANCE PARAuE TERS 


PYYTTT it dd 


FuD OF SMH TNTETCOTIAW OF AL ongte D 


6? 2s 


V.9PIUIIEMEChO ABTA BATT he 


5-102 


aaCu WINNER oPanTesser cet Vac THRUST CHEFEICTENT ol GLEPHEKE OO 
SAFERUME (9S10) oP PESESaE oN? > varius SPECIFIC reer Se 1SEC) 026000713 28°93 
VELOCITY (ET SEFC) o? Fens ° 
A sewer ty (ease tts 2 1 A@RRHONE OO? 
¥E4IG WAT VME (NFReR) ebta7AaANF COO InftEGaatine Patan€TEas 
TyTea_ov CAties_Ay wo tAQSR1 166004 
BOF HWE AR VETANT oF O79 46060F 002 gree SIZE o$8000080E-02 
mF at CaPaCe@ly (Ris Ho A) 0674716 26E 060 _ontaL ONSET ION . o 2 7SO1 19 9E 06h 
RAtius 01 PH16P4 76001 PEACENT ENTHALOY CHANGE 2967605 360E 6} 
CUseENT anaves -§$OQ9n790F 6 C2 PERCENT MASS FRACTION CHansé 0 187700098009 
ALEC TOG WE TANT IOTES 2 PROPAI RAC CO? 
Cantar CHPOKTT ION 
wh, wectes maee FRacrign ere FRaCT ION uh, SPECIES wags Faactinn “OLE FRACTION 
. O eLLGANQVSeehe =, BNF 34A2SF OHO 2 ALCL LONGESESEWAY —_ 0509280 7E 67 
? accu? eTRAPESIG aA] «—«._-, @NONNE TOE OOP r arc. eHPORTEETE OND «—_, 851 99031 E09 
§ ae S54 741 OF oO eSISTINOZOT 6 a ayo - 10661 OTE SRS 9 oSS108eet~es 
? ay ocr ASIAN RE «=—«s_,. BHA VS TEL ONE e aLow o3*. MOVtIEoas «=. B3H9AI VED 
e 0,02 eDOORTPTIFR&ES «=. 1:71 917 760K 984 te aLO2e eOPAIVSNBSEAD «=—«-_-. 4 99 9.998993 
Wi 079 ePRIATIEREF OFA =. 198] OAT Wee ? 6.203189 eLOOPaTOLECOR «=, 08 9BROOAE CS) 
a3 Oy 70R0u9 eM PAN6089F of} oP PPERSSOE <O1 ¥e €o 9 PATHE TARE CHA eBevvesreceoand 
a8 c%? 9 2ORCCLISE og B91 WHO RF OO) 16 <5 eISaPREPESHS , BAOTEHLEEOF 
1? eye oh2O7RA0GF. OF  ,3100971SC007 TL) c.2 e2tenerarstees «_,, PEN ZOZITE 06 
\ en. RT A ee a eTROFOAATE CRS 1 30OHEERE OE 
PT RTITS (i a) ts on oh | sPUMALEDESL «2000 75097008 
23 “90 MOOS IPIGFOG} «—«_. 1: 9 9 PSADIE CRO 6 no eQOPiI7POERD §«—_ . 0 20002 0S E989 
nS eABPACNEEF CH} «—«_ , FTIANVZNTE WAI |) a) -THOTSEPPESng =, } OT ASHEE 96 
2? fe etAGARPENF OOS = =—«_ , SATAPPSSF CON 2a 02 VS1a20VIoNE «, 2H00 907 1E 08 


PPPPTTT TT TTT TTT TTT et alate ate 


@ SeARPOPECOO 


= . 
: a: 
fod 
5 CRPANG INN FAUNTTIONS KINETIC FrOANSTON amEa RATIO 6.0n9 
® 
: F.0e PRH*FATIES _ SFRFORMANCE PaRawE TERS 
° wate WeRte oPTPPIANAF CAL Valine THmysT COEFFICIENT o $A80339SE 009 
° eRecKuer (eStal 0 9 5A990 296 002 Vacuum SOECTEIC $M SE (SEC) o27243280E +03 
VELICITY Ce’ 7480) o PARITA WEE CHS 
NENGITY (Lasers otha redragebl 
TEMRENAT IRAE CNEROD) 2 BHAISARE SHO InvVeEGration PaRame tras 
Fulualey CETERA) oo tO?) YKOF0gs 
BAS MALECIR AG WET RHT o PASAT AAC cE? evew Size o SROROOOOE SOP 
aFAy CaracrTy c4tis, Meh es + 90027146000 avfat POSTION o ROFEOATRE 004 
aae~a 0} 7nad7 AE eH) SFREENT ENTHALPY CHaNGs @,52967073SEe0} 
CUMOENT Bheveg o SAAD PEATE C2 OFACENT “ASS FRACTING CHahGe o 2674 7ROOE o6O 
UU CCHAR GETIMT CER o FRO VVORRE COZ 
Puret sa, COMPOSITION 
wo, serctes MASE FRACTION = =6—«MILE FRACTION Ww, SPECIES wass Faacting <O.€ Faactiow 
’ a BISTROT LNG =—«_, MAN HN PF anS ? arc. ol P2VOO{SEony =—«._. ON KOEEE MHZ 
3 acue o LPUARTEPS. 9} eAeAnrl SUF a0? r acer othCaretEmag .21%60305E°03 
3 aye ePONGWUEPFeAB «—«-_ , PRO OE 78TH 96 4 ayo eTAQAATTITECOR§ , 0THSOREEE COS 
? ay Cn eATAOCUGI AAP «, P7P3AIZ IEF OND a a On oD7MAAOTEECHD § «=—«.- . 21 95 1 $01 Eo 89 
e aod ePSPOIRURG AAG =. 2: 91 0 TFORE NE yn a a] oFAPeOPATEoNG =, 03 060 90SE ORD 
18 a 20 ePOSSIADD FOR «=, , OMA PLEO TF ANG ” OL 2090S oPPTODIATECNG «=, PERV OOAEE HO) 
19 A. 203¢L) fe 4, 16 te) e2O™VOZZiEeag =. 2075040 9E 908 
8 Cre 03001 POFS eA, «=—s_ , 1 79D DeFeet 16 CL o"PSazeTekuge 270669508 002 
a? Cao o2OPHPLEH AP = L BPPEPORE WOT 18 cy? TAPLPOTEens «=—, 00 PAB TEL 86 
7e a tPA TPORGoA? «=, B4ALI1OGF OP a onal 5 el A97OOQIEone §=«, 8 349 08REE SOO 
FY) “ay el OQnPROPG.4Q =, SH TTEREF oA >? bel d »Zd7e10S2Eong = =—_ , FRR OREE CEO 
22. ue eARTOVMPEF EHD «=... 2 1 ONT E98 » wo oA ATE DVPELCgS «=, 62190 9 00S E OD 
re ay eAWARCOEF eh} =—_ , *TANAOP AF oD} 6 0 eSB aMyareeng 804 PARROTS 
a? fe ePROGARNTF ACH = , wa RADIVF oO? aa ey J eHOOnQeaREons =, 199URTEPERSS 


5-103 


<n cctus - ue << exesmuaasimneanont_catguumaanatn. ct anes na ain oa s.na bhp adele. sec sn maet Sistas itn obi sith. vonmaamabhielay navidsldiddi, Mihi ede tell 


FLde Surerelies 


24Cu WICAFR 

OMECQ ee (esta: 

wFLWOETy (elses ry 

ATNCITY tasty ys 

VEUSE RATION (16 Red) 
Cates, ad 

SO6 GN CON 22 WE LRNT 

at CaPacrly ATU Rong) 


TutusLeov 


panes 


CUMOENT BOeVeR 
CAFC AR wT GeTiers) 


wh, SPRCIES 
»  & 
e accup 
© Mr 
v £4, Cr, 
© are 
i 6,20 
4% M& 2000, 
is cw 
i? ea 
36 ro] 
a J 
pr es] 
~s © 
Te) 


FEPANRTOe CHAATTIONS IVETIC ExPANSTON ana Ratio 


0 WPMAWIIF COL 
o Sha SUNORE oO} 
ohe tPbO32F eas 
oAPATS 1 PAHO? 


oR4ROTOE1F 0046 
=o 21 SAPROAF COS 
o PAUSIANW 0 O2 


oSPa7earfee? 


wate FRACTION ame FRACTION 


oF 0A28C0. 06 
Oe hh dee | 
oOCTRI21S¢—17 
080847) 007 —02 
ot AARM) 29 ONG 
e8 DTT Per ogo 


039701 POF on) 
2482679000 208 
oFOSenzeas 2% 
odtavenyar ice 
ESS Ch Pe | 
oA WPCORCOF og! 
0 160700) 0083 


2 PRO PCear oO? 


e 


16.008 


OFAFHAUMANCE PaRauEt ras 
VACIN® THRUST COEFF IC TENT 


" vat SPECIFIC Ieeinse (See) 


InTeGMat ion Patan€tges 
Sree S$tTZE 


avtaL SOSTT ION 


OFRCENT ENTHaL PY Cranee 
PFECENT “OSS FRaCTING Crane 


PuFeiCa, COMPOSITION 


oSMPRP TS 2F 306 
o389n) 0667 one 
0 PODIOF LF 
0 PAAD 7O4ConP 
oAPnete d7F 006 
oAAe 00) TOF 60S 


018001299001 
oPeDVIAl Peoon 
2 7001 7S0RF oO? 
eA PASADOIF COO 
o8Pta7260C eee 
0 PING PO PF oO} 
085226601550) 


WA, 
? 
e 


SPECIES 
ace 
arcu3 


aL On 
8.27931$) 
co 

c. 


5-104 


MASS Faact inn 
0} PZAONP Eons 
0 7784095 7E oes 
2bb$ne Tart oan, 
o 39209391 OL ons 
0067) 2590E ong 
2772330 7Eond 
0 Fara 3200" 00g 
o ETRE RECAP 
o 317A8660E coe 
0 1% 7 198eg oan 
0390686 Mons 
oP) 79S ED 
© 990631942009 
ob} 3eS1 2eCoag 


o 3 FE6RS9OKE CO) 
o&RS21007E9O9 


COLE FRACTION 
040060038C 62 
80s ?29SrEoed 
2 2600G9CAL<05 
0} 9978200C-03 
0438949928093 
« 70700990C°0) 
02070 70062 206 
0847000 36C008 
od Besgdzegors 
SR 
o 306) 68068008 
0421070120009 
oS SPL O7 7ECoRG 
9 99 090030C~06 


2 + 

| 4 

4 = 
Reproduced from oO : 

| cone available copy. 

yi 

: 

' : 

. ‘ 3 

: c ERPANC TAN FAUNITISNS <INETIC Ex®ANSION anEa BaTIO 20-800 

FLMs SNgPERTT =< PFETORUANCE Fa@auE EAS 

4 wats WueaeD 0 A860 7F 003 — varuue THRUST CHEFE TC TENT o 1aSne 202k oe) 

PAFSSIOF (OTs) 0 WH9GO34F oO} vacuum SPECIFIC T4maSF (SEC) 02009153 3€009 

VELXCITY (FT/SEC) o OR INGAAAE OH 

YENCITY (Laser ePenTrAgnrea2 

VEUDE IAT IME IW Ao8s e287 IL OE2F ORO TwTFGRAs 10y Pahang teas 

Fut <oLev Cay tsa we P IAP PONE 

BAS OY FN'H AR UF I RHT o P0184 PRAK ON? SvFP SIZE 2 3e000RCOE OP 

wfar CaPacrlviMfuAenFG) 009809985 008 antec POSTTiNN o8O2Q6300€ 64 

Rates PFILSIEY Lal OFACENT ENTHal Oy Cranes @e 

Cueetv? Bneves oSO 70700 78 002 OFSCENT MOSS FRACTION CHanee o6e02TPO2E 6G 

CL FCN OO WERENT COLE) 92?Se300RRF 062 


PurMESs, COMPOST TION 
em, seecirs MOSS FOACTI ON er_F FRacrinn wa, SPECIES wats reaction woLe pRactioN 
 ] 1. 0 FORT TAD Or 0% oP WA PORT PE 08S 2 ares otdest Oregoes 09376200 3E oO? 


9 o.cLa oSOAROG1SFoAL =, PH 7CKAESE ROP Y acces oI407E0STECRD «=, FOLOP LATE SO 
. 8 eeterarz~.ny =, ORaaTEsFoan & ao olTPaieSeenes =, 107313 00EoES 
y 00 ohOPPRISRF OO} §=—-, PVVTIEZAF OAD 9 aon sPeReCeerecen 1 Peandereced | 
® ang POV IASG ROK «= _, BITTEN OF OLA ie o.O2w Ree eT eee 203 = 
19 #LR0 eDVPIANGTSOR «= T17G7SS7E00G = 1B AL PNTESD o2770D10Tong = PoveNRECESaS 
al aenmas 2 ye 0 c2anpITTERoan = ,2eaSeTaeeeee 


1% cae 9 TPP 70CR PF AA} 2 MOONS EF oH} te Ci WIAs ISeens BDA PeseRTeeR 
‘? oe POR IYS TS oo REL TAL ee | 1e cue S61 PGCE ony =, 141 OPTSEE COP 
= Sa aa a I) a UO MPLS 
t 2 | aap oDPUARLOF OL @ PSU NOO OF 206 ry} “a .7 ase ieee) .00neerences 
= = ePIBNNTNAGOM, = tIAOISPECOA = fe NA ooSTSFAM CMD «= 02ND AOR ETSED 
a ue oAOP°OF 0G 08} aTTOevyeowoaAt ry 0326200108008 ereraricees 
29 Mes eVUAPEARLF ee = =—«_ , DTP PESTLE LAY pa 0? Any rerzegonn — , aganenneroes 


i] 


Se eer SS a Pere ee 


FxPANRTOn FOMRTTIONS 


Fite PnrogRti«s 


e 


CEVETIC FuPanszOn AREa Patth 


PFOFORGANCE Pahavt Tres 


30.000 


uote Quenge 2 Waneersr o6i wecdus Venus CHEFF IT. rent Pe | 
ORESRUPE (ORTOD 9 tanererer cet vacuum SPECIFIC 1eP.08E SEC) 0 306286 30€ 063 
vELcite eet setes et@dVAAGR iE 204 
rewesty ¢ashydy oe f@anMai Foe? 
TEaSE RET INF (HERS) e 10 90F 0 04 STF GRation Pahan€ tres 
Cet MALPy CRTs an we Pet ‘ 
308 COR AA UP Tan? PRP waren? oOF gver Siz€ o FSh50GGK OOP 
aay Cavacrtviatus.sowei , pea gener aula POSE ION o 1840661 30009 
sates mToi ch 2 Brertut Entra. ev Chaves =. $269087 7E-8t 
PUBTE NT OReves 2 pheg oO? PFOTENT “ASS FRACTING CHaNeE oSt?ETO3IECOD 
WEFT NAR WE TENT ess) 2 Phe Pegnas c 62 
euesicn, cOuPestT Il 
SP ated LL | maee FEACTICN ance FRACTION wh, SPECIES wags Faactinn “Wil FRACTION 
| a. ot090ee aT e8G =, 10801 1S OF ONG arc. 089CQRCHEE—y = _ L ST PA PIAECOZ 
2 CUP eWERSUGO1F SP 8 =—«_,. 997 LAN BOEugP auc.s e8BP2TOMeRD «—_ , 00H FF DOTH 046 
$ me QUaPePEEa™™ ,FEK0 NEW aA se Re 
9 aPC colt PQNTAsF oy = =—_ , SPD OAEORT OF e a0 o27h08141Eon9 «=. , BOP ONNSETH OT 
& a CP 0900000 78 GH 0-261 0ORIGEE COA te dL Oem eVeRaT7EFoey =. 0 479 BP. OLS 
ut 6.79 olGPEPTNCF HG =. 69.9.1 19 290 COS TJ a4, 203089 e2?PedierZene§ = §«,, 9079 79S0E 08 
B 64, PORT @. @. ve ce POOL OSOPECRG «©, 806098 90E 00 
8 cap o BPO 17 806F 28} 91098683 IE est 6 +5 oB8Sripeekene 8 «—_ , 840044 0C 02 
7 6A RY SES Me so) en) ee: © sP9veaeLony 08005 67 L080 
30 208 | = e804 Seong §«, 1949 00688908 
t A) wap Alita te te. eS ba] a “a « 2IHGRINGT CRY 9 907591498000 
a A i ee |. a) cA? S72 GS «=. O2NEEAOPE CSD 
aa od o8BPOCQCEF eH} =—_-. VPFOON EPC eSh a ° mIYY, i) ¢ MIT Ss. yc. 
9? Owe etPariengtets «_ Hao FH FOTIC 4 a) cf . yapade WMieeq 8 ,F4029 7468 8d 
FO wat af ACutA 


5-106 


ashe 
7, anene 
Se ANACE 
Ae AHHOS 
18, neNnEn 
Fa, ANKEA 
30, @neee 
30.8800@ 


GUeuaty ATRUT Fane 1Oq WNNLE 
1S? (Amery 


PtH. 82RH0 
Pah, 6aTaa 
PTF 24393990 
ue. Plane 
290,33500@ 
904, 27000 
30h. 2Pang 


EVACRIND 


5-107 


6. References 


Gordon, &, and McBride, B.J., “Computer Program for Calculation of Com- 
plex Equilibrium Compositions, Rocket Performance, Incident and Reflected 
Shocks, ani Chapman-Jouguet Detonations, " NASA SP 273 (1971). 


Nickerson, G. R., Coats, D. E., and Bartz, J.L., “The Two-Dimensional 
Kinetic (TDK) Reference Computer Program," Enginzering and Programming 

Manual, Ultrasystems, Inc., Dec. 1973, prepared for Contract No, NAS9- 
12652, NASA JSC. 


Nickerson, G. R., and Kliegel, J.R., "“Axisymmetric Two-Phase Perfect 
Gas Performance Programs, Volume 1 Engineering ard Programming Descrip- 
tion" prepared for NASA MSC under Contract NAS9-4358, April 1967. 


Weingold, H. D., “The ICRPG Turbulent Boundary Layer (TBL) Reference 
Program," p-‘epared for ICRPG perf. Std. Working Group, July 1968. 


Barron, J. G., Jr., Cook, K. S., and Johnson, W. C., “Grain Design 
and Internal Ballistics Evaluation Program, Program No. 64101," Hercules 
Powd=* Co., Bachus Works, July 1967, AD 818321. 


Harry, D. P., Price, C. F., Small, kK. R., and Tavior, D. E., “User's 
Manual For Nozzleless Rocket Motor Interna! Ballistics Computer Program, “ 
AFRPL TR-73-20, March 1973, 


McBride, B. J., and Gordon, S., "Fortran IV Program for Calculation of 
Thermod ynam'c Data” NASA TN-D-4097, August, 1967, 


Stull, D, R., and Prophet, H., "JANAF Thermochemical] Tables, Second 
Edition," National Bureau of Standards NSRDS-NBS 37, June 1971, 


Pieper, J. L., “IC8PG Liquid Propellant Thrust Chamber Performance 
Evaluation Manua.," CPIA No. 178, Septembe. i962. 


Thiokol Chemical Corm., Elkton Division, Elkton, Maryland "TE-M-364-4 
Extended Delta Rocket Motor" Final Report NASA Contract NAS7-678, 1972. 


6-1 


Appendix A - Control Card Set Ups 


Control care set ups for two typical cases are shown for both CDC 6000 
(SC@PE 3.3) and Univac 1100 (EXEC 8) seires computers and operating systems. 
In both cases absclute versions of the SPP code and the thermodynamic data files 
are assumed to be on permanent mass storage data files. 


The first case considered is when one or more of the computational modules 
are executed and the :inkage data is written (punched) out for subsequent runs. 


The second case considered is when linkage data from the first case is 
used in executing subsequent modules for the same motor. 


The bracket symbol (}) is used to explain or comment on the function of 
the control cards. Since the first, and most Itkely the second, control card of 
each set varies from installation to installation even for the same computer and 
operating system. The accounting information control card is left off from each 
set and the job specification contro! card (normally the ¢econd control card in 
the set and the first shown here) should be checked for consistancy at each in- 
stallation, 


CASE I: CDC 6000 Series Computers 
J@B, CM10006, CLiz6000, T700*, 1D1000. 
ATTACH, SPP, SPPABS. | get SPP absolute 
ATTACH, TAPE 25, SPPTHERMG equaie thera:c “ate to logical unit 25 
RFL, 126000. |request core to execute SPP program 


SPP. execute the SPP program 
7/8/9 

program date 
6/7/8/9 ond of file 


CASE I: Univac 1100 Series Comput. 
@ASG,AX SPP. | assign program iiie 
@ASG, AX SPPTHFRMG.~ | ass. thermo date 
@USE 25., SPPTHERM@.* | equate logical unit 25 to thermo data 
@ASG,U_ LINE XATA,,F | assign a file for the linkage data 
@USE 8., LINKDATA, | equate logical unit 8 to the linkage Jata file 


@XQT SPP.SPP execute the absolute element SPP from the SPP file 
program cata 
@FIN 
@@ remote joh entry end of file 
CASE II: CDC 6000 Series Computers 
J@B, CM10000, CL126000, T700*, 191000. 


ili REA FMEA TORT OS CUI SITEMAP FOL a PSI IMAP NC er EIU RT ESLE AAI ATTN 


ATTACH , SPP, SPPABS, 
ATTACH , TAPE::5 ,SPPTHERMG@. 
C@PYBR, INPUT, TAPE3. copies previously generated linkage data to TAPE3. 
REWIND (TAPE3) 
RFL , 126000. 
SPP. 
7/8/9 
linkage data from Case I 
7/8/9 
program dita 
6/7/8/9 
CASE II: Vnivac 1100 Series Computers 
@ASG , AX SPP, 
@ASG, AX SPPTHERM@, * 
@USE 25, ,SPPTHERMG@. * 
@ASG, AX LINKDATA .** 
@USE 3.,LINKDATA,** 
@ASG,U LINKDATAZ2, ,F assign and equate further linkage 
@USE Pp. ,LINKDATA2, | até to unit $ 
@XQT SPP .SPP 
program data 
@FIN 
@€@ 
NOTES: 


* CPU time lim:t depends on modules to be executed. 
+ Only needed if the @DE and/or @DK modules are to be executed. 
** If this data has been punched, the following control cards can be used: 


@ASG,T 3.,F4 
@DATA,I 3. 

linkage data from Case I 
@END 


